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6.0 KU-BAND COMMUNICATION SYSTEM ANALYSIS 

The studies and investigations of the Ku-band communication system 
were focussed in the following areas: 

(1) Shuttl e/TDRSS si delobe avoidance 

(2) Shuttle/TDRSS acquisition and tracking performance 

(3) Ku-band TDRSS ground station, two-channel, Costas-loop tracking 
performance for unbalanced quadriphase-shift-keying (UQPSK) with arbitrary data 
formats 

(4) Ku-band bent-pipe mode 2 BER performance 

(5) Degradation due to dither of the antenna during angle tracking 

(6) Potential Ku-band sidelobe reduction. 

Analyses of Shuttle/TDRSS sidelobe avoidance as well as Shuttle/TDRSS 
acquisition and tracking performance are presented in sections 6.1 and 6.2, re- 
spectively. Analysis of the two-channel Costas-loop tracking with UQPSK sig- 
nals with arbitrary data formats is presented in Axiomatix Report No. R8010-8, 
dated October 30, 1980, which is included in this section. The analysis m 
this report was kept general, so the results can apply to a variety of Shuttle 
payload data formats. Specifically, the mean-square, phase-jitter performance 
is determined for the TDRSS ground station, two-channel Costas loop with active 
arm filters, synchronous or asynchronous symbol clocks, and arbitrary data 
formats. 

The analysis of the Ku-band return-link mode 2 BER degradation is 
presented in Axiomatix Report No. R8102-4, "Ku-Band Bent-Pipe Mode 2 BER Per- 
formance Degradation," dated February 18, 1981. There are three inputs to the 
FM modulator in mode 2. First, on channel 3, there is the bent-pipe wideband 
signal which is amplitude (peak) regulated to provide 22 MHz peak-to-peak (p/p) 
carrier deviation. In addition, channel 1 has operations data from the network 
signal processor, and channel 2 receives data from the payload recorder, opera- 
tions recorder, payload low-data-rate link or the payload interrogator narrow- 
band bent-pipe signal. Channels 1 and 2 together provide for a p/p FM frequency 
deviation of 12 MHz. The specific case of concern in this report is the one 
wherein the signal of channel 1 is a 192-kbps biphase-L (Manchester) baseband 
data stream without thermal -noise corruption, channnel 2 has an NRZ baseband 
data stream plus thermal noise, and channel 3 is composed of two subcarriers, 
one at 1.024 MHz and one at 1.7 MHz, plus thermal noise. 

PRECEDING PAGE BLANK NOT FILMED 
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Table 6.1. Si del obe-Discrimi nation Budget [1] 



Weakest 

Mainlobe 

Strongest 
Si delobe 

Cumulative 

Difference 

Pointing Loss, dB 

-3.0 

-22.0 

19.0 dB 

Flux Density, dBW/rn^ 

-126.9 

-113.5 


Antenna Gain, dB 

37.1 

37.1 


Effective Area (omm), 
dBm 2 (at 13.775 GHz) 

-44.2 

-44.2 


2 

Effective Area, dBm 

-7.1 

-7.1 


Polarization Loss, dB 

-0.2 

0.0 


Received Power, dBW 

-137.2 

-142.6 

5.4 dB 

System Temperature, 
dB°K (at antenna 
terminals) 

30.6 

28.9 


Boltzmann's Constant 

.... 

-228.6 

-228.6 


C/N 0 dB-Hz (at EA-1 
input) 

60.8 

57.1 

3.7 dB 

Threshold Setting, dB 

-0.3 

0.3 


Noise Immunity Margin, 
dB 

-0.6 

1.2 


Effective C/Ng, dB-Hz 
(at Detector ) u 

59.9 

58.6 

1.3 dB 
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Another aspect of the Ku-band investigations performed was determining 
the degradation in the Ku-band communication system performance due to the in- 
troduction of dither during antenna angle tracking. The dither is introduced 
to overcome the effects of stiction and the resultant limit cycling of the an- 
tenna during angle tracking. The analysis for this performance degradation is 
presented in Axiomatix Technical Memorandum No. M8105-1, "Degradation Due To 
Dither," dated May 11, 1981. It is shown that the dither has no significant 
degradation on either angle-tracking or communication performance. However, a 
minor improvement in steady-state angle-tracking accuracy is required to main- 
tain the 0.3-dB allowable pointing error with dither present. 

Finally, Axiomatix Technical Memorandum No. M8302-1, "Potential Ku- 
Band Sidelobe Reduction Recommendation," dated February 16, 1983, presents an 
approach to modify the Ku-band antenna-feed-support structure in order to re- 
duce the sidelobes. The relatively high sidelobes found in radar and communi- 
cation systems have caused problems with acquisition throughout the development 
stage. The high (~ 20 dB) sidelobes in the communication system continue to be 
of concern during TDRSS acquisition. From the antenna measurements, the domi- 
nant cause of the high sidelobes appears to be due to antenna-aperture blockage 
by the feed-support structure. By using a wedge-shaped ridge along this struc- 
ture, polarizing characteristics can be achieved. The linear polarization, con- 
sisting of both RHCP and LHCP components, will then contribute some offsetting 
aperture illumination versus none with the existing configuration. Axiomatix 
feels that this modification can achieve sidelobes of at least -22 dB. 

6.1 Shuttl e/TDRSS Sidelobe Avoidance 

6.1.1 Introduction 

Some concern has been expressed about the possibility of sidelobe ac- 
quisition of a TDRS signal, particularly in light of revised flux-density speci- 
fications at the Orbiter. Note that Table 6.1 [1] gives the minimum flux density 
as -126.9 dBW/M 2 and the maximum flux density is -113.5 dBW/M 2 . The resultant 
margin is given as 1.3 dB. Several assumptions were made [1] to arrive at this 
margin, and the specified dynamic range of received flux density will be slight- 
ly increased in the revised specification. In the following sections, we devel- 
op a new pair of sidelobe-discnmmation budgets. One is based on the new re- 
vised flux-density specification which includes a margin at the low level, and 
the other is based on actual expected values of flux density. In both cases, 
a more realistic assumption is made concerning polarization loss. 
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6.1.2 Derivation of Sidelobe Discrimination Budgets 

The following derivation of the sidelobe-discnmination budgets fol- 
lows the format of Table 6.1 [1] so that the results may be compared. This 
includes the threshold settings and noise-immunity margins. In Table 6.2 [2], 
we present a conservative estimate of the sidelobe-discrimination margin based 

p 

on a minimum received flux density of -126.0 dBW/M and a maximum of -110.5 dBW/ 
2 

M . The minimum flux includes a margin of 2.5 dB added by Rockwell to the esti- 
mated -123.5 dBW/M 2 , based on the Orbiter ICD [2] which includes polarization 
loss. Table 6.3 [2, Table 4-IX] is reproduced herein and the flux-density val- 
ues are calculated based on a minimum loss of -162.0 dB M and a worst-case 

O 

loss of -613.5 dB M . In addition, polarization loss on a sidelobe should be 
greater than on a mainlobe: typically, polarization loss increases monotomc- 

ally off boresight. For our conservative estimate. Table 6.2 [2], we assume 
that the sidelobe and mainlobe polarization losses are equal. The net margin 
is -0.6 dB which, as we have stated, includes a 2.5-dB margin added by Rockwell 
to the flux-density specification. 

Table 6.4 [2] presents a more liberal estimate of the sidelobe dis- 
crimination budget without a margin added to the expected flux density. In 
this budget, we assume a polarization loss at the first sidelobe down 1 dB 
from the mainlobe, which we feel is realistic. The result is a positive mar- 
gin of 2.9 dB. 

6.1.3 Conclusions 

Based on realistic estimates of received flux density and polariza- 
tion loss, sidelobe acquisition should not be a problem. Moreover, the maximum 

p 

flux density of -110.5 dBW/M is predicated on the TORS pointing directly at 
the Orbiter. The TDRS Ku-band antenna full beamwidth is 0.28° and is open-loop 
pointed at the Orbiter during acquisition. We feel that the most likely event 
which would cause the maximum flux during acquisition is an attempted reacqui- 
sition immediately following loss of track while the TDRS is still accurately 
pointed toward the Orbiter. This should not be a problem but, in order to more 
accurately quantify the performance, the axial ratio at the most significant 
si delobes should be measured. 



Table 6.2, Sidelobe-Discrimination Budget [2] 


> 

Weakest 

Mamlobe 

Strongest 
Si delobe 

Cumulative 

Difference 

Pointing Loss, dB 

-3.0 

-22.0 

19.0 dB 

Flux Density, dBW/m 2 

-126.0 

-110.5 


Antenna Gain, dB 

37.1 

37.1 


Effective Area (omni), dBm 2 
(at 13.775 GHz) 

-44.2 

-44.2 


Effective Area, dBm 2 

-7.1 

-7.1 


Polarization Loss, dB 

0.0 

0.0 


Received Power, dBW 

-136.1 

-139.6 

35 dB 

System Temperature, dB°K 
(at Antenna Terminals) 

30.6 

28.9 


Boltzmann's Constant 

-228.6 

-228.6 


C/N 0 dB-Hz (at EA-1 Input) 

61.9 

60.1 

1.8 dB 

Threshold Setting, dB 

-0.3 

0.3 


Noise Immunity Margin, dB 

-0.6 

1.2 

< 

Effective C/Nq 5 dB-Hz 
(at Detector) 

61.0 

61.6 

-0.6 dB 
L 



Table 6.3 TDRSS-To-SSO Ku-Band Forward-Link Mode 1 Interface Characteristics [2] 


Information 

Time-Division- 

-Multiplexing 

Spectrum 

TDRSS Transmission 

SSO 

SSO 



Operational 


Spreading 

Carrier 

Carrier 


Polari- 

Min. 

Required 

Channel 

Rate 

TDM c 

216 kbps TDM 

d 

Modulation 

Frequency 

EIRP 

zation 

G/T 

P^„ / N n 
rec 0 

Voice 1 

32 kbps 


Bit Rate: 

PM Sequence: 



Acquisition: f 



64.0 


Rate: 72 kbps 

216 kbps 

Gold Code 



40.0 - 49.5 dBW 



dB-Hz 

Voice 2 

32 kbps 

± 0.001% 










Sync Pattern: 


Chip Rate. 

PSK 

13.775 


RCP 

6.0 

i 

Command 

6.4 kbps a 

76571440g 

Format: 

(31/141024) 
x Cxr Freq 

8 = -nr/2 rad 

GHz 

Tracking: g 

Axial 

dB/K 


Sync 

1.6 kbps^ 

Biphase-L 
Sync Pattern: 

±3° 

e 

Normal power: 
46.5 - 49.5 dBW 

ratio 
<1 dB 

h 

65.4 

dB-Hz 


Format: 









Forward 

Data 



26772706354g 

NRZ-L 



High power: 

g 



128 kbps 

N/A 

Period: 



48.5-51.5 dBW 



J 

Sync 

16 kbps 

N/A 


1023 chips 








Legend: a 

b 
c 
d 
e 
f 

g 

h 

i 

j 


Encoded and permuted command data 

1.2 kbps synchronization pattern plus 0.4 kbps arbitrary fill bits 

Same as S-band high data rate uplink and forward link operational TDM; encrypted for transmission 
Unspread capability exists 

Doppler compensated in TDRSS for 17.5 kHz at SSO 

Effective EIRP toward SSO during acquisition with open-loop TDRSS antenna pointing and allowing for SSO 
of 6 dB; not to exceed upper limit 
With Ku-band return link acquired 

Based on SSO antenna temperature of 10 K(no Earth or Sun intercept) 

SSO link acquisition threshold 

For BER of 10"^, subject to specifications in 4. 4. 2.1 




Table 6.4. Si del obe-Discrimi nation Budget [2] 



Weakest 

Mainlobe 

Strongest 
Si delobe 

Cumulative 

Difference 

Pointing Loss, dB 

-3.0 

-22.0 

19.0 dB 

Flux Density, dBW/m 2 

-123.5 

-110.5 


Antenna Gain, dB 

37.1 

37.1 


O 

Effective Area (omm), dBm 
(at 13.775 GHz) 

-44.2 

-44.2 


Effective Area, dBm 2 

-7.1 

-7.1 


Polarization Loss, dB 

0.0 

-1.0 


Received Power, dBW 

-133.6 

-140.6 

7.0 dB 

System Temperature, dB°K 
(at Antenna Terminals) 

30.6 

28.9 


Boltzmann’s Constant 

-228.6 

-228.6 


C/N 0 dB-Hz (at EA-1 Input) 

64.4 

59.1 

5.3 dB 

Threshold Setting, dB 

-0.3 

0.3 


Noise Irrenunity Margin, dB 

-0.6 

1.2 


Effective C/Ng, dB-Hz 
(at Detector) 

63.5 

60.6 

2.9 dB 
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6.2 Shuttle/TDRSS Acquisition and Tracking Performance 

6.2.1 Introduction 

The current scenario for Shuttle/TDRSS acquisition consists of having 
the TDRS open-loop point toward the Shuttle and radiate a minimum of 40.0 dBW 
in the direction of the Shuttle. TDRS pointing-loss considerations are dis- 
cussed in [3]. The Shuttle Ku-band high-gain antenna then executes a spiral 
scan to acquire the forward link. Once the Shuttle Ku-band system detects the 
forward link and re-radiates on the return link, the TDRS can initiate closed- 
loop tracking. When the TDRS is accurately closed-loop tracking, the minimum 
received flux at the Shuttle is increased due to the more accurate TDRS point- 
ing and the Shuttle, in turn, can maintain more accurate pointing toward the 
TDRS. 

Hughes Aircraft Company (HAC) had recently expressed concern that 
the Ku-band system cannot accommodate the wide variation in expected signal 
strength (dynamic range) and still maintain accurate pointing. This concern 
stems from observations that the variation in received flux density at the 
Orbiter, combined with umt-to-umt variations, thermal effects and detector 
weak-signal suppression, add up to provide the communication- track servo with 
an error signal having an estimated 39-dB dynamic range. The incident flux 
density variation was taken from the Rockwell specification. Rev. B. 

6.2.2 Dynamic-Rate Limitations 

Hardware limitations affecting Ku-band performance are discussed in 
Section 3 of this report. Briefly, a combination of manufacturing variability, 
AGC performance and temperature effects combine with the input-signal flux var- 
iability to produce a 39-dB dynamic range at the angle-track servo input. The 
servo was not designed to accommodate this wide dynamic range and, in fact, the 
capability of the servo is unknown at the present time. 

HAC system engineering has partitioned the problem between the servo 
and the detector: Investigate whether or not the servo can handle a 15-dB dy- 

namic range and explore detector design changes and specification-relief options 
which can reduce the dynamic range from 39 dB down to 15 dB or less. Hughes has 
compiled a list of 10 candidate options to improve the detector performance, 
one option is simply 11 dB of incident flux-density specification relief with 
a slightly narrower track IF bandpass filter. The flux-density relief is need- 
ed at the low end, e.g., the minimum signal level for accurate tracking would 
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have to be boosted up to 11 dB. While this seems like a prodigious amount, an 
examination of the Rockwell specification indicates that it is possible. This 
is discussed in the next section. 

6.2.3 Flux-Density Specification Considerations 

HAC dynamic- range calculations using received flux density are based 

on the Rockwell Rev. B specification. Apparently, Table 6.5, which appears on 

page 29 of the specification as Table 3-1, was used to determine minimum flux. 

This table is reproduced herein and gives "minimum EIRP required for autotrack" 

as 35 dBW at the TDRS. Using a worst-case space loss of -163.5 dB/M and a po- 

2 

larization loss of 0.5 dB, we arrive at a minimum flux density of -129 dBW/M . 
This is the worst-case flux density used by HAC as a requirement for accurate 
tracking. 


Table 6.5. TDRS RF Interface Requirements [1] 


Parameter 

Return Link 

Forward Link 

Center frequency 

15.0034 GHz (nominal) 

13.775 GHz ± 10 kHz 

Bandwidth (3 dB) 

225 MHz (transponder 
throughput) 

50 MHz (transponder 
throughput) 

Antenna gain (55% eff.) 

52.6 dB 

52 dB 

Antenna polarization 

RHC 

RHC 

EIRP 

N/A 

48.8 dBW end-of-life for 
data transfer 

40.0 dBW minimum for de- 
tection in the direction 
of the Orbiter 

Maximum-allowable EIRP 

60 dBW at Orbiter 

55 dBW at TDRS 

Minimum EIRP required for 
autotrack 

30 dBW at Orbiter 

35 dBW at TDRS 
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O 

A flux density of -129.0 dBW/M is overly pessimistic since the 
Shuttle ICD [4] guarantees a mininum of 40.0 dBW toward the Orbiter, including 
polarization loss. This is equivalent to a flux density of -123.5 dBW. More- 
over, this flux represents open-loop TORS pointing. When the TDRS is illumi- 
nated by the Orbiter return-link signal, the flux density will increase to a 

o 

minimum of 46.5 dBW for an equivalent flux density of -117.0 dBW/M . The maxi- 
mum TDRS EIRP toward the Orbiter is 51.5 dBW. With a minimum space loss of 
-162 dBW/M 2 , the maximum flux density is -110.5 dBW/M 2 . The EIRP values come 
from the Shuttle ICD. From Table 6.3, note that, when the TDRS is closed-loop 
pointing, the flux density at the Orbiter is 12 dB higher than the worst-case 
flux used by HAC in the dynamic-range calculation. The maximum flux density of 
-110.5 dBW/M 2 is 3.7 dB less than the -106.8 dBW/M 2 employed by HAC m their 
dynamic-range calculation. 

6.2.4 Summary and Conclusions 

The revised flux-density dynamic range provided by using the EIRP 
values from the ICD, which are more current than those in the Rev. B Rockwell 
specification, plus about 3-dB relief from a narrower bandpass filter at the 
detector input, will allow the Ku-band system to provide adequate communica- 
tion tracking. This conclusion assumes that one servo will indeed accommodate 
a 15-dB dynamic range. Axiomatix has made this recommendation to Rockwell and 

NASA. As a consequence, Rockwell will be revising the flux-density specifica- 

2 

tion. The maximum flux density will be specified at -110.5 dBW/M , and the 
two minimum flux-density numbers, including a 2.5-dB pad, will be specified at 
-119.5 dBW/M 2 and -126.0 dBW/M*\ The former applies when the TDRS is accurately 
closed-loop pointing, and accurate SSO pointing will be required at this flux 
density. The latter will be specified for acquisition only and assumes open- 
loop- TDRS pointing. 
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1.0 INTRODUCTION 

In a recent pair of papers [1,2], the authors investigated the 
carrier-tracking performance of a two-channel Costas-type loop for unbal- 
anced quadriphase (UQPSK) signals. Such a loop evolves naturally as a 
closed-loop system motivated by the maximum a posteriori (MAP) phase esti- 
mation approach. Although not explicitly stated, the performance results 
presented m these papers pertain only to the case where both data channels 
have NRZ formats. In many applications [3-5], the data formats on the 
two channels are different; for example, the high data rate channel might 
be Manchester coded, while the low data rate channel might be NRZ, or vice 
versa. Alternately, both channels could be Manchester coded. 

In this report, we generalize the results given in [1] so they 
will apply to any of the mixed-format cases mentioned above. Specifically, 
we examine the mean-square phase jitter performance of two-channel Costas 
loops with active arm filters, synchronous or asynchronous symbol clocks, 
and arbitrary data formats on the two channels. Also considered will be 
the performance of such loops with hard limiters in the m-phase channels. 
Here again, we shall follow the approach taken in [7], with the exception 
that the data formats on the two channels can be arbitrary. 

2.0 ' PERFORMANCE ANALYSIS OF TWO-CHANNEL COSTAS-TYPE LOOP 

WITH ACTIVE ARM FILTERS WITHOUT HARD LIMITERS 

2.1 System Model and Loop Equation of Operation 

To avoid duplication of effort, we shall draw heavily upon the 
notation, definitions and results given in [1], Consider the two-channel 
Costas loop with Integra te-and-dump arm filters illustrated in Figure 1. 

The input x(t) is an unbalanced QPSK (UQPSK) signal-plus-noise, i.e.. 


x(t) = d-j (t) cosC(o 0 t+e^ + x /2P^d 2 (t) sin(u 0 t+ a] + n(t) (1) 




Channel i is NRZ (i=l,g) 

! 1 , 0 < t < T i 

0 , otherwise 


Channel i is Manchester (i=l,2) 

( 1 , 0 < t < T^/2 
P n (t) = <-l ; T/2 < t < T 1 
( 0 , otherwise 


rsa 


Figure 1 Two-Channel Costas-Type Loop with Active Arm Filters 
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where ojq is the radian carrier frequency, and P 2 represent the power 
in the two quadrature components, d-j(t) and d 2 (t) are 1 *' wo independent 
binary waveforms (arbitrary format) with respective data symbol rates 
R-j = 1/Tj and R 2 = 1/T 2 , and e is the carrier phase to be tracked by the 
loop. The additive channel noise n(t) has the usual narrowband process 
representation, i.e., 

n(t) = >^{N c (t) cosC^gt+e} - N g (t) sinC^gt+eDl (2) 


where N (t) and N (t) are approximately statistically independent, sta- 

v J 

tionary, white Gaussian noise processes with single-sided noise spectral 
density N Q W/Hz and single-sided bandwidth B^<Wg/2Tr. 

The input signal is demodulated twice by the quadrature refer- 
ence signals 


r s (t) sin(u 0 t+ e} , r c (t) = yj2K^ cos^t + e} (3) 

to yield the phase detector ^gain ^i<~j outputs 
* c (t) - x(t) r c (t) 

d, (t) + N c (t)] cos <p(t) + ^(l 2 (t)-N s (t)J sin </>(t)j 


V. 


S s (t) s ^x(t) r 5 ( t) 


m 


j^j" d-^t)* N c (t) sin <?(t) - |^~ d 2 (t) - N s (t)J cos <p(t)j (4) 


where $ = Q - 6 is the loop phase error which shall be assumed constant 
(in time) over the duration T 2 of a symbol corresponding to the lowest 
data rate d 2 (t)*. 

The signals e c (t) and e s (t) are each passed through integrate 
and dumpswhich are matched to the pulse shapes p-j(t) and p 2 (t) respectively 
of the two quadrature modulations d-j(t) and d 2 (t) The outputs of like 


As in [1], we arbitrarily assume that d, (t) is the higher 
rate channel. 
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OF POOR QUALITY- 

integrate and dumps are multiplied, and the difference of these 
multiplications is used to form the loop error signal. This error signal 
is then passed through the loop filter F(s) whose output drives the VCO , 
generating the pair of demodulation reference signals given in (3). With- 
out belaboring the details, we may follow the identical approach taken in 
[1] and immediately write the linearized loop (sin 2 = 2cp) equation of 
operation in Heaviside operator notation as* 


- V p) 



(5) 


where 


t-\ - S'(0) F(s) 

% (s) ' s+ S* (O)F(s) * 


( 6 ) 


S ' (0) is the slope of the normalized loop S-curve at the origin {<p - 0} 
given by 


S'(0) = K-jT-j' 


P 1" P 2 d 21 (t) 


+ k 2 t 2 


P 2 -P 1 d 12 (t) 


(7) 


and N(t) = N-j(t) - N 2 (t} is the equivalent zero mean additive noise with 


M ‘> - i N u < t > ! 

1=1 


-N in (t) = K i |N ci 2 (t) - N s1 2 (t)J sin 2 2K, N si (t) N c . (t) cos 2 cp 



sin 2cp+2 ^/P^ P 2 d li^ d 2i ( t ) cos 2 < p[ 


The linearization of the equation of operation as characterized 
by [1 (28) - (31 )] is incorrectly scaled. 
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In (7) and (8), the overbar denotes statistical expectation and 


d Ui(t) 

6 T / kTl d « (s) M s) ds 
(k-nq 

; (k-l)T. 
i=l ,2 


a r kT ,- 

- | N c (s) p^s) ds ; 

(k-l)T. 


(k-DT, 

i=l ,2 

N si<t> 

- J 1 N $ (s) p 7 (s) ds ; 

(k-1 )T 1 


(k-l)T. 

i=l,2 


( 9 ) 


Note that (9) differs from the corresponding expressions in [1] in that 
the channel pulse shapes p-j (t) and p 2 (t) are included as weighting fac- 
tors in the integrands. 

Assuming that the bandwidth of H 2< p(s) (loop bandwidth) is much 
narrower than the bandwidth of N(t), the variance of Zcp can be approximated 
by 



V \ 

Cs'(o )] 2 


( 10 ) 


where is the single-sided loop bandwidth and Ng' is the single-sided 
spectral density of the equivalent noise N(t). 

The next two sections are devoted to an evaluation of S * ( 0) and 
Ng 1 for the various data combination formats on the two channels 

.2.2 Evaluation of the Loop Signal Component S'(0) 

From (7), we observe that evaluation of S’(0) reduces to the 
evaluation of the two mean-square crosstalk quantities d 2 ^ 2 (t) and d-j 2 2 (t) 
for the four combinations of NRZ- and Manchester-coded data on the two 
channels. Such evaluations have already been performed in [6] under the 
guise of a different notation. In particular, [6, Table 1] evaluates what 
is herein referred to as d 2 -j 2 (t) for the case of synchronous symbol clocks, 
and [6, Table 2] gives the corresponding results for asynchronous symbol 
clocks. In the former case, the ratio of data rates R-j/R 2 is constrained 
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to be an integer. Translating these results into the notation of this 
report, we obtain the following tables. 


p 

Table 1.* Evaluation of d 21 (t); Synchronous Symbols Clocks, 
R-j/Rg = n (n an integer) 


d 1 (t) 



Manchester 

NRZ 

Manchester 

0 ; n even 

1 ; n odd 

1 ; n even 
; n odd 

NRZ 

0 

1 


2 

Table 2. Evaluation of d 2 -j (t); Asynchronous Symbol Clocks; 
Ri /R 2 = n (n need not be an integer) 


d-j(t) 



Manchester 

NRZ 

Manchester 

l In 2 - £■ + 12 - 6n ; 

o n 

R 2 < R] < 2R 2 

in ’ R 1 > 2R 2 

„ 1 n 2 , , 1 
n " 6 n “ 1 3n 5 

r 2 < Ri ^ r 2 
1 - I ; R, > 2R 2 

NRZ 

1 

6n 

>-ir 


Also, from the definitions of d 12 (t) and d 21 (t) in (9), it is 
straightforward to show that 


*Note that [6, Table 1] contains an error for the Manchester- 
Manchester case. 
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d 12 2 ^ " n d 21 Z ^ 

for all entries in both Tables 1 and 2. Thus, further introducing the 
notation [1] 

A A Po A P Tl 

p i = p i + p 2 * x - pf • 'i - « 2 > 

the loop signal component of (7) simplifies to* 


S’( 0 ) 



P(l+x) -1 ^ -X d 


21 



n 



(13) 


2 

where 6 ^ is obtained from the appropriate entry of Table 1 or 2. 

2.3 Evaluation of the Equivalent Noise Spectral Density Nq' 

The single-sided noise spectral density Nq ' of the equivalent 
noise N(t) whose components are defined in (8) is obtained from 


N 0 ' - 2 p R n (t) dt (14) 

—oo 

where 

R n (t) = < N(t) N ( t+r ) > = < R N (t,t+x) > (15) 

and < > denotes the time average. Using (8) and taking the statistical 
average required in (15) gives [1] 


* 


Herein we eliminate the notational dependence of d^-j 


^(t) on t. 
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hih’h) ■ K « 


4 ^ 2 _ 2 
a £ +a £ T £ 


(v p a7) 


+ T £ 4p i P 2 d ££ Z 


Ijj, .tg) 


£=1 L 
h 


“ J 



" 4K 1 K 2 Z| R N c 1 N c2 (v**) + T 1 T 2 




X 

_ PlRd ll d l 2 ( trt i) +l>2Rd 21 d 22 

( t £ ,t £) 

Vc*' 

( t £* t £) 

+ T iV p i p 2 

R d d (^f M ^'£ 
a 21 22 V * - 

)J (16 > 


where we have introduced the simplifying notation 



1 ; l = 2 

2 , £ - 1 


Also in (18), 


2 A N 0 T £ . 
ff £ 2 


£ = 1,2 


and 

■A.tjP 


1 ; if t^ and tg are both in the interval 

K> < k+i > t «) 

0 ; otherwise 

£ = 1,2 


(17) 


08) 


09) 


Clearly then, we must evaluate, for each of the mixed data format combin- 
ations, the various nonstationary correlation functions required in (16) 
before performing the time-averaging operation which finally results in 

^|( T ) • 

We begin by presenting some general results which relate these 
correlation functions to one another, namely, it can be shown that analo- 
gous to (11), 
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R d d ^1 *^2^ = rT d »to3 
a H a 12 1 L n a 21 22 1 1 


where again n is the ratio of data rates. Furthermore, 


R M N = ^ d C^l iXn) 

N cl N c2 1 6 1 d 21 a 22 1 1 


The above two relations are valid independent of the data formats on the 
two channels and whether or not these data streams are time synchronized. 
Using (20) and (21), we can determine the products of correlation functions 
needed to evaluate (16). Thus, we obtain 


r n n ^i ^ 2 ) 
”cl w c2 


a l R d d ^1 *^ 2 ) 
I d 2l a 22 1 L 


R d d ^1 ’^2-^ N 
a ll a 12 1 L "cl"c2 1 c 


n ^2^22^1 


R d d ’^2-^ N ^1 5 ^ 2 ) 
a 21 a 22 1 * N cl N c2 1 L 


^d^V 1 ^ R d 21 d 22 ^l 5t 2^ 


a l d C fc i ^ 2 ) 
1 a 2l a 22 ! L 

w R d 21 d 2 /vy 


from which we observe that it is sufficient to evaluate only Rd2]d 2 2^ t l 
for the various data format combinations on the two channels. These spe- 
cific relationships will now be given in summary form. 

Channel 1 is NRZ, Channel 2 is Manchester; Synchronized Symbol Clocks, 
n an integer 


! 1 ; 0 < t-j ^ Ti , 0 ^ t 2 ^ nT-j 

1 ; T, < t, < n Tl , nT-| < t, « 2nT, (23) 

0 ; otherwise 



original pag£ is 

OF POOR QUALITY 


10 


Channel 1 is Manchester, Channel 2 is NRZ; Synchronized Symbol Clocks; 
n an integer 


R d d Ct r O = 0 
a 21 a 22 1 d 


all t.j ,t 2 


( 24 ) 


Channel 1 is NRZ, Channel 2 is Manchester; Synchronized Symbol Clocks; 
n an integer 

n even 


R 


d 21 d 22 Ctl 


t,) 



0 < t 1 < T 7 , 0 < t 2 < nT 1 

T 1 < t ] < (J+ l) T 1 , nT 1 < t 2 < 2nT 1 

(f + ] ) T 1 < t l < nT l ’ nT l <t 2 < 2nT l 
otherwise - (25) 


n odd 


R 


d 2 i d 22 ^ t l st 2^ 



;0<t,<T 1 , 0 < t 2 < nT, 

; T 1 < t-j < (^rjT-| > nT ] < t 2 < 2nT 1 

’ ( i r I ) T i < *i < (ir) T i ■ nT i <t 2 < 2nT i 

; (^)t 1 < t 1 ^ nT-j , nT-j < t 2 < Zr\T } 

, otherwise (26) 


Channel 1 is Manchester, Channel 2 is Manchester, Synchronized Symbol 
Clocks; n an integer 


n even 



n 


n odd 


R 2 

d 21 d 22 
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0 ; 0 < t-j < T| , 0 < t 2 < nT-j 

0 ; T-j < t, < , nT-, < t 2 < 2nT-| 

1 ; (t 1 ) 1 ! < t l < ( I ^) T 1 » nT l < t 2 < 2nT 1 
0 ; (¥) T 1 < t] < nT ] » nT 1 < t g < 2nT 1 . 

0 ; otherwise (28) 


Letting t-j = t, t^t+r and averaging over t yields the following 
results for the stationary correlation function 


A(X) = - ]r~ r T} R 2 • (t.t+f) dt 

ni l^ d 21 22 

Channel 1 is NRZ, Channel 2 is NRZ; Synchronized Symbol Clocks, 
n an integer 


(29) 


a(t) 


nT-j CT 1+ x) , 
^[(20-1)7,-,]; 
0 ; otherwise 


< t < (n-l)T-, 

(n-1 )T-| < t < (2n-l)T 1 

(30) 


Channel 1 is Manchester, Channel 2 is NRZ; Synchronized Symbol Clocks, 
n an integer 


A(t) = 0 ; all t 


(31) 


Channel 1 is NRZ, Channel 2 is Manchester; Synchronized Symbol Clocks, 
n an integer 

n even 


A(t) 



CvO ; -T, 

(2n-l)T r 3; 

otherwise 


< t < (n-1 )T ] 
(n-DT-j < r < 


(2n-l )T-j 
(32) 



n odd 
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A(t) 


nT-j 0 

i + T ) * 

- T 1 < T < (^h 

n-1 

2n 

* 

9 


T 

nT i 

* 

9 

(t 1 ) 1 ! 4 * 4 (n - ,)T l' 

i 

nT l 

(2n-2)T 1 -x 

; (n-DT, < t < (^)t i 

n-1 

2n 

9 

( 3 T 3 )T, < * < ( 3 2 ’) T 

1 

nT l 

(2n - 1 )1 } - t ; < t < (2n-l )1 


0 ; otherwise 


(33) 


Channel 1 is Manchester, Channel 2 is Manchester; Synchronized Symbol 
Clocks; n an integer 

n even 

a(t) =0 ; all t (34) 

n odd 



yields 


Making the same substitutions in (16) and averaging over t 
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R N (t,t+T)dt 


a i A+a i J i( P i +P l d ££/ +T £ 4p l P 2 d ££ 


£= 1 



4K-,K 2 


4 2 

a l + cf l T 1 T 2 




P 1 P 2 



(A(t)+a(-t)}(36) 


where we have also made use of the relations in (22). Finally, perform- 
ing the doubly infinite integration on t gives the desired results for 
Ng' of (14), namely 


i— 

v - 8 I k A 


£=1 


16 K-j K 2 


a i * a i J l ( P £ + P £ d ££ ) + T £ 4 P 1 P 2 d ££ 2 


CT 1 4 + °1 2 T 1 T 2 ( 1T + P 2 ) + T 1 2j 2 2 P 1 P 2 W 


ICa(t)3 (37) 


where 


iCMtJD 4 f 

—CO 


A(r)dr 


(38) 


All that remains to put (37) into its final form is to evaluate I(a(t)} 
for the various mixed data format combinations. 

Substituting (30)- (35) into (38) and carrying out the integra- 
tions enables us to put the results m a surprisingly simple and compact 
form, namely. 


I&M3 ■ nT, 


(39) 


where d 2 -| 


2 


is obtained from Table 1. 
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Finally, then, substituting (39) into (37) and making use of 
(11), we obtain, after some simplification. 



Although the result in (40) has been derived for the case of 
synchronized symbol clocks, with some effort, it can also be shown that 
it al so a pplies for the unsynchromzed symbol clock case with the caveat 
that dgi 2 is now found from Table 2. 


2.4 Squaring Loss Evaluation 

The "squaring loss" of the two-channel Costas loop of Figure 1 
is defined by 


4 * 


4 4N 0\ 


where is found from (10). As is traditional in suppressed-carrier 
types of tracking loops, this loss represents the degradation relative to 
the mean-square phase jitter performance of a linear loop caused by SxS, 
SxN and NxN distortions. Using (10) m (41) allows to be rewritten 
in the form 
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Substituting (13) and (40) into (42) and simplifying using the definitions in 
(12) gives the desired result, namely, 



or using the optimum* ratio of the gam constants 



p 2 t, _ 
P 1 T 2 ' 


X 

n 5 


(44) 


we have 



(45) 


Comparing the above result with [1, (38)], we observe that the only modifica- 
tion necessary to make the previous result applicable to the more general data 

2 

format situation is to calculate the single parameter d,^ from either Table 1 
or Table 2 and use it where appropriate in (45). Also observe that by letting 
l<2=0 in (43), we get the result for a single-channel Costas loop, namely. 



which agrees with [6, (. 32 )] if the appropriate changes in notation are made. 
* 

The word "optimum" is used here in the same sense as that used in 
[1], namely, the ratio of gains motivated by the MAP estimation approach. 
Note that this gain ratio does not necessarily minimize squaring loss or 
mean-square phase jitter, as demonstrated in the next section. 
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2.5 Optimization of Channel Gain Ratio to Minimize 

Wean-Square Phase Jitter 

In the previous sections, we developed a theory for calculating 
the mean-square phase-tracking jitter based on a linear loop model as a 
function of such system parameters as signal -to-noise ratio in each chan- 
nel, the ratio of powers in the two channels, and the ratio of data rates 
in the two channels. Although the general results were derived as a func- 
tion of the ratio of the gains in the two channels, specific results were 
presented only for the "optimum" choice for this ratio based upon MAP 
estimation theory. 

Unfortunately, as we shall see shortly, selection of this gain 
ratio as motivated by MAP estimation theory does not necessarily optimize 
the tracking performance. In fact, selection of this gain ratio to directly 
optimize the tracking performance in the sense of minimum mean-squared 
phase-tracking jitter can yield considerable improvement in this perfor- 
mance relative to the MAP estimation choice of gain ratio. 

To illustrate this point, we shall begin by rewriting (41) in 

the form 


where 



4N o b l 



V A i( 

k 2 \- 

K 1 L 


O + A) 

B 0 +B l 

/ K 2\ n (hf 



(47) 


(48) 
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Clearly, from (47), minimizing the phase-tracking jitter a Z(p is directly 
equivalent to maximizing S L (minimizing the squaring loss). From the form of 
(47), it is a simple matter to differentiate this equation with respect to 
K 2 /tC| and thus find that value of channel gam ratio which maximizes 
After some routine algebra, the solution to the equation 




(51) 


Note that the first factor in (51) is the "optimum" value of as motivated by 

the MAP estimation theory, namely, that given by (44), which is independent of 
signal -to-noise ratio. Only when the channels have equal energy, i.e., X-l/n, are 
the results of (44) and (51) equal, i.e., K^/K-j = A . 

Substituting (51) into (47) results in an expression for the minimum squaring 

loss, namely. 



where A Q , A-j , B Q , B-, and B 2 are found from (48). 
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As an example, when both channels have NRZ data and synchronous 



symbol clocks, from Table 1, we then have that d ^ =1 and (52) simplifies 
to 



(53) 


2 

On the other hand, letting d^-j = 1 in (45) gives the result obtained in 
[1], namely. 



Figure 2 illustrates as computed from (53) and (54) versus 
the power ratio A with n=10 and ^ = -3 dB and 10 dB. It is clear from 
this figure that considerable reduction of squaring loss is possible by 
choosing the channel gam ratio so as to minimize this quantity. 

Another illustration of the squaring loss improvement obtained 
by optimizing the channel gain ratio is given in Figure 3 Here we con- 
sider a mixed data format where the high data rate channel (channel 1) is 
Manchester coded, while the low data rate channel (channel 2) is again 
NRZ. Furthermore, we now assume asynchronous symbol clocks, which implies 
-that the value of d^ 2 to be used in (48) and, hence (52), is obtained 
from Table 2, namely, l/6n= 1/60 (for n=10). Comparing Figure 3 

w ith Figure 2, we observe that the small value of cross-modulation loss 
(^12^) ^ or t * ie rn ' ,xec * f° rrna t case (Figure 3), as compared with a unity 
value for the case where both channels have NRZ data (Figure 2), results 
in a considerable reduction in squaring loss for both the optimum and MAP 
estimate selections of the channel gam ratio. 
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Figure 2. Squaring Loss versus Power Ratio for Unbalanced QPSK, 
Two-Channel Costas Loop; n = To/T 1 - 10 
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3.0 PERFORMANCE ANALYSIS OF TWO-CHANNEL COSTAS-TYPE LOOP 
WITH ACTIVE ARM FILTERS AND HARD LIMITERS 

3.1 System Model and Loop Equation of Operation 

Here again, we draw heavily upon the notations, definitions and 
results given in £1]. The loop is still represented by the block diagram 
of Figure 1 except that now G(x) = sgn_x. Thus, the dynamic phase error 
z(t) now becomes 


z(t) 


yioj* z sl (t) sgn jyiGj" z cl (t)J - yicf z c 2 ) sgn j/icT 
yKfVlOO sgn [* cl (t)] - z c2 (t) sgn jz $2 (t)J (55) 


where 


z si (t) 


C kT 

J e s (t) P 7 (t)dt ; 

(k-l)T. 


kT n < t < (k+l)T- 
1 = 1,2 


z ci (t) 


J e c (t) p-(t)dt ; 

(k-l)T. 


kT < t < (k+1 )T 

i=l ,2 


(56) 


are the integrate-and-dump outputs. Splitting z(t) into its random and 
nonrandom parts, i.e., 

z(t) = E{z(t)foj + H(t,<p) = S (cp) + N(t,<p) (57) 

where S{<p) is the loop S-curve and N(t,<p) is the equivalent zero-mean 
additive noise process, then the loop equation can be written in Heavi- 
side operator notation as 


P <P = - K„ F(p) z(t) = - K 


F(p)[s(<p) 


+ N(t ,<p) 


(58) 


where K y is the VCO gain in rad/s/v. Finally, linearizing the loop by 
letting 
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S(<p) = S'(0)<p ; N(t,<p) = N(t,0) 

we get an equivalent form analogous to (5), namely. 


v - M p) f¥7^ 


where 


S'(0)K v F(s) 

*V s) = s + _ s r (oJTc v "FO 


Assuming, as before, that the bandwidth of H<p(s) is much narrower than the 
bandwidth of N(t,0), the variance of <p is approximated by 


V B L 
ls'(o)Y 


where Nq'is the single-sided spectral density of N(t,0) and B^, the 
single-sided loop bandwidth, is simply 


\ = hf~ KM 2 d " 


3.2 Evaluation of the Loop Signal Component S 1 (0) 

. Comparing (55) with (57), we immediately observe that the loop 
S-curve is defined by 

S M = ^S } (cp) - yi(JS 2 (<p) (64] 


where 


Sjfa) = E|z sl (t) sgn jz cl (t)J 

S 2 ((p) - e | z c 2^^ s 9 n 


( 65 ) 
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Substituting (4) into (56) and carrying out the integration gives 
([assuming KgK m = l} 


hiM - 

T i ^7 d n (t ) + N ci(t) 

sin cp - 

Ji'/ i T d 2i' t > + N s 1 < t >] 

z ci (t) = 

T i7 p 7 d ii (t) + N ci {t) 

cos cp + 

_T i VPi‘<l 2l (t) + N sl (t)] 


i-1,2 (66) 

Performing the statistical expectations first with respect to the noise 
components in - { t ) and z cl (t), keeping d-j n (t) and d 2l -(t) fixed, gives 
the conditional S-curve components 

S i(^| d n ,d 2i) = e n c1 ,n s1 { z sl sgn ( z cl) K d ll ,d 2l} 

- rfi d n sin ^ “ T i 7^2 d 2 i C0S(p ] 

x erfc j - T-j yF j" d-^ cos (p - T-| d 2 -| sin (p / /2 o-jj 

- \ T i d ll sin ^ - T ] v/P^ d 21 cos ^ 
x erfc j T-j cos <p + T-j yp[J" d^ sin <pj/ /Z 

(67) 

where erfc x is the complementary error function defined by 

exp(-y 2 )dy (68) 


erfc x = 


/7 


/ 


Since 

i 


erfc(-x) - erfc(x) = 1 + erf x - (1 - erf x) = 2 erf x 


(69) 


then 
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s l(<p| d ir d 21 ) - T 1 jjfP^d 11 sin (p- JP^d 21 cos ^ erf 

(70) 

where ^ is given by (18). Similarly, 

S 2 ^| d 12 , d 22 ) = E N c2 ,N s2 | z c2 sgn ( z s2)i^ sd 12’ d 22| 

= I 2 ^[^2 cosg?+ ^/^ 2 d 22 S1 " n ^] erf 

(71) 

Since evaluation of a ^ from (62) requires only the slope of the loop 
S-curve at the origin, we next proceed to differentiate (70) and (71) with 
respect to cp and evaluate these conditional results at <p - 0. Thus, 


KP- 0 




r i|/7 d n sin <p| 


Similarly, 


- MVll erf 


Aijr d n\ 

l / 





' T i/7 d n erf (i^ d n) 






(72) 



(73) 
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Since, from their definition in (9), and d 22 take on only 
±1 values with equal probability, then (72) and (73 ) simplify to 



Finally, then, averaging over d^ and d-j 2 and recalling (11) gives the 
desired result, namely, 



For the ratio of channel gains as specified by (44), (75) further simplifies 
to 
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Evaluation of the Equivalent Noise Spectral Density N Q ' 

To evaluate Nq 1 , we must first find the correlation function of N(t,0), 


R m (t) = <N(t,0) N(t+r,0) > (77) 

then integrate it between -■» and «, as in (14). Since the value of the loop 
S-curve S{<p) at (p = 0 is zero, from (57), we then equivalently have 

R n (t) = <iTtlF0 T z(t+x|(P=0) > (78) 

Letting <p - 0 in (66) and substituting in (55) gives 

z(t|¥>=0) • +N $1 (t)] sgn Jt, ^ d,, (t) +N cl (t)] 

+ \Fj> [t 2> /FY d 12 (t) + N c2 (t)J sgn [j 2 -v/PJ d 22 (t) + N s2 (t)l (79) 

Substituting (79) into (78) and carrying out the statistical averages over the 
various noise components and data process results in, after much simplification, 


(t-, |^0)z(t 2 ^ 0 j = K-j 


T 2 d . 2 X 2 
T 1 P 2 d 21 +CT 1 


~A l 


T r^2 R d 21 d 22 


i ( l i ’‘ 2 ) + 


T 2 d , 2 X 2 
T 2 P 1 d 12 + a 2 


A, 


l^l ^ 2 ) 


2a / T 2 P 
+ — — expf ^-|B 


/ ^ \ 20, 


(*1 ^ 2 ) 


r 2 V^T R dndi 2 (*1 • t z) ert V^J- 


T-i fP7\ 2a ? / T^P, 

™ U — L exp /_ 


/2irn \ 2a 




T l \ fl ^ 2 ^ 22 ^ ^ 2 ) 


t zM, 2 «i 


1 

T 2 2p 2! 


2 rn ’ 2 


erfl-^-^ 1+ — — exp|- -^MB 
■J2 a 0 1 /2irn 




T i /TV Rcr 2 

+ exp 


2 a, 


t, 2 p 


(vs) 


v^rrn \ 2 a 


- ^ "(v 

2a. < V ^ 


( 80 ) 
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where 


A i (h ft z) 
B (*1 }t 2 ) 


1 ; t-j and t 2 are both in the interval 
(kT. ,(k+l)T.j for some k 

0 ; otherwise 

1 = 1,2 


d 21 d 22 



(81) 


Using (11) and (20), we are able to simplify (80) whereupon, averaging over t-j 
with t 2 = t-j + t, we get 


R n (t) = •<! 


T 

I 1 P 2 





x 


T 2VT erf 





(82) 


where a(t) is defined in (29) and evaluated in (30)-(35). Finally, integrating 
t as required in (14) gives the desired result, namely. 



x 


(83) 




or for K^/K-j as in (44), 
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3.4 Squaring Loss Evaluation 

Analogous to (41), the squaring loss of the two-channel Costas 
loop with hard limiters in its in-phase arms is given by 



Vl 

•V 2p 


where a, 


W 


is found from (62). 


Making this substitution gives 


( 85 ) 


= [S'(O)]^ (86) 

p(n 0 '/n 0 ) 

Finally, substituting (75) and (83) into (86) gives the desired result, 
namely. 




(87) 





where 
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For Kg/K-j = X/n as in (44), (88) simplifies in the same manner that (75) 
becomes (76) and (83) becomes (84). 

3.5 Optimization of Channel Gain Ratio to Minimize 

Mean-Square Phase Jitter 

As for the double-channel Costas loop without hard limiters, here 
we can choose the channel gain ratio to minimize mean-square phase jitter. 
.From (87) and (88), we observe that the squaring loss can be written in 
the form 



( 89 ) 



where 
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From (40) and (50), it is obvious that the gain ratio that minimizes the 
phase-tracking jitter is 



with the corresponding minimum squaring loss 


(91) 


J L 


For the case of NRZ data on both channels with synchronized sym- 
bol clocks. Figure 4 illustrates as computed from (92) versus the power 
ratio X with n = !0 and ^ = -3 dB and 10 dB. For comparison purposes, ^ 
as computed from (88) with K^/K-j^x/n is superimposed on these same curves. 
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Clearly, optimization of the channel gain ratio in accordance with (91) 
offers a significant reduction in squaring loss relative to the values 
obtained by choosing this ratio in accordance with the MAP estimation 
relationship of (44). 

To further illustrate the effect of power unbalance on squaring 
loss performance, we consider the case of equal bit energy in the two chan- 
nels, namely, 


P 1 T 1 • P 2 T 2 ' E b 


(93) 


or, equivalently, from (12) and the definition of n, 


nA = 1 


(94) 


Substituting first (94) into (87), then letting = A/n = A 2 gives 



erf /? - — v'ca ( d 2i 2 ) ex P( - s) 

YTt _ 

2 

~M 2 ~ 

S0+A)( d 21 2 ) + 1 - 2/X ( d 21 2 ) 

i/c erf /z. + 
_ 

— exp(-c) 
/rr 

,2 

2(1 + A) 



(95) 


where 





p i T i 




(96) 


Equation (95) represents the generalization of [1, (70)] for the case of 
arbitrary data formats on the two channels and either synchronous or asyn- 
chronous symbol clocks. Figure 5 is a plot of as determined from (95) 
versus A with E^/Nq as a parameter and NRZ data on both channels. 

The corresponding curves for the optimum channel gain ratio are 
obtained by substituting (94) into (93) and also are illustrated in Fig- 
ure 5. Unlike the two-channel Costas loop without hard limiters, the two 
results here are not, in general, equal when the channel energies are made 
equal. In fact, for large unbalance in power and large values of E^/Nq, 



8 


0 


Figure 5 


Squaring Loss versus Power Ratio for Unbalanced QPSK; Two-Channel 
Costas Loop with Hard Limiters, n = 1/A (Both Channels have Equal 
Bit Energy), Both Channels are NRZ Data with Synchronous Symbol 
Clocks 




34 


much benefit is derived by optimizing the channel gain ratio. Only when 
both the power and data rate ratios are equal to one (i.e., balanced QPSK) 
do the two results agree exactly as one would expect. 

Figure 6 is another example of the effect of channel power 
unbalance on squaring loss performance when the bit energies in both chan- 
nels are equal. Here we assume that channel 1 has a Manchester data format 
while channel 2 is still NRZ data. The symbol clocks, however, are now 
assumed to be asynchronous. Thus, the appropriate cross-modulation loss 
factor is determined from Table 2 as d^Z = (6n) = x/6. For large power 

unbalance (small X), the cross-modulation loss is small and, thus, we observe 
an improvement in squaring loss performance relative to the previous case 
where d^ = 1 for all X. 



Figure 6. Squaring Loss versus Power Ratio for Unbalanced QPSK; 

Two-Channel Costas Loop with Hard Limiters, n = 1/x 
(Both Channels have Equal Bit Energy), Channel 1 is 
Manchester, Channel 2 is NRZ, Asynchronous Symbol Clocks 
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INTRODUCTION AND DISCUSSION OF PROBLEM 


The purpose of this report is to determine the bit error rate 
(BER) degradation for the Ku-band return link mode 2. Figure 1 illustrates 
the overall Ku-band return link. The typical return link data sources are 
indicated in Table 1. In mode 2, there are three inputs to the FM modula- 
tor. First, there is the bent-pipe wideband signal which is amplitude 
(peak) regulated so as to provide 22 MHz peak-to-peak carrier deviation, 
which is channel 3. In addition, there is channel 1, which has operations 
data from the network signal processor, and channel 2, which receives data 
from the payload recorder, operations recorder, payload low data rate link 
or low data rate payload interrogator. Channels 1 and 2 together provide 
for a peak-to-peak FM frequency deviation of 12 MHz 


Table 1. Ku-Band Return Link Data Sources 


Source 

Type 

Rate or Bandwidth 

CHANNEL 1 (MODE 1/MODE 2) 



Operations data— Network 
Signal Processor (1,2) 

Digital 

192 kbps (biphase) 

CHANNEL 2 (MODE 1/MODE 2) 



Payload Recorder 


25.5 - 1024 kbps (biphase) 

Operations Recorder 


25.5 - 1024 kbps (biphase) 

Payload Low Data Rate 


16 - 2000 kbps (NRZ) 

16 - 1024 kbps (biphase) 

Payload Interrogator (1,2) 
Low Data Rate 

Digital/ 

Analog 

16 - 2000 kbps (NRZ) 

16 - 1024 kbps (biphase) 
1 .024 MHz (PSK or FM) 

CHANNEL 3 (MODE 1) 



Payload MUX 

Digital 

2 - 50 Mbps (NRZ) 

CHANNEL 3 (MODE 2) 



Payload Interrogator (1,2) 
High Data Rate 

Digital/ 

Analog 

16 - 4000 kbps 
0 - 4 5 MHz 

Payload 

Analog 

0 - 4.5 MHz 

Payload High Data Rate 

Digital 

16 - 4000 kbps 

Television 

Analog 

0 - 4.5 MHz 






Bent-Pipe 
Wideband <; 
Signal 

NSP Data , 


Ch3 

Tpf 

lBW-5 MHz 


Chi 


LPF 

BW=5 MHz 


Bent-Pipe 
Narrowband 
Signal 


Ch2 
LPF 
BW=5 MHz! 


Amplitude 



LPF BW 


iNLruanu i i ariom i l 

Regulator 


r 



6.7 MHz' 












Limiter 







i 





— 

20% 

QPSK 

Subcarrier 
Modul ator 


BPF 
BW = 

8 MHz 

y 


Carrier 

Frequency 

Modulator 







j 


1 



H 







_2sl 



FM Receiver BPF 


Frequency 

Demodulator 
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Tracker 


Bent-Pipe 

Narrowband 

Signal 

Processor 




I Tmi fpr 

Filter 
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Figure 1. Ku-Band Bent-Pipe FM Return Link 
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Both channels 1 and 2 are hard limited by TTL receivers 
(essentially limiters) and QPSK modulated onto an 8.5 MHz subcarrier. 

This quadriphase subcarrier is added to channel 3 at baseband and frequency 
modulated onto the carrier. Figure 2 gives a more detailed block diagram 
of the modulation process. In Figure 2, both modes 1 and 2 are illustrated 
however, in this report, we are concerned with mode 2 only, the FM mode. 

The specific case of concern in this report is the case wherein 
the signal of channel 1 is a 192 kbps biphase-L (Manchester) baseband data 
stream without thermal noise corruption, channel 2 has a 2 Mbps NRZ base- 
band data stream plus thermal noise, and channel 3 is composed of two 
subcarriers, one at 1.024 MHz and one at 1.7 MHz, plus thermal noise. 

The power spectral density (PSD) for the case in which channel 1 
is 192 kbps Manchester, channel 2 is 2 Mbps NRZ and the 1.7 and 1.024 MHz 
tones are unmodulated and are sketched in Figure 3 along with the respec- 
tive lowpass and bandpass filters (into the summer) of Figure 2. 

2.0 CHANNELS 1 AND 2 QUADRIPHASE SIGNAL 

In this section, we discuss the Hughes Aircraft Company (HAC) 
implementation of the quadriphase modulator. A digital subcarrier oscil- 
lator (SCO) running at 8.5 MHz supplies a clock and its complement. 

Both clocks are passed through a fixed-time delay (t) corresponding to a 
fixed-phase delay 2e Q in the square-wave fundamental. When not filtered, 
the result is a constant-envelope four-phase signal, as shown on the I-Q 
diagram of Figure 4 that has a ratio of I Q of 4:1 The actual signal 
wave is filtered so that only the first harmonics are applied to the FM 
modulator. 

In order to understand the power relationships of the I and Q 
- channels, we establish the form of the filtered quadriphase signal follow- 
ing the work of Cager [3] in the case of no noise. The digital subcarrier 
phase modulator is described by 



where 


<f>(t) = ed 1 (t)d 2 (t) - J Cd 2 (t)-l3 


( 2 ) 



4-100 Mbps 



-e» 


Figure 2 Ku-Band Return Link Transmitter-Modulator 
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Figure 3. Sketch of the FM Input Spectra Along With the Associated Separating Filters 
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8.5 MHz 
(Square wave) 


Q (Channel 1 ) 



Time Delay ►Phase Delay 2e Q 

Figure 4 QPSK Subcarrier Oscillator Modulator 
for Mode 2 Ku-Band Return Link 
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with d-j (t) the channel 1 baseband modulation and d 2 (t) the channel 2 
baseband modulation so that the four phases of <j>(t) are described by 
Table 2 below. 


Table 2. Square-Wave Oscillator Phase Relationship 


f(t) 

d ] (t) 

d 2 (t) 

6 

1 

1 

-0+1T 

1 

-1 

0+ir 

-1 

-1 

-G 

-1 

1 


The value of e Q is the value of e such that the I and Q rela- 
tionship yields a power ratio of 4:1 in the first harmonic. Expanding (1) 
in a sme-Founer series yields 



where E denotes summation on the odd terms. Expanding (3) yields 


s(t) = 


• n=l 


co . 

^ cosfncKt)] sin nu) sc t + £ ^ sin[n^(t}] cos nu^t 


nir 


n=l 


nir 


( 4 ) 


Using (2} for <j>(t) yields 


cos[n(j>(t)] 
sin[ri(f) (t) ] 


cos 

sin 


nd.jd 2 

nd l d 2 



since sin 


in ndd-p-O^/^ 


( 5 ) 

( 6 ) 


equals zero. Simplifying (5) and (6) yields 
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cos[n$(t)] = d 2 cos ne n odd (7) 

sin[n$(t)] = d^ sin ne n odd (8) 

Therefore, the quadn phase signal is of the form 


d 2 (t) | cos(ne) sinCn» sc t} + d^t) j, sin(ne) cos(n» sc t) 

n_T n ’' (9) 


U(t) » E cos ne s1n(nu t} 
n-1 niT 5 


w - 

V(t) = 2 sin ne cos(mo t) 

n=l nir sc 


s(t) = d 2 (t) U(t) + d-j(t) V(t) 


To identify the waveforms expressed by (9), note that 


A 

U(t) = 2 zr- cos(ne) sin(nw t} = 

n*i m 


i(a_ r t- 0 + ^in(V r t + e] 

— 2 — 03 ) 


+ 5™C“ c r t + 


where 


inQo t+e) is a hard-limited sine wave of the argument (w t+e}. 


Since 


Wt-03 + Wt + e)= f ■£- cos(-ne) sinCna) 0 + sin(-ne) cosCnoi t) 


sc “-J ntr 


+ E ^ cos(ne) sinCnw sc t3+ sin(ne) cosCna> sc t] 
n=l 


( 14 ) 
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It can also be shown in the same manner that 


V(t) = 2 ^ sin ne cosCna) Qr t) = 


SinQo sc t+e3 - [] inCu sc 


t»e} 


n=l 


sc 


05) 


Now the right-hand sides of (13) and (15) can easily be plotted as shown 
in Figure 5. Thus, we see that the digital quadriphase signal can be 
broken down into two orthogonal basis functions that span the space of 
periodic functions Note that the sum of U{t) and V(t) yields a square- 
wave signal that is early (from t = 0) by 2e/co seconds, whereas U(t) - 

b U 

V(t) generates a square wave with a delay of 9/w seconds, etc. 

b L 

Note that the total power in the waveform of (9) is given by 


P s = i I Ccos 2 (ne) + sm 2 (ne)3 = 4- f X- = 1 
s * n=l nV it n=l n 


(16) 


Furthermore, the power in channel 1 is given by 


P V = P,= X i-^n 2 ».-£ 


0 > 0 


07) 


and the power in channel 2 is given by 


P = P ? = f \~T? cos2 ne = 1 - — 
U 2 A 2 .2 2 tt 

n=l n ir 


\> e > 0 


08 ) 


Notice that, when 6 = 0.464 radians (26 58°), the power split is 
P-j = 0.295 (19) 

P 2 = 0.705 (20) 

and P 2 / p i 15 2.38:1 However, since the BPF at the output of the MUX 
(quadriphase modulator) passes only the fundamental, the expression for 
the respective channel powers become 




Figure 5 Two Basis Functions for QPSK Channels 1 and 2 Signal 
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= A- sin 2 e 

(21) 

Tf 


o p 

= ^ COS 0 

(22) 


ir 


and, for 6 = 0.464, we have that P 2 /P-j i$ 4:1, as was desired. 

We note from (9) using only the first harmonics that the signal 
can be expressed as 

s f (t) = — d^(t) sine cos(w sc t3 + ^ d 2 (t) cose sinQu^t} (23) 

where d-j (t) is the output of the channel 1 limiter and d 2 (t) is the output 
of channel 2. Thus, we see that we can model the bandpass process in terms 
of its baseband processes. 

3.0 IMPERFECT CARRIER REFERENCE LOSSES 

Since we have a quadri phase signal with both a low rate NRZ wave- 
form and a high rate Manchester waveform, we must consider two cases of 
cross-channel interference: (1) a high rate NRZ interferes with a low rate 

Manchester data stream, and (2) a low rate Manchester data stream inter- 
feres with a high data rate data stream. The two data streams are illus- 
trated in Figure 6. Besides the cross-channel interference, the imperfect 
carrier reference loss also exists and is, in general, not negligible. We 
essentially follow the approach of Osborne [8] in the following. 

First consider the case in which the low rate Manchester signal 
interferes with the high rate NRZ signal. We model the received RF quadri- 
phase signal by 

NRZ Manchester 

r(t) = d 2 (t) sinCo sc t+e 0 3 + n | 2P^' d^t) cosC^^+Sq} + n(t) (24) 

where is the channel power in the ith channel, d^(t) is the ith data 
stream (i=l,2), and tu sc and 8 Q are the angular frequency and phase of the 
8.5 MHz subcarrier frequency, and finally, n(t) is white Gaussian noise. 

The pulse shapes P 1 (t) are either NRZ or Manchester symbols. The demodu- 
lation circuit is shown in Figure 7. 



Manchester 


Figure 6 


NRZ -R 


Data Sequence Model for NRZ and Manchester Cross-Channel Interference 


n.| ones 


■*+ 


ones 


t - 


ORIGINAL PAGE IS 
OF POOR QUALITY 




Figure 7 Quadri phase Demodulator Model 
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In channel 2, after demodulation by i/2 sinG^c^o}* we have 
r(t) t/2 sin (“ sc t+ ®o) ™ 7^2* ^2^ cos< i> " #7 d-j(t) S1 * n ^ + n '( t ) (25) 


where (f> = 0 q - and is the phase error s and n'(t) is an essentially white 
Gaussian noise process. Assuming that the bit transition probability is 
50%, the conditional bit error probability is given by {[assuming <j>(t) is 
essentially constant over max (T-,,T 2 )D » 



In order to evaluate (25), it is necessary to specify the phase 
error probability density. To do this, we assume, for lack of a better 
model, that p(tj>) (for a CW loop) is Tikhonov, i.e.. 


p(4>') = e ^! a i C 0 °(l) ,) ' » -it < <j>‘ < it (28) 

To convert to a loop having an S-curve that has a period of it 
radians (such as a Costas loop or a quadriphase with power unbalance of 
4 1), we let <j>' = 2<f> so that 


PU) 


exp(ctcos2<f.) , 
irr I Q (ot) 




(29) 


To determine the meaning of a in this modified Tikhonov density, let a be 
large so that we can use the asymptotic expression for (28), producing 
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P(+) = 



. {30} 


implying that the linearized phase error variance is given by 



(31) 


In the same manner, loops having S-curves periodic over ±ir/4 radians, 
such as balanced quadriphase tracking loops, have a density given by 


P(i>) 


2 exp(g cos4<t >) 
IT IqU) 



(32) 


and the linearized phase error variance satisfies 

2 „ 1 
°L 16a 


(33) 


The power ratio and data rate ratio for channels 1 and 2 

satisfies 




10 


(34) 


Using (29) to average over (26) produces the results shown m Figure 8. 
The parameter is the inverse of the linearized phase error variance 
so that 

a E = (35) 
a L 


We now consider the case where the high rate Manchester signal 
interferes with the low rate NRZ signal. 

From the received signal of (24), after demodulation by 
/2 cos(u) sc t+ Sq}, we have 
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(HR) ' (LR) 

r(t) t/2. cos^ sc t+e) = JP~ d 2 (t) sin<j> + JP ~ d^t) cos<j> + n"(t) (36) 

where n"(t) is an essentially white Gaussian noise process and is inde- 
pendent of n"(t). Again we notice that, if <j> a 0 identically, there is 
no cross-channel interference. 

When the transition probability is 50% after matched filtering, 
the detected signal is given by 


s 


d 



n 



* 


= jp7 t i cos< fr + 




,/P^T.j sm$ (37) 


where n^ denotes the number of ones in the first half of the Manchester bit 
and n 2 denotes the number of ones in the second half of the Manchester bit 
and n is the total number of ones in a Manchester bit (see Figure 6). It 
follows that the conditional average bit error rate is given by (n even) 


R(eU) = 


l l 
2 n/2 2 n/2 





(38) 

However, by direct numeration, it can be shown that a single summation form 
of (38) is given by 


P(EU) «• -L 
2 n 





(39) 


Since this is simpler, we shall use this. 

In our present case, n=10 and P 2 =4P 1 Evaluating (39) by 
averaging over <j>, i e , 



P(E|<}>) P(<i>) d$ 


PE 


(40) 
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yields the results of Figure 9. Note that these curves assume that <j>(t) 
is essentially constant over max(T-| »T 2 D • 

By comparing Figure 8 and 9, we see that the case of the low rate 
Manchester signal being interfered by the high rate and high power NRZ 
wave has more degradation than the case of the low rate Manchester inter- 
fering with the high rate NRZ signal. 

We note that both cross-channel effects and an imperfect carrier 
reference have been degradation sources for Figures 8 and 9. 

4.0 CHANNEL 3 NOISE INTERFERENCE EFFECTS ON CHANNELS 1 AND 2 

We assume m this section that noise dominates channel 3's input 
and that the amplifer is linear (neglect clipping) since the degradation 
is expected to be quite small anyway. The filter preceding the amplitude 
regulator is a three-pole Butterworth filter so that the noise spectral 
density shape is given by 



|H(f)] 2 


2 i + Cf/ffP 6 


(41) 


where f Q =5 MHz It is known that the lowpass filter (LPF) of Figure 3 has 

a rejection of 21.6 dB at 8.5 MHz. Furthermore, the 5 MHz LPF is down 

10 1 og [1 + (8 5/5)^] = 14 dB at 8.5 MHz so that the total noise is down by 

p 

about 35.6 dB. Since channel 3 is only roughly (22/12) = 5.3 dB larger, 

we conclude that the noise is down about 30 dB m the band of channels 1 

and 2. 

Consider this in greater detail We have, since channels 1 and 
-2 together deviate the carrier 12 MHz, that 

P CH1 + P CH2 = K ( 12 ^ ( 42 ^ 


where K is the constant relating signal voltage to frequency deviation. 
Now 


P CH1 _ 1 

P CH2 4 


(43) 


so that 



Average 
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or 


\Fcm - 5 - 37 K ■ VV% 


and 


V P CH2 " 10,7 K -^S 


2 +p n 2 


Now the new noise degradation in channels 1 and 2 due to 
nel 3 is given by 


DEGR-, 


\ + P N 3 L 
P N, 


and 


DEGR^ 


Pm +Pm L 
N 2 N 3 


where L is the reduction in channel 3 noise at 8.5 MHz. From (45) 
(46) , we have that 

p „ (5.37 )¥ 


and 


1 + 


PSi 

PNi 


(1Q.7) 2 K 2 


1 + 


52 

P N2 


From (49) and (50), and (47) and (48), we obtain 


( 44 ) 


(45) 


(46) 
chan- 

(47) 


(48) 
and 

(49) 


(50) 
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DEGRpn-v » lOlog <1 + 16.8L 1 + 


= 101og< 1 + 4.23L 1 + p 


Recall that L - -35.6 dB so that, at minimum noise conditions, i.e. 


S 

= 20 dB 


SNR, ■ -ff 1 - lOlogCB,!,] 


SNR 2 = 1j~ " 101 °9 C B 2 T 2^ 


Using (53) in (54) and (55) yields 


SNR-j 

= 5.9 dB 

(56) 

snr 2 

= 16 dB 

(57) 


Using (56) and (57) in (51) and (52), respectively, yields 


DEGR CH1 

= 0.1 

dB 


DEGR CH2 

= 0.2 

dB 

(58) 


We conclude that the degradation in channels 1 and 2 from the 
noise in channel 3 under worst-case conditions is essentially negligible. 
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5.0 UPLINK AND DOWNLINK NOISE COMBINING 

In this section, we briefly discuss the effect of uplink noise 
on the received signal. The signal flow in the bent-pipe mode goes from the 

user to the Shuttle, from the Shuttle to the TDRSS, and from the TDRSS to 
earth. 

The signal power to the total noise power is given by 


N o, +b o 2 


+ 


(59) 


where S denotes the signal power at the ground, No-j is the Shuttle noise 
spectral density received on earth, and No 2 is the receiver noise at the 
earth receiver. Equation (59) can be written in the form 




1 


+ 


M lV w °2 \ 

S/\S+N if 


( 60 ) 


where N-j/S is the noise power to signal power ratio at the Shuttle (K-band) 
transmitter. Notice that the right-hand side of (60) is expressed in terms 
of S/N 0l on the user to Shuttle link, the SNR at the Shuttle transmitter 
Cs/N-jD , and ratio of received downlink power to the earth receiver noise 
spectral density C(S+N] )/No 2 D * From (60) or (59), we see that, when S/Nq 2 
is large, (60) approaches S/No-j > as would be expected ([according to [10, 
pages 4r7 and 4-8][). 


'6.0 HARMONICS OF CHANNEL 3 TONES 

We now use Appendix III to evaluate the level of the harmonics 
of both tones from channel 3 in order to see if they will degrade channels 
2 and 3 BER performance. We accomplish this task by comparing the ampli- 
tude coefficient h Q10 with h Q7(r See Appendix III for details, particularly 
equation ( II 1-9) . The term h 07Q corresponds to the seventh harmonic of 
1.024 MHz and, as such, corresponds to an interfering tone at about 7.2 MHz 
which lies in the band of channel 2 when operating at 2 Mbps. 
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From Appendix IV, equation (46), we have for K= 0, m-] = 1 and 
m 2 =0 (see Appendix III, equation (9)} that 


*010 




2 . D 
p n + P 2 


-i 1/2 


, b /2 v p i 

V p i e £ f 0 £! (m+£)! 


£/2 

C 

3-2£-l 3 b 2 

p i 

1 F 1 

L 2 »2»2J 


r (l^i) 


( 61 ) 


where 


a. output clipping level of clipping amplifier 
d* input clipping level of clipping amplifier 

b: b = - 



, ■>„ = N o Be 

G: power gain of AGC in Appendix III, Figure 1 


A, 


A 2 

P 2 = , power in the 1.7 MHz tone 

P, 


p 1 : P 


1 


’ (% + p 2 ) 


input SNR of tone 1 regarding tone 2 as noise 


^F-j (a,b,z)- confluent hypergeometri c function 


p l 


{% + I P 2 ) 

P2= (^R 


To account for the AGC effects on the clipper amplifier (see 
Appendix III), we note that the RMS power into the AGC is given by 


P 1 + P 2 + *n 


RMS 


(62) 



If we let 
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then it is easy to show that 


(63) 


(64) 


2 SNR-j 

a l = 1 + SNR-j + SNR 2 

and 

2 SNR 2 

“2 = 1 + SNR-j + SNR 2 ( 66 ) 

2 2 

where SNR-j = P-j/o^ and SNRg = P ^ /Op • If we assume a representative tone 
power to noise power ratio of each tone to be -17 dB in the 5 MHz bandwidth, 
then. 


P 1 ~ p 2 = O' ^ 


(67) 


Set 


Pi + p 2 + a/ = 2 2 = 4 watts (68) 

and set the clipping level at the output of the clipper to be 4 volts 
(.two-sigma clipping level) so that a=d=4. We have by calculation 

2 2 

a-j - «2 - 0.183 

a n 2 = [1 - 2(0.0183)]2 2 = 3 854 watts (69) 


0.073 watts . 
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Evaluating (61) produces 

h Q10 = 0.195 (70) 

and the corresponding tone autocorrelation function is given by 


R-j (t) = 2(0. 195) 2 cos uj-j t (71) 

where 

= 1.024MHz (72) 

Now we will determine hg 70 We have from Appendix IV, equa- 
tion (46) that 


h 070 


b 



+ A 2 2 /2 
2 



7/2 



p / /2 



(73) 


Evaluating (73) after some algebra, we obtain 

h Q70 = 3 13 x 10~ 7 (74) 

so that 

7 2 

R-, (t) = 2C3 13 x 10"O cos (75) 

By comparing (71) and (75), we see that the seventh harmonic is down 116 dB 
from the original tone. We conclude that harmonics are totally negligible. 
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7.0 BANDPASS FILTER LOSSES: CHANNELS 1 AND 2 

In this section, we estimate the losses associated with filtering 
the quadriphase channels 1 and 2 signals and matched-filtering them. In 
Figure 10, a sketch indicates the block diagram including the bandpass fil- 
ter under consideration. The multipole channels 1 and 2 BPF has (upper 
left portion of Figure 10) a 3 dB bandwidth (specified) of about 8 MHz. 

For channel 1 at 192 kbps Manchester data (approximately 384 kbps NRZ 
data), we have that 


BT = 21 

and 

(BT)" 1 = 0.05 


(76) 

(77) 


By estimating the predetection filter losses compared in [9, Figure 7] 
with the transmission filter degradation from the transmission filter 
degradation of [9, Figure 6], we conclude that, for channel 1, the BER 
degradation is approximately 

DEGR-j = 0.06 dB 

However, the 0.5 dB ripple (filter specification) causes about 0 25 dB 
degradation based on [9, Figure 9]. The total degradation to channel 1 
is then given by 

DEGR 1 = 0.31 dB (78) 

We now consider channel 2 filter losses. Channel 2 has a 
2 Mbps NRZ waveform so that 


(BT)" 1 = 0.25 (79) 

so that the BER degradation, using [9, Figure 6] (seven-pole curve) and 
[9, Figure 7] to relate predetection and transmission filter losses, is 
about 0.4 dB. Additionally, the 0.5 dB ripple in the filter bandwidth 
could produce about 0.25 dB degradation. Consequently, the total chan- 
nel 2 bandpass filter loss is about 



Channel 1 
Channel 2 


Channel 3 


Figure 10 
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DEGR 2 = 0.65 dB (80) 

Now we estimate the BER degradation when the additional BPF and 
LPF shown in Figure 10 are used to reduce channel 3 tone interference of 
channels 1 and 2. These filters are identical in filtering characteris- 
tics to that of transmitter bandpass and lowpass filters of Figure 2. 

From [9, Figure 6], it is estimated that the cascade of BPF's increases 
the BER degradation less than 0.05 dB for channel 1 and less than 0.2 dB 
for channel 2. We conclude that the overall BPF BER degradations are 
given by 


DEGR ] ' 

= 0.36 dB 

(81) 

degr 2 ' 

= 0.85 dB . 

(82) 


8.0 RF-FM MODULATION/DEMODULATION LOSSES 

The mam thrust of this analysis has been at baseband; however, 

the FM modulation/demodulation process is not perfect, so BER degradation 

-5 

at PE=lxlO has been assumed to be 

DEGR = 0.5 dB (83) 

This degradation estimate includes nonlinearities as well as 
the effect of thermal noise in the FM demodulation process. 

9.0 DEGRADATION OF CHANNEL 2 DUE TO CHANNEL 3 TONES 

From Figure 10, we see that the two tones of channel 3 will be 

"seen" by the integrate-and-dump matched-filter detector for channels 1 

and 2. 

Therefore, if the additional filters are not included as shown 
in Figure 10, the cross-channel interference of channel 3 onto channel 2 
will not be negligible. In Figure 11, the two tones of Figure 3 are 
illustrated along with the channel 2 matched-filter response. Since chan- 
nel 2 has a much higher data rate, its response is affected more by the 
channel 3 signals modeled m the report as two unmodulated tones. The 
channel 2 matched-filter response is given by (in both time and frequency) 




Figure 11. Response of Channel 


Channel 2 
Matched-Filter 



the Two Channel 3 Interfering Tones 
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with corresponding frequency response 

H(f) = sinUfT) g-jiuT/2 

irfT 
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(84) 


y(t) 


Define the channel 2 signal -plus-received noise as 

& A d(t) sin togt + /L n g (t) sin cogt + v4T n (t) cos Wgt 

+ /l A t sin to-jt + /l A t sin (85) 

i i 2 2 


where n c (t) and n s (t) are independent, bandlimited, baseband, Gaussian 
noise processes. When this signal is coherently demodulated by v^sinugt, 
there remains 


e D (t) 


A d(t) + n (t) + A t 
s i-, 






cos 




( 86 ) 


where 

f Q = 8.5 MHz 

f ] = 1.024 MHz 

f 2 = 1.7 MHz (87) 


by 


We see that the interfering signals and noise have powers given 


P T = VarCn 3 + P T + P| + P T + P' 

i 5 ! 2 l 2 


( 88 ) 


where 
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Var&g = \ 


T C T 


0 0 


dt du 


VarCn s } ■ £ 


with Ng/2 the two-sided noise spectral density and Pj-j and corre- 
sponding to the Ai-j interference terms in (86) and Pj 2 and Pf 2 being the 
Aj 2 interference terms in (86). 

We now evaluate the channel 2 interference terms. We have 


sin 2 f7.48TTxl0 6 /2x 10 8 ^ 
(7.48tt/2) 2 


P 1.024 0-0039 P-| Q24 


sm (9. 52?r/2) p = 0 0021 P 

(9.52V2) 2 P l- re4 T- 024 

sin 2 (6.8ir/2) p 0.0079 P, 7 

(6.8u/2) 2 1 ' l * / 

s in 2 (l 0 ; - 2ir/2) p = 0.00037 P, 7 
(1 0. 2ir/2) 2 l *' 1#/ 


In order to compute the total degradation, we must evaluate the 
tone power levels Assume that channel 3 is at maximum SNR (no noise) so 
that the sum of both tones modulate the carrier to 22 MHz peak-to-peak. 
Therefore, for channel 3, 


2[A 1 +A 2 3 = K 22 x 10^ 


where A-j and A 2 are the tone peak amplitudes of channel 2. Since 


( 93 ) 
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g 

^jT- = A ] = 11 K - ^ 10 ' = y- K' 00 = KxlO 6 ] (94) 
So 


P T = 30.25(10 ) 2 (95) 

'3 


and the individual tone powers are 

P T = P T = 15.13 (K') 2 (96) 

1 31 1 32 


Now in channels 1 and 2, we have, since the filtered quadri- 
phase is in quadrature, that 

2 ^Aj 2 ^ A Q 2 = 12K x 10 6 (97) 


where Aj is the peak channel 1 subcarrier amplitude and A ^ is the peak 
channel 2 subcarrier amplitude. Also, 

P T Ii(! - I *Q* • fl Q ’ 2fl I ' 98 > 

Solving (97) and (98), we obtain 


Aq = 5 37 K’ 

Aj = 2.68 K 1 

P T = 18.1 (K 1 ) 2 

I »Q 


CK* = K x 10 6 3 


(99) 


Now, if the SNR in channel 2 is 10 dB, then the noise power in the Q chan 
nel is given by 
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(14.4)(K*) 2 


Therefore 9 


C N o B ^ CH 2 = 1 - 3 ( K ') 2 


( 100 ) 


( 101 ) 


where B is the equivalent 4 MHz lowpass process from the 8 MHz bandpass 
process. It follows that 



1.3(K*) 2 

B 


(102) 


Therefore, (89) 


can be written using (90), (91) and (102) as (j=l/2x!0 6 ) 


( p X N ) » 2 Tb 1.3(K') 2 + 0.0039 (15. 13)(K 1 ) 2 + 0.0021 (15 13)(K‘) 2 
CH2 

+ 0.0079 (15.3)(K‘ ) 2 + 0.00037 (15.3)(K’) 2 (103) 

This can be factored as (b ~ 4 x 10 6 ) 

P T = 0.325(K') 2 [1 + 0.66] (104) 


We conclude that the degradation is given by 


DEGR ch2 = 2.2 dB 


(105) 


Now we consider the same calculation utilizing the benefit of 
the BPF and LPF in Figure 10 at the receiver to chop out the I&D filter 
response near the channel 3 tones Figure 12 illustrates this point. 
Equation (103) can now be written as 


P T n = 2TB H. 3 )^') 2 + 0-0039 (15.13) (K 1 ) 2 (^-f-) 2 

+ 0.0021 (15 3)(K') 2 (^~) 2 + 0.0079 (15.13) (K 1 ) 2 (^^) 2 
+ 0.00037 (15 13)(K’) 2 


(1061 
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where the factor 

/ 0 . 02\ 2 
\ 0.8 / 

is the reduction in filter voltage response (squared) due to the presence 
of the ADL BPL. 

Evaluating and factoring as in (104), we obtain 

P T = 2.6(K’) 2 [1 + 0.0004] (107) 

*N 


We conclude that the degradation due to the channel 3 tones is less than 
0.01 dB for channel 2 when the transmitter BPF is also used on the receiver. 


10.0 CHANNEL 3 DEGRADATION DUE TO CHANNEL 2 NOISE 

This source of degradation is based on the fact that channel 2 
noise can get into channel 3, thereby increasing the noise level. The 
transmitter BPL, however, limits the noise interference as well as the 
5 MHz lowpass filters into the channel 2 limiter. It should be noted that 
the additional BPL and LPF at the receiver does not affect this phenomena. 

From the ADL filter response data, the filter reduces the chan- 
nel 3 noise by 


L = 



-32.1 dB 


(108) 


When channels 2 and 3 are near the same SNR's, we see that this additional 
noise will be essentially 32 dB down, plus or minus a few decibels due to 
the different FM drive levels. It follows that this noise has negligible 
effects on channel 3 performance 

Now consider the degradation to channel 1. Equation (88) still 
applies, but now Var(n) denotes the variance of channel 1 noise. The 
response of the Manchester matched filter is given by 


sm Z (nfT/2) „-ja>T/2 
— ' 0 


H(f) 


(irfT/2) 


(109) 
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Equation (91) becomes, for this case. 


p = sin 4 fir7.48xl0 D /0.4xl0 b ^ 
T 1 (7.48ir/0.4) 2 


6 


P 1.024 " °* 00012 p -[ 


sisi ( - ^y° ; 4) P •. 0 . 000021 P 

(9.52ir/0.4r ,U24 1,024 


sin 4 (6.8ir/0.4) 
(6.8 tt/ 0.4) 


2 P T.7 " 0 


P L 


p- _ = 0.00016 p- 7 
(10.2ir/0.4)^ 1 * 7 


From (99) and assuming that channel 1 is at 10 dB, the 
power is then given by 


jj (3 59) (K 1 ) 2 

Therefore, 

C n 0 b 3cH1 = 0.33(K* ) 2 

It follows that 



0.33 

B 


( 10 ) 


2 


Hence, using (88), (90), (110) and (112), 

( P T ) = 2TB (0-33)(K* ) 2 + 0.00012(15 

' Nj/ chi 

+ 0 0001 6(15. 13) (K* ) 2 


we have 

13) (K 1 ) 2 + 0.000021(15. 
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.024 


(no) 

noise 

( 111 ) 

(112) 

(113) 

13){K') 2 

(114) 


This can be factored as 
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0.0825(K') 2 [1 + 0.055] 


We conclude from (115) that the degradation is given by 


DEGRqj^i - 0.23 dB . 


015 ) 


016 ) 


11.0 BER E /N n DEGRADATION SUMMARY 

D D 

In this section, we summarized the Eg/Ng BER degradation for 
channels 1, 2 and 3. The link budget for channel 1 is given in Table 3 
and the link budget for channel 2 is given in Table 4. 


Table 3. Channel 1 BER Degradation 


Component 

No Receiver 
Filters 
(dB) 

Receiver 

Filters 

(dB) 

1. Imperfect carrier reference including 
crosstalk (assumes cr, = 4.5° @ BER - 
1x10-5) L 

0.2 

i 

0.2 

2 Channel 3 noise interference 

0.1 

1 

0.1 

3. Harmonics of channel 3 tones 

0.0 

0.0 

4. Lowpass filter losses 

0.31 

0.36 

5. RF-FM modulation/demodulation losses 

0.5 

0.5 

6. Matched-filter response of channel 1 
to -channel 3 tones 

0.23 

0.0 

7. Channel 3 degradation due to channel 2 
noise 

0.1 

0.1 

Totals 

1.44 

1.26 
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Table 4. Channel 2 BER Degradation 


Component 

No 

Receiver 

Filters 

(dB) 

Receiver 

Filters 

(dB) 

1. Imperfect carrier reference including 
crosstalk (assumes a. = 4.5° § BER = 
1x10-5) L 

0.7 

0.7 

2. Channel 3 noise interference 

0.2 

0.2 

3. Harmonics of channel 3 tones 

0.0 

0.0 

4. Bandpass filter losses 

0.65 

0.85 

5. RF/FM modulation/demodulation losses 

0.5 

0.5 

6. Matched-filter response of channel 2 to 
channel 3 tones 

2.2 

0.01 

7. Channel 3 degradation due to channel 2 
noise 

0 1 

0.1 

Totals 

4.35 

2.36 


These tables do not include bit synchronizer degradations or 
relay losses or take into account that the modulation BPF centered at 
8.5 MHz is not narrowband, so that additional losses may be present. In 
addition, they do not include the C/Ng loss out of the limiters in chan- 
nels 1 and 2 (see Appendix II). Also, numerous assumptions on SNR were 
made to obtain the estimates, with most of them based on worst-case con- 
ditions ' Therefore, the two BER degradation tables represent, for the 
most part, a worst-case estimate for each entry. 

Channel 3 BER degradation was found to be negligible due to 
additional thermal noise from channels 1 and 2 

12.0 CONCLUSIONS 

Estimates of BER degradation have been obtained, in most cases, 
by assuming worst-case conditions. One exception is the four-phase loop 
which demodulates the quadri phase signal of channels 1 and 2. There an 
RMS phase error of 4.5° was assumed as a reasonable guess. 
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It is seen that the additional filters of Figure 10 reduced 
BER degradation in decibels almost in half for channel 2 and about 10% 
for channel 1. 
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APPENDIX I 

BASEBAND DIGITAL SIGNAL-PLUS-NOISE THROUGH A HARD LIMITER 

In this appendix, we determine the output spectral density of 
a baseband digital signal plus band! invited Gaussian noise. This calcula- 
tion allows one to determine the output C/Nq value as a function of the 
input C/Nq value. The derivation of the autocorrelation function follows 
that of Painter [1] with minor corrections in the final results. The 
model is shown in Figure 1 below. 


Am^ (t) 


+ n(t> 


LIM 


d-j (t) 


Figure 1-1 . Baseband Limiter Model 


where m-j (t) is an NRZ or Manchester binary-valued digital waveform, n(t) 
is bandlimited, white Gaussian noise, and d-j(t) is the output binary-valued 
waveform. Let 


x(t) = Ad(t) + n(t) 
y(t) = sgn[x(t)] 


(I-D 

(1-2) 


The first- and second-order density functions are given by 

(1-3) 


P x & t .t3 * j= exp 

/2i r 


" A C x t^ t 3 

2a 


P x & t .x T D = 


1 

-i rw ri 

. i , 

,Cx T -ii T ] 2 2pCx T -u T 3Cx t -p t 3 i (x t -w t ) 2 ) 

CAP 

ZwN^-p 2 

20-p 2 D 

2 2 2 ( 
a a a 1 


(1-4) 




1-2 


where 


or 


x t = 
x = 

T 

V t = 
P T = 
a 2 . 
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x(t) 

x(t+x) 

E[x(t)] « Ad (t) = s(t) 
E[x(t+x)] = Ad (t+x) = s ( t+x ) 

E [ & t _,1 tP 2 ] = E^-tQ*] 

1 


P = 


2 


iUx) 

E fn(t+r)n(t)] = - n 9 


(1-5) 


First let us compute the mean value of our process. We will 
temporarily drop the time subscript. 


v = 


/ 


sgn(x) p x (x) dx 


( 1 - 6 ) 


which can be expressed as 


V = 

/*" 1 

# — exp 

- \ (x-s) 2 

f 0 i 

dx - / — — exp 

- -V (x-s) 2 

J 

q o 

L 2a 2 J 

J o/2tt 

— CO 

2a 2 


dx (1-7) 


where s = s(t) = Ad(t). Letting - v produces 


2_ 

/ir 


/ 


JZa 
2 / Jlo 


e v dv = erf 


s(t) 

i/2 a 


( 1 - 8 ) 


Now let us determine the output autocorrelation function defined 


by 


/ / sgn[x t ] sgn[x T ] p x Cx t5 X T )dx t dx T (1-9) 


— 00 — oo 


R y (t+x,t) = E[y(t+x)y(t)] = 



Equation (9) can be written 
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1-3 


00 -CO Q Q 

V ttT ' T) ■ J J PxCV\3 dX t dX T * f f P X &f X T^ dX t dX T 


0 0 


f f PxC» t .X T 3 *< t dx T - / / P x Cx t ,x T }dx t dx T (1-10) 


Letting x^. - = u<r, x - v = vo produces 


y. v 
t r 


R y (t+T,T) 


/*” f m c /*u 

Jv t J v g ^ u,v,p ^ dudv + J J g( u ’V,p) 


dudv 


-00 -00 


y T y t 

- J ° g{u,v,p) dudv - J ° f\i 9 ( u ’ Vs 


p)dudv 
( 1 - 11 ) 


where g(u,v,p) is the density function of the random variables u and v 

induced from the variables x. and x . 

t r 


Let 



k = 



a 


then 


R yy (T,t) = L(-h,-k, P ) + L(h,k, P ) - L(-h,k,-p) - L(h,-k,-p) (1-12) 
so using [2] 


L(h,k,p) = Q(h)Q(k) + f p n+l 

n=0 (n+l)l 


(1-13) 



and 
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Q(h) 



1-4 


(1-14) 


we arrive at 


“ n, t,\,(n), ,x n+1 

R w (t,r) = Q(-h)Q(-k) + I ^~- h)z ^ 

yy n=0 (n+1)! 

+ q(h ) Q ( k ) + f A^tk^ 

n=0 (n+1)! 

. Q( . h)Q(k ) . f ^(-hiz^Vu -pr 1 

0 (n+1 ) ! 


where 


- Qth)Q(-k) - 


f z (n) (h)z (n) (-k)(-o) nt1 
n=0 (n+1)! 


z (n) (x) 


dx n 



(1-15) 


(1-16) 


Now 


R y (t,t] » CQ(h)-Q(-h)][q(k)-q(-k)] 


+ I -^-|z (n) (-h)z n (-k) + z (n) (h)z (n) (k) 
n=0 (n+1)! ( 


- (-D 


n+1 


: (n) (-h)z Cn) (k) + z (n) (h)z (n) (-k) 


(1-17) 



Now note that 
Also note that 

So 

Ry(t,x) - erf | 

+ 4 

where 5 denotes 
indices yields 

R y {t,x) = erf 
+ 4 
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Q(-x)-Q(x) = 


z^ n ^(x) = z n (-x) , n even 

z^(x) = -z n {-x) , n odd 


'h_) erf (x\ + 5 >, Z (n) (h)z (n) (k)p n+1 
,/2/ \/27 n=0 (n+1 ) 1 


^ z( n >(h)z (n) (k)p n+1 
n=l (n+1)! 


sums over n odd and denotes suras over n even. 



00 


1 

m=0 


P 2m+1 .(2m) 
(2m+l)« 



1-5 

( 1 - 18 ) 

(1-19) 

( 1 - 20 ) 

Changing 

, 2m+2 

M — 

(2m+2) 1 
( 1 - 21 ) 


We have, since s(t) = Ad(t), 
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1-6 


Ry(t,T) 


= d(t)d(t+x) erf^ 


+ 4 d(t+x) d(t) 2 
m=0 


■W) Lie 

N 


2 ^(Tn 2 ^ 2 


'(2m+2) ! 


+ 4 l 

m=0 


,(2m )LjS 


2 r R n (x)“ r2m+1 


( 2m+1 ) ! 


(1-22} 


where we have let 


A 


2 ~ N 


(1-23) 


Now consider d(t) ; 


d(t) = I d p(t-iT) 

•?= 1 


1 =-' 


where p(t) = l,te(0,T) and d^ = ±1 with equal likelihood. Hence, the time 
and ensemble average of d(t)d(t+x} yields 


or 


R d (x) = <(d(t)d(t+x})> =/ f p(t-iT) p(t+x-iT)' 


vi --<*> 


R d (x) 



NRZ waveforms , 


so that we have 


C- SL 
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= R d (x) err 


+ 4R d (x) l 

m=0 


,(2m+l ) 



£ 

N 


*n (t n 


2m+2 


'(2m+2)! 


+ 4 I 
m=0 


rt2m)A/I 



“ ,2 ^(tJ! 


N, 


(2m+l ) 


^(2m+l ) ! 


(1-24) 


The first term is the SxS term, second is the SxN terms, and the last term 
is the NxN term. Notice at high SNR 


R y (x) s R d (x) 


(1-25) 


At low SNR 



lim erf Vtt - 0 


S n 


1 im 

1+0 L 



2 

if" 

(2m)! 

2tt \ 

2/ 

m! 


lim 
S n 
fT° 


.2m+lL IS 


n2 


- 4 - 0 


(1-26) 


Therefore, (1-24) converges at S/N=0 to 


lim R (t) 
N’ V ° 


4 I 
m=0 


iiyf 

2m n 

2 1 


v*)\v 

_m'_ 

(2m+T)T 

2 


2m+l 


(1-27) 


or 


R y ( T) 


arc sin 


ry*Ji 


(1-28) 


which is a well-known result. It is easy to show that R d (-r) is given by 
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i -JiL 

1x| < T j 


1 

> NRZ 


0 

M>t J 

(1-29) 


R hW = 


i -11: 


T 

< 2 




} < 


Manchester 


• < T 


(1-30) 


Therefore, at the input to the FM summer due to channels 1 and 2, we have 

s(t) = A ~ d^(t) cose sln G° sc iO + A ~ d-j (t) sine COS G J|1 SC '0 (1-31) 
where d 2 (t) has the spectral density given by 



R„(t) 



(1-32) 


with^(-) the Fourier transform and where d-j(t) has the spectral density 
given by 


jf d-j (f) = ^(f) 



R n (x) 




(1-33) 




1-9 


and where >e^(t) denotes the data spectrum of that respective channel. 

Evaluation of the SxN and NxN spectral densities is made in 
Appendix II. 
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APPENDIX II 

SPECTRAL DENSITY OF SxN AND NxN TERMS OUT OF THE LIMITER 

The purpose of this appendix is to obtain the analytical form 
of the NxN and SxN terms at the output of the limiter when the input is a 
baseband signal plus bandlimited Gaussian noise. We follow Helgeson’s [2] 
work in this section and slightly extend his results to cover fifth-order 
NxN terms. We model the lowpass filter into the limiter as an ideal LPF 
of the form illustrated in Figure II-l below. 


PAGE IS 
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(B = 5 MHz) 


► f 


-BOB 

Figure II-l. Baseband Noise Spectral Density Model 


We therefore model the limiter input noise spectral density by 



f [u(f+B) 


u(f-B)] 


(II-l) 


where denotes the original process and n^ is the kth-order convolution. 
Now, using Figure II-2, we have 


n 2 (f) 


J */ n -jC f i3 J n C f -f]D df -| 

— 00 

(r) / [C f ! +B D - uCf r B] 


uCf-f-j+B] - uCf-j -f-j -Bl) 


df 


(H-2) 
(II -3) 


Evaluating 
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1 1 —2 


V f > 

{c — 


5 MHz 


-B 


r* “ 


Vn-j & - f-| D 

1 

I 

1 

i 

i 

i 


! 

1 

(// // 
K/// 


, i 

fss/ 


f 

- B (> 

f + B 


Figure 1 1 -2 Spectra for Determining 


so that 



2B-f f > 0 

2B+f f < 0 


^ 2 { f ) = (t) [2B " l f l ] I f l <2B 


= 0 


|f| > 2B 


- Now determine 



Using Figure II-3, we have 



/ - f.p df ] 


So that 


-*4, 


(f) 


r » 

-CO 


-I^D (uCf-fT + B} + 



(II -4 ) 


(II-5) 


(II-6) 


(H-7) 



1 1 -3 
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Vn 2 Cf-,3 



Figure I 1-3. Spectra for Determining ^n 3 Cf-|D 


In the region 0 < f < B, we then have 


n*(f) r° 


/ 0 /-f+B 

[2B- |f-, Dd^ + / C2B-|f ] Ddf 1 

f-B 0 


(II -8) 


Since >/n 3 (f) is an even function, we obtain 


Vn,(f) = (~) C3B 2 - |f! 2 ] 0 < Jfj < B (II-9) 


Now in the region B < f < 3B, we have 


j/n-(f) 


/*2B r 

~ j [2B - |f| 

f-B 


df 


Since t/nJf) is an even function of f, we have 


yn~(f) q 9 f 2 

1ZTT = |B 2 -3B|f| + ^ 


B < { f | < 3B (11-10) 


3 



To summarize, we obtain 
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*/n 3 ( f ) = (^) = C3B 2 - f 2 D 


,(f) - 


- (r) (f ® 2 - 3 B l f l + t ) 


0 < If | « B 


B < f < 3B (II-H) 


Note that the two equations are equal at f = B. Now consider 
Figure I I -4 illustrates the spectra involved in the convolution. 


J*t^) = J J n s ~ df ] 



Figure I 1-4. Spectra Used for the Fourth-Order Convolution 
of the Noise Spectra 


Consider the case 0 < f < 2B. We have 



(11-13) 


or, using symmetry. 
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In the region 2B < f < 4B, we have that 


0 < [f[ < 2B (11-14) 



/ 3B 

. 


| B 2 - 3B|f|+ ~ 


df 


2B < f < 4B (11-15) 


After evaluating, we obtain 



^ B 3 - 2Bf 2 + J- jf| 3 


0 < |f| < 2B (11-16) 



-y- B 3 - 8B 2 |f j + 2Bj f j 2 - 



2B < |f } < 4B (11-17) 


Now compute 


A( f) 


Figure I 1-5 illustrates the spectra involved. Now 



Figure I I- 5 


Spectra Used for the Fifth-Order Convolution 
(0 < f < B) 
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Evaluating for 0 < f < B, we have 



so that, after some algebra for 0 < f < B, we have 



115 R 4 

TT B 


- ~ B 2 f 2 



(11-18) 


(11-19) 


( 11 - 20 ) 


Now consider the region B < f < 3B illustrated in Figure I I -6 



Figure 1 1 - 6 Spectra Used for tne Fifth-Order Convolution 
(B < f < 3B) 




y- B 3 - 2Bf 2 + l 


, \ rf+B 

’) df + / 

*2B 


3 2 

B J - 8B^ 


f + 2Bf 



(11-21) 


Evaluating (11-21), we obtain by symmetry 
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= B 4 + | B 3 if I - 5B 2 f Z + | B jf | 3 



B<|f|<3B (11-22) 


Now consider the case when 3B < f < 5B is illustrated in Figure II-7. 



Figure 1 1 - 7 Spectra Used for the Fifth-Order Convolution 
(3B < f < 5B) 


We have 




—oo 


n 

f-B 


B J - SB 2 [ f j + 2B | f | 


2 lfL 

6 


df 


(11-23) 


(11-24) 


After some tedious algebra, we obtain by symmetry 



1875 B 4 _ 1|5 B 3 f + 25 R 2 f 2 - | Bf 3 + ^ f 4 3B< |f| <5B (11-25) 




In conclusion, for the fifth-order convolution, we have 


II-8 


n 5 ( f ) 


2 

K(f) 


B 4 + |- B 3 f - 5B 2 f 2 4- 4 Bf 3 
6 3 3 


115 „4 5 R 2 f 2 r 

TT B " 2 B f + T 




0 < [ f 1 < 3B 


0 < [fl < B 


/Ac(f) 


5^ . 1875 b 4 _ 125 ^ x 25 R 2^2 


,5 72 


f + Bf - f Bf° + r 3B- < If I < 5B (11-26) 


From equation (1-22), we have that the spectral density out of 
the limiter is given by 


4 lf] - 





R dW R n‘ 






n 2 


1 

10 


Jr\Af) 


(11-27) 


-where x/^(f),j^ 3 (f) and^^f) are given in equations (II-l), (11-11 ) and 
(11-26), respectively. 

The second and third terms are SxN terms and are yet to be eval- 
uated Consider the spectral density at f * 0. 


*4xN^ \ 

a 





(11-28) 


or 



2 
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>} 


+ 




,(D 


n 2 


n 0 2 /b 


+ 0.056 


.( 3 ) 




n 0 4 /b 


(11-29) 


(11-30) 


The output signal power first term of equation (11-27) to NxN 
and SxN spectral density is evaluated in Table II-l using equations (11-27), 
(II-l), (11-11), (11-26) and (11-30). 


Table II-l. Input and Output SNR Comparisons 


©IN 

(dB) 

/C \ 

K 

' u IN 
(dB-Hz) 

iw) 

X U 'OUT 
(dB) 

/ c \ 

N - r ) 

Vo / 

\ u / 0UT 

dB-Hz 

10 

77 

27.3 

94.3 

5 

72 

10.4 

77.4 

0 

67 

- 1.3 

65 7 

- 5 

62 

- 6.2 

60.8 

-10 

57 

-12.15 

54.85 

-15 

52 

-17.4 

49.6 


Notice that the output C/Nq ratio is greater than the input 
C/Nq ratio at 5 dB SNR^ and higher. Due to this apparent improvement 
(which parallels the bandpass limiter I/O ratio improvement [1]}, we con- 
clude that knowing C/Ng is not sufficient to estimate BER degradation. 
Consequently, we leave this component out in the BER degradation tables. 
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TWO TONES PLUS THERMAL NOISE INTO A CLIPPER 
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The clipper model is shown in Figure III-l. In this section, 
we follow the work of Davenport and Root [1], Helgeson [2], Shaft [3] and 
Shaft [4]. We extend the analysis to the two-tone case using an approxi- 
mation due to M. Gyi [5] and an extension of that approximation due to 
the author of this report [6], 

The input to the clipper is modeled as the sum of two tones and 
band! united thermal noise: 

1.024MHz 1.7MHz 

x(t) = cos ( a) sc 1 t+ e l) + V^2 cos (“sc 2 t+ e 2) + ( nI -D 


where P-j and P 2 are the respective powers of the two tones. We use the 
indirect approach [1] and form the autocorrelation function 


RyCt-j »t 2 3 



(II 1—2 ) 


From [1, section 13-3], we can write that 


R yC t T t 2-) ^^2 y* f C w i3 daj i j da) 2 M x^“l“2^ ( In * 3 ) 


where M x C“i >“ 2 3 is the joint characteristic function of xt-j and xt 2 * 
.Now for independent noise and signal processes, we have 


M x fr>i ^ »co 2 3 M n C Wl ,0) 2 } M s s w 2 3 M s (w-, ,to 2 }.. . (Ill -4) 


From Shaft [3], we can write 



LPF 

(5 MHZ) 



y(t) 


g(x) 



(a) Block Diagram 


(b) Nonlinearity Transfer Function Model 


o o 

-n 33 
-0 G5 

o a 
o 

X r 3 

<° IS 
c ^ 

3? 6> 

r ^ 


Figure III-l Clipper Model Used in the Analysis 


i — I 


J 
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M xh ^ 


where 


CO 00 


I I ... I 


e e 


m l =0 m 2 =0 

m^O m l m 2 

m n 

I—. (MO i 

nii v I i- / m 0 2 2/ m 

l 2 n 

GW 

r* 

O 

(/) 

3 

</> 

O 

u 

■ cos^co^t) 

. . . cos 

e = 1 

m.. = 0 


m i 



= 2 

nt. f 0 



9 





Now I mi -(x) is the modified Bessel function of the first kind and of 
n^. Hence, we obtain 


^yC^i >t 2 3 


if/ f f 


oo 00 


f&O Z ' S e m "s_ I CAiuO I 

2 m^=0 m n =0 m l m 2 m l 1 1 \ 


k =0 


This can be written as 


CO CO oo 


••• I ni 1 CA 1 M 2 D ... I^CVjD • exp 
“ R n k (T)u 1 k M„ k 

x 2 


n / 2 , 2 


2 i o)-| 


! ) 


k! 


* cos m 


-. 0 ) t\ ** cos/m (o t\ 
1 SC 1 / \ n sc n ) 




^yC^i Z Z Z " "r. | — e e * • • e h £ , 

y k =0 m n =0 m „=0 k m i m 2 % kn, l m 2 m - 


"1 u n 


cos 


( m l“sc 1 ' r ) cos ( m 2“sc 2 j •" cos ( m n“sc n T ) 


m-3 


(II 1—5 ) 


order 


CA n «-[D 


(III-6) 


( 1 1 1 ~7 ) 


where 
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Wm.,. . .nr 
l i n 


m f fM “ k 'm/M 


2 2 
a a) 


• e 


III-4 

dto 

(in-8) 


where c is an appropriately chosen contour in the complex plane. Shortly, 
we will discuss this contour more precisely. 

An expansion for two signals is given by 


Ry * T ' k=M,... nt^=0 nijr=0 ^ km l m 2 


cos ( m 2“sc 2 T ) 


m-j+m 2 =0,2,. . . 


+ ^ S ^ m cos 

k=0,2,... m^O m 2 =0 111 1 m 2 km l n 2 


m l +m 2 = l ’ 3 > 


( [ "l“sc 1 T ) cds ("2“sc 2 t ) 

( III— 9) 


Now consider the clipper amplifier shown in Figure III-l(b) and 
the expression in equation form in equation (III-10), 


-a 


g(x) = (a/d)x 


x < -d 
| x [ < d 
x > d 


(III-10) 


From ( II 1-8) , the h km m are given by 


km ^mg 


2 2 
0 ( 1 ) 


2ttj /W I m 2 ^ A 2°^ e ^ 2 


dw (III-ll) 


where A-| = /2P^ , A 2 = /2P 2 , the A^'s are peak voltages and a = N Q BG, with G 
the AGC power gain. 
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1 1 1 -5 


Let 


g.(x) = 


g(x) x > 0 

0 x < 0 


; g_(x) = < 


0 x > 0 

g(x) x < 0 


and also let 


fU) = f + (w) + f_(aj) 


with 


f>> - f 


g (x) e"“ x dx ; f _(w) = / g_(x) e _a)X dx 


/ 


For an odd transfer characteristic (such as our clipper), we have 


f n U) = 


f Q (a)} - -<J>(-0)) 


Hence, from equation (III-ll), we have 


2 2 
a oi 


4, (“)“ k e 2 d “ 


2 2 
O Oi 


/ u tu 

e 2 dffl 


(in- 12 ) 


(in-13) 


(III-14) 


(III-15) 


(III-16) 


where c + and c_ are shown in Figure III-2. Let -<»> = w' in the second term 
of ( II I -1 6) so that we get 



m-6 



2 2 
0 to 


km^m 2 


ST j I m ] 1 m 2 ^> e 2 d “ 


2 2 
0 0 ) 


1 /* r, ITU — y 

- EFT 7 (-1) H) T * &>'A n X e V) da) 1 


*V' m 2 ^ 2^ 


(III -17) 


k+m-j+nig 


So when (-1) =1, i.e , k+m-j+m 2 is even, then 




(k + eyen) 

2 2 
0 £i) 

h km 1 ni 2 ' K7 / 'm^V 1 m 2 ^ e ? d “ 

C + 

(k + + m 2 odd) (III-18) 


We now compute <j>((d). We have that 
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-f 


4>(w) = / g(x) e “ x dx 


•A 


/ 


x e'“* dx + S a e dx 


or finally. 


♦w ■ 


Notice that we can write 


$(w) = limf 
S+G 


e -6a)_ e "do)~| 


0) 


So 


I 


( 1 ) 


km -j nig 


2 2 

/ -So] a M 

4- e y - 0) k I QoA-,3 C w ^3 e ^ dco 
d 2 m 1 L m 2 2 

c + 


> 


( 2 ) 


2 2 s 
CT„ (0 


" 2uj y d 6 2 “ k l mS^ exp \ n 2 

y to I Z 


dco 


and 


for k + m-j + m 2 odd 


h, m = 0 for k + m, +m 0 even, 
km i nig I c 

Note that h. = 1^- I^\ Rotate axis tt/ 2 radians clockwise b\, 
letting w= iz so that (see Figure I I I -3 for the new contour) 


(III-19) 

(III-20) 


(III -21 ) 


(II 1-22) 


(II 1—23) 



m-8 


OF 




Figure III-3. The Contour c + After Rotation 



IT 1 



a e 


-iSz 



(1z)k 


2 2 
a z 


idz 


or 



Let 



(iz) 




(z) 


We finally obtain that 


I 


(D 



Let 



2,2 

CT Z 


-2 k-2 -i5z , 

z e J 


m n 


J m CA ? z) e 


dz 



(x) 



(iz/2) 2£ 

£* r (nii+£+l3 


( I II -24 ) 


(III -25 ) 


(III-26) 


(iri-27) 


(nr-28) 
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III -9 


so that 


( 1 ) _ 


/A \ m- i + 2£ 4 

.k+m-|-Hn2-2+2£ ( rt ll /. / 9 \ 

4- 1 ^ f z k+2+m ’ +21! 0 n CA 2 z3exp -ISZ-a 2 z \ d2 

1,1:1 £=0 £!rCm,+W) J , "2 2 V " 2 / 


(III -29) 


To evaluate this expression with n^O, we approximate JqQ^z} by 

2 


J 0 CA 2 z} = e 


[V z 2 

\Z 2 ) 


(III -30} 


to get 


A, 


\irW-2£ 


T (l) A. f 

1 ird Z 


1 k+m 1 +m 2 -2+2£p_ 


£=0 £'rCm-| + Jt+0 



From (III-31) and (III-28), we have 


Ck+m- | +m ? -2+2£3 /A 1 \ m l +2£ -k+2-m 1 -2£ 

/ ,, a « l l-sW iri 

T (1 J _ _a_ V ______ V 2 / 

1 " ird z 


£=0 


£irCm 1 +£+l3l — 


2 A 2 
ct + — o 


2\ k ' 1+ni l +2£ 


1 F 1 



k+m-|+2£ ^ _ fi 2 


3 - k ~ m-j - 2£ 



(II I -32) 


where iF-|(a s b,z) is the confluent hypergeometric function. 





III-10 


or 





(II 1-33 ) 


Hence, from (III-23), we have for the coefficients 



(II 1-34 ) 



ORIGINAL PAGE IS 

OF POOR QUALITY 


m-n 


Note that, since r(-k) = ±», k>0 (k an integer), when £>0, the first 
denominator will be either ±°° and hence, h^Q= 0. Also, when £= 0 and 
m.j+k=l, then -jF-j (0,b,z) = 1 , and again the first term drops out. Since m-j+k 
must be odd, we see that the first term drops out. 

We therefore have 




/0 m 

m-,/2 ~ P 1 
p i £ “q £!r(m-]+£+lJ 


^2£+k+m-j 
1 F 1 2 




x 


r 


( 2 - 21 - k - m-j} 
2 


(III -35 ) 


Now note that 


-|F-|(y- x, y, -z) = e" z -j F-j (x , y, z) 


(III -36 ) 


so that letting 


b 



we finally obtain 
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= 0 otherwise 


( I 11—37 ) 



In ( 1 11-29) 3 we now consider the case and all m-| and all k. 

the approximation 



Consider 


(III —39) 


See Figure I I 1-4 for a comparison of the approximation. From (III-29), we 
get 



(III —40 ) 






in-14 

Going through the same calculations as before, using equation (III-28), we obtain 
the result 






(III -41 ) 




and 



( I 11-42) 


we then have 
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P 1 


in-15 


£= 0 £!r£m-|+£+l 


fa 1 F 1 


k+m-j+2£+l 0 -6^ 


2 9 2 ’ . 2 

9 2 4 . 2 
2 <j„ + — 5— 

n 3 


A * 

, 2 "2 

o„ + 


'n 3 


-k+1 -m-j -2£ 


(III— 43) 


Hence, using 


h km-j 1 Jjj h km-|l " h km 1 l 


(II 1-44 ) 


8 = d 


along with the fact that 


lim 1 F-, {a, b, 6) = 1 
fr >0 1 1 

yields from (III-42) 



(II 1-45) 
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III-16 


Since 


Jr(n) J = ® n=0,-l ,- 2 ,, . . 


Hence, when £> 0 r(*)-*° and contributes nothing to hkm-j-j. Furthermore, 
when £= 0 and m-j+k= 0, the first term 


1 



= 0 


Also, when m-|+k=2, the gamma function is ® and again we get zero, 
we obtain for h. 


so that 


'hr 


A c 

2 a 2 
+ 


-k/2 


km-jl 


a n mi/2 | n 3 
d \T P 1 \ 2“ 


Pi 


/2fi 


1 F 1 


k+mi+2£+l g 


£=0 £ ! r C m -| +£+0 
-« 2 

’ 2 ,- 7 


4 2 4 


v. 


A 

2 A 2 
a n + T 


'-k+1 -m 1 -2£l 


(I 11-46) 
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EVALUATION OF AN INTEGRAL USED IN THE CLIPPER PROBLEM 


We start with a result due to Middleton [1]: 
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Using equation (IV-5), we obtain 
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Hence, we finally obtain 
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'1.0 INTRODUCTION 

In order to reduce the effects of stiction in the Ku-band 
deployed assembly and the resultant limit cycling, dither is being 
introduced. This dither is an autonomous sine wave, located at a fre- 
quency sufficiently removed from the mechanical assembly resonance to 
preclude excitation, which keeps the antenna in virtually continuous 
motion. An heuristic description of dither is given in [1]. 

Not only must a possible resonance excitation be considered, 
but the dither frequency cannot be too high or low. At very high fre- 
quencies, the mechanical assembly cannot repond and stiction is not 
affected. At low frequencies, the mechanical assembly excursion is 
excessive and introduces pointing errors. Hughes has picked a dither 
frequency of 17 Hz, which does not fall close to any known mechanical 
resonance. In the following paragraphs, we derive the mechanical motion 
due to dither and the resultant AM on the sum and difference channel. 

It is shown that worst-case AM on the sum channel, within the 3 dB beam- 
width, is about 0.4 dB; AM on the sum-plus-difference channel is about 
5 dB. 


2.0 DERIVATION OF ANTENNA RESPONSE 

The two axes of the deployed assembly antenna are excited by a 
17 Hz sine wave at sufficient current to produce the equivalent of two 
foot-pounds of torque. The response of the antenna can be derived from 
the angular equivalent of F=MA; L= Ie. The quantity L is torque, equiv- 
alent to force, I is the moment of inertia, equivalent to mass, and e is 
angular acceleration. Values of the parameters follow. 
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The angular excursion. 

8, is given by the double integral of 

e, as shown below: 


„ _ 2L cos tut 

1 2 cos 34 tt 

2 T 
<u I 

O X 

(34rr) I 


The peak-to-peak excursion is 


0 : — _ x _ 

(34irr I 

Thus, = 0.000876 radians, or 0 0502°, and e g = 0.000565 rad- 
ians, or 0.0324° The vector sum of these two gives a worst-case peak-to- 

peak angular excursion of 0.06°. The vector sum is not an unrealistic 
quantity since the two dither signals are presumably being derived from 
the same source, and the resultant motion is probably close to being 
in phase. 

3.0 SUM CHANNEL DEGRADATION 

In this section, we derive the sum channel degradation based 
on worst-case 0.060° peak-to-peak antenna motion due to dither. These 
calculations are based on the approximation to the antenna mam lobe pat- 
tern given below 

G(e) = 201og^^- , dB 

Taking a nominal 3 dB beamwidth of 1.6°, we can solve 
201og S1 q = 3 to find k = 1.74, with e in degrees and the sin(*) 
function in radians. Now the peak-to-peak gain variation can be computed 
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from 


AdB 


20 log 


sin 1.74 9 
1.74 6 


- 20 log 


sin 1.74 ( 9 - 0.06) 
1 .74(9 - 0.06) 


This is plotted in Figure 1. The points plotted, which we calculated in 
0.06° increments, appear to fall on a straight line. Since the computa- 
tions were done on a calculator with built-in linear regression, it was 
a simple matter to compute the best straight line fit to the points in 
the region of interest, 0.06° to 0.80°. In fact, the points are very well 
correlated to a straight line, with p = 0.999. The equation is 

AdB = 0.5940 - 0 030. 

If we want to maintain the total pointing loss less than 0.3 dB, 
the maximum allowable e from the (sin x)/x approximation is 0= 0.26° minus 
one-half the dither amplitude, or 0,23°. That is, the maximum steady 
state pointing error must be less than 0.23° to maintain the required 
0.3 dB maximum loss. 


4.0 DIFFERENCE CHANNEL DEGRADATION 

The AM on the difference channel can be estimated using the 
autotrack scale factor and the worst-case antenna motion of 0.060°. The 
scale factor is defined as 


Representative parameters taken from ADL test data are given below: 


^- 0.5 

A0 = 0.060, from prior calculation 
e 3 = 1.6, the 3 dB beamwidth 

G = 79.4, the antenna gain (38 dB) in volts/volt 

2i 




8 (degrees) 


Figure 1 Sum Channel AM due to Dither versus Angle Off Boresight 
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aG^ is the difference in A channel gain, G 2 - G-| » in volts/volt. Using 
these values, we find that G 2 - G-j = 1.49 volts/volt. In order to convert 
this to decibels, we need a reference value. From the test data, the 
difference channel gam at the point of steepest slope appears to be 
about 6 dB, or 2 volts/volt. Thus, G 2 = 2+ 1.49 = 3.49, or 10.86 dB. The 
AM near boresight on the difference channel is then 10.86 - 6= 4.86 dB. 

5.0 CONCLUSION 

The 17 Hz AM on both channels is insufficient to degrade track- 
ing or communications performance significantly since the angle-tracking 
filter is so narrow and the 17 Hz is well below the data frequency. How- 
ever, a minor improvement in steady state angle-track accuracy is required 
to maintain the 0.3 dB allowable pointing error with dither present. 
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FROM: R. Iwasaki FILE: 16067'^" Distribution 

SUBJECT: Potential Ku-Band Si delobe Reduction Recommendation 


During the recent Ku-band meeting at Hughes, sidelobe acquisitions 
by both the radar and communication systems were major problems which required 
many modifications to circumvent limitations in performance. The relatively 
high (-20 dB) sidelobes are an antenna characteristic that Axiomatix believes 
may be reduced by relatively simple means which may or may not affect the pro- 
gram schedule, but should at least be brought to the attention of Hughes man- 
agement as a possible alternative to consider for future STS missions. 

First, it is to be emphasized that no major modification is being 
considered. A simple test is proposed during a currently scheduled antenna 
pattern measurement that will quickly determine whether or not the change is 
even feasible. Only one additional pattern cut in the plane of the highest 
sidelobe is required, which should not impact the tight schedule. 

This modification is the second item recommended during Axiomatix 's 
antenna study. The primary fix— the RF fence to minimize mutual -coup! mg ef- 
fects in the monopulse feed — reduced the sidelobes from 13-14 dB to the pres- 
ent 20-dB level. At the time, this improvement was considered acceptable but 
recent experimental tests indicate that a wider acquisition dynamic range is 
desirable. Although antenna design is still somewhat of an art and therefore 
unpredictable, based on our antenna design experience which corrected the ear- 
lier unacceptable sidelobe problem, Axiomatix feels that an improvement to 
-22 dB is not unreasonable. This would increase the measured acquisition dy- 
namic range from 17 to 19 dB, thereby eliminating the acquisition threshold 
setting problem that will exist for missions beyond STS-7. 
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Axiomatix proposes a test in which a wedge-shaped ridge of heavy-duty 
aluminum foil is placed along both sides of the feed-support structure, as shown 
in Figure 1, to minimize the antenna aperture blockage effects which now appear 
to be the dominant causes of the high sidelobe levels. The rationale for this 
modification is to reduce the deleterious asymmetrical obstruction of the aper- 
ture by allowing some radiation to propagate past the feed, thereby allowing 
more destructive interference to occur with radiation impinging on the opposite 
side of the antenna. 

The wedge-shaped ridge can be pictured as having polarizing charac- 
teristics similar to those of the septum which allows electric fields that are 
oriented perpendicular to the ridge to be transmitted, whereas the parallel 
electric fields see an electrical short and are therefore reflected. This be- 
havior is valid even for structures larger than the wavelengths of the electro- 
magnetic radiation and therefore does not depend on any diffraction arguments. 
The linear polarization, consisting of both RHCP and LHCP components, will then 
contribute some offsetting aperture illumination versus none with the existing 
configuration. 

Incidentally, this technique is used in radome structural support 
designs to minimize blockage effects of metallic ribs. 

The exact taper of the wedge is somewhat arbitrary because the stow- 
age envelope of the antenna must accommodate the modification but, as a general 
rule, the longer the taper, the better the performance. If a metallized fiber- 
glass shell is used to shape the wedge, this structure would also serve as a 
thermal shield which will help stabilize the thermal environment surrounding 
the sensitive RF mixers, thus reducing the need for the electrical heater. 
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Frontal View of the 
Parabolic Antenna 


Side View of the 
Parabol ic Antenna 
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Cross Section of the 
Feed Support 


Modified Wedge-Shaped Cross Section 
of the Feed Support 


Figure 1 The Wedge-Shaped Blockage-Minimization Technique 
for Si delobe Reduction 
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S-BAND SYSTEM INVESTIGATIONS 

Task 2 of Exhibit B was oriented toward the S-band network communi- 
cation equipment which had already been through the major development stages 
before this contract effort started, but on-orbit tests with the TDRSS have 
taken place only in the last few months due to the delayed launch of the TDRS 
and the initial problems of the satellite achieving synchronous orbit. There- 
fore, analyses of the S-band network overall system performance were limited 
to acquisition times with the wide data filter and data off as well as with 
Shuttle G/T and EIRP calculations. These analyses are presented in this sec- 
tion. It should be pointed out that additional analyses related to problems 
which developed in ESTL and KSC tests are documented as part of Exhibit A, 

Task 15, which is presented in Section 4. 

Analysis of the Shuttle S-band acquisition times with the wide data 
filter and data off shows that the acquisition threshold is approximately C/N Q 
50.5 dB-Hz. 

A Shuttle S-band G/T of 29 dB and an EIRP of 15.6 dBW was calculated, 
as documented in Technical Memorandum No. M8206-3, "Shuttle G/T and EIRP Calcu- 
lations," dated June 21, 1983. This memorandum was generated in order to elim- 
inate any possible differences of opinion regarding G/T and EIRP. 
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July 9, 1980 



TO: Gaylord Huth 

cc: Q. Tu 

E. Dodd 
B. Batson 
B. Teasdale 

FROM: J. K. Holmes 

641-8600 

DATE: July 9, 1980 

SUBJECT: S-Band Shuttle Acquisition Times With The Wide Data Filter And 

Data Off 

SUMMARY : 

The S-Band Shuttle PN code despreader mean acquisition time 
is estimated for the case the data is turned off and the high data rate 
bandpass filters are employed in the despreader. 

Based on the available data, acquisition should occur at and 
above C/No = 50.5 dB-Hz. 

ANALYSIS 

The intent of this memo is to estimate the mean acquisition 
time of the Shuttle S-Band and PN code despreader with the high data rate 
filter switched in and the input data turned off. The only available 
acquisition time data, however, is for the high data rate mode with the 
high data rate signal being received Thus to estimate the desired 
acquisition time we must compute the filtering losses that occur in 
mhe pre-detection filters of figure 1 for the case of the high data 
rate signal in high data rate pre-detection bandpass filters. 

After the filtering losses are computed, it is then a 
simple matter to adjust the measured acquisition times from reference 1 
by the filtering loss to yield to the desired results. 

The S-Band Shuttle forward link PN code length is 2047 chips 
and has a rate of 11.232 mega chips/sec. The data is Manchester encoded 
so that if we assume a random data stream with the probability of a data 
one and a data zero being 1/2, independent of the previous bits, then the 
power spectral density is given by 



DAT~ 

rate 

CMP 


il 3 cSm 
AT 31 MHz 



DESPREAD CHANNEL 


Figure 1 5-Band Spread Spectrum Processor Functional Diagram 
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sin 4 (nfT/2) 

(nfT/2)2 
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where P is the Manchester data stream power, and T is the Manchester 
symbol duration. 

For zero carrier and code doppler we have that the filtering 
loss in the bandpass pre-detection filter (assumed to be ideal) is given by 


L = 



sin 4 (nfT/2) 

■ '~m/rr 


df 


where BW is the RF bandwidth which is 480 KHZ for the high data rate 
mode. Numerically evaluating the integral yields L=0 72=-1.43 dB. We 
conclude that with zero doppler and no data the despreader sensitivity 
increases by 1.43 dB. Because of the inherent bias in the code loop 
during acquisition the S-Band Shuttle despreader biased the code loop 
frequency to -200 chips per second so that a true code doppler of 
+400,0, and -400 chips per second respectively appears to the acquisition 
circuitry and the code loop as +600, +200, and -200 chips per second, 
respectively. 

Now at a code doppler of 400 cps, corresponding to 
about 60 KHZ carrier doppler, the filtering loss is given by 


+60 KHZ 


-/ 


PT 


sin (nfT/2) 
(JlfT/2) 


df 


— |^ + 60 KHZ 


which when evaluated numerically is given by L = -1.52 dB 



By replotting figure 9 of reference 1 (which is the average 
code acquisition time versus C/No (dB-Hz).} with the filtering loss 
subtracted from the abscissa for each acquisition time we obtain figure 2 
which is our desired result. Notice- from figure 2 that threshold appears 
to be around 50.5 dB-Hz. 
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TO: G. Huth 

FROM: P. Nil sen 

. SUBJECT: Shuttle G/T and EIRP Calculations 


DATE: June 21, 1982 

COPIES: J. Johnson 
M. Zeal on 


On June 17, 1982, I met with Mike Zealon, Glen Yates and Dick Davis 
to review the RI methodology for Shuttle S-band G/T and EIRP calculations. 

Also, we reviewed the RI calculation of RF feedline insertion loss. The RI and 
Axiomatix methodologies for these calculations were essentially the same, with 
the minor exception that RI did not include contributions to effective noise 
temperature beyond the preamp. The formula used by Axiomatix for effective 
noise temperature, namely, 

O-o - l}T n L 1 

T eff = T A + 0-! - O T 0 + C NF PA " ^ T 0 L 1 + ^ pA ~ 

A ^ NF XPDR " 4 T o L i L 2 

g pa 


where 



G PA 

NF XPDR = 


antenna temperature (75°K) 

loss between antenna and preamp (4.34 dB) 

loss between preamp and transponder (0.5 dB) 

preamp noise figure (2.6 dB) 

preamp gain (20 dB) 

transponder noise figure (7 dB) 


did include the temperature contribution for the cable between the preamp and 
the transponder as well as the transponder 's contribution. 

However, we see that, after making these calculations with and with- 
out these last two terms, the difference in G/T is only approximately 0.15 dB, 
with the more precise result being G/T = 31 - 2 = 29 dB for 80% coverage. 
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The EIRP calculation is given by 


EIRP - P t *G ant -L t 


where Py is the transmitter output power (20 dBW), G^y is the antenna gain 
(80%) (2 dB), and Ly is the transmit circuit loss (6.4 dB). 

The calculated EIRP is 15.6 dBW for 80% coverage. It should be 
noted that, in both the EIRP and G/T calculations, the transmission line loss 
due to VSWR was included in circuit loss values. The VSWR loss was found 
according to the attached nomograph. 
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VSWR Nomograph #2 
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ORB ITER ANTENNA STUDIES AND INVESTIGATIONS 


Task 8, Orbiter Antenna Patterns, was composed of two separate areas 
of investigation. First, were analyses to improve the performance of the S- 
band quad antennas. Axiomatix was heavily involved in this effort and provided 
many significant contributions to the improved coverage performance. The sec- 
ond area investigated under this task was to develop techniques by which to de- 
termine the actual Orbiter antenna gain and coverage from on-orbit measurements. 
After an initial investigation, the results were presented in Axiomatix Report 
No. R8009-1, "The Concept of Maximum-Entropy Filtering and Its Feasibility as a 
Candidate for Shuttle Antenna-Gam Processing," dated September 5, 1980. Fol- 
lowing the review of this report by NASA, Axiomatix was directed on June 1, 1981 
to terminate further analysis m this area as there seemed to be little use for 
on-orbit measurements derived from several ground stations which were calibrated 
differently. 

In order to investigate methods by which the S-and quad antennas 
would meet the relaxed gam-coverage requirement of 4 dBci over half the sphere, 
Axiomatix extended many of the concepts used in the switched two-beam quad an- 
tenna, taking advantage of the Rockwell -developed basic hardware. The outcome 
of this study has been the evolution of a new generation of variable-gain, 
mul tiple-switched-beam antennas, the principles of which can be applied to the 
improved Shuttle Orbiter S-band quad antennas. 

The switched two-beam quad antenna developed by Rockwell was an im- 
provement over the earlier single-beam quad antenna built by Watkins- Johnson 
but, due to over-optimistic computer analysis projections of the antenna pat- 
terns, the actual performance results were disappointing. The primary limita- 
tion of the gain peformance was the selection of only two elements and two beam 
positions, the minimum required for switched-beam operation. Using the same 
type of hardware configured in a slightly different manner, Axiomatix devel- 
oped a two-dimensional cross array that could be electromechanical ly scanned 
in both the roll and pitch planes as 2^x2^ discrete-beam positions, where M 
and N are the number of doubl e-pol e/double- throw electromechamcally latching 
switches in each plane. The advantages of using this type of switch rather 
than electronic-phase shifters are: (1) high power-hand! ing capabilities, (2) 

low insertion losses, (3) absence of electrical "holding" power, and (4) the 
binary nature of the double-throw switch which can readily be converted to a 
binary code for a specific beam position. 



As part of Axiomatix participation in the development of S-band quad 
antennas, Axiomatix personnel attended the design reviews and technical exchange 
meetings. In Axiomatix Report No. R8007-2, "S-Band Quad Antenna PDR Evaluation," 
dated duly 23, 1980, recommendations are made to improve the gain-coverage per- 
formance since the two-beam approach did not achieve the specification. This 
report attempted to illuminate the dilemma by suggesting, first, some concave- 
lens design concepts for the existing randome which would broaden the dipole 
beams in the array plane, an essential feature to correct the present two-beam 
coverage deficiencies. If the proposed three-beam scheme was to be used instead 
of the two-beam design, Axiomatix recommended that a switched four-beam config- 
uration be seriously considered as an alternative system since the hardware and 
software would be extensively modified for the three-beam scheme. Finally, a 
novel RF fence design usign the total -reflection concept of critical angle was 
outlined, which should substantially reduce the amount of mutual coupling and 
deleterious interaction effects between the two dipole-array cavities for all 
switched multiple-beam configurations. 

Further investigations of the four-beam configuration to simplify 
the design are presented in Axiomatix Report No. R8106-4, "Simplified Orbiter 
S-Band Quad Antenna Planar Array Candidate Configurations," dated June 25, 1981. 
Also, improvements to the lens concept required to broaden the beam of the 
two-beam design are presented in Axiomatix Report No. R8112-1, "A Protruding 
Wedge-Shaped Dipole Element Modification For The Orbiter S-Band Quad Antennas," 
dated December 9, 1981. Finally, the concept of the two-dimensional cross ar- 
ray that would be electromechamcally scanned m both the roll and pitch planes 
at 2 M x2 N discrete-beam positions is developed in Axiomatix Report No. R8112-4, 

"A Large Switched-Beam Planar Array," dated December 15, 1981, and Axiomatix 
Report No. R8209-3, "A Variable-Gain Multiple-Switched-Beam-Antenna Configura- 
tion," dated September 22, 1982. To implement a variable-gain multiple- 
switched-beam array, the beam-switching logic must be simple. This switching 
logic is presented m Axiomatix Technical Memorandum No. M8209-4, "Binary 
Code Designation of Multiple-Switched-Beam Positions," dated September 30, 

1982, as well as in Axiomatix Technical Memorandum No. M8209-5, "Simplifi- 
cation of Antenna Beam-Switching Logic," dated September 30, 1982. 
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1.0 


INTRODUCTION 


Maximum entropy filtering has been suggested by John Tranter as 
a candidate approach for filtering the raw antenna gain measurements made 
during early Shuttle flights. The motivation for this suggestions is the 
fact that this method has been shown to yield power spectral density (PSD) 
estimates which have higher resolution than traditional periodogram (Wien- 
er) or weighted periodogram estimators. This advantage is especially true 
when the measured data is sparse. 

The following discussion offers an explanation of maximum entropy 
filtering and presents several examples. Most important, an explanation 
of why this technique is not suited to the determination of the Shuttle 
antenna gain pattern is presented 

2.0 DISCUSSION 

The expression "maximum entropy" is a misnomer m the sense 
that one generally maximizes the information, not entropy. However, this 
expression is used in connection with the channel capacity of a communi- 
cation system where maximum entropy corresponds to maximum information 
transmission rate. In statistics and probability theory, the entropy 
H(X) of a random variable X with probability density function (PDF) f(x) 
is defined by 


H(X) 



f (x) £nf(x) dx 


0 ) 


The entropy H(X) is closely related to the disorder or uncertainty associ- 
ated with making a realization X. For that reason, maximizing the entropy 
is a method for finding the distributions that represent high uncertainty 
or, equivalently, a state of high ignorance Thus, given a set of con- 
straints, the maximum entropy method provides a candidate distribution 
which confirms the constraints without assuming anything more. 

For example, let the constraints be 



f(x) dx = E(X) 


m r 


o 


( 2 ) 
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Then the distribution that maximizes (1) subject to (2) is 


f(x) 


I_ ~ x/m 0 

m o 


(3) 


Notice that, since constraint (2) specifies mg, the maximum 
entropy distribution contains no additional information and is therefore 
maximally random. 

As another example, consider maximizing (1) subject to 



x f(x) ds 


E(X) - m Q 



Cx - E(x )} 2 


f(x) dx 




(4) 


In this case, the distribution that maximizes (1) subject to (4) is 



Notice that, since and a Q are specified by constraint (4), the maximum 
entropy distribution given by (5) for this case contains no further infor- 
mation on the underlying PDF of X and is thus maximally random. 

Similarly, extending our concepts to a random process, it can be 
shown that the maximum entropy random process with constraints on only the 
first and second moments is a Gaussian process. Also, the entropy for a 
stationary Gaussian process may be shown to be given by 


f r, 

% / 


log f ( A) dx 


H 


( 6 ) 
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where 2) is the domain in the frequency space and X^ is the spectral 
measure of @ . 

Let X(t ) be the stationary Gaussian process in question, with 
zero mean and with covariance matrix of observations X^t-j}, X(t,p,..., 
X(t n 3 sampled at times t-j, tg,...,^ given by 




where, for brevity, we write X(t^) = X^ , i=l,2,...,n and 


r 


k 



e j2irXK f dx 

Discrete 
Time PSD 


( 7 ) 


The maximal entropy spectral density estimate subject to the 
constraints 


r K = r K° * K = 1 ,2, ’ ,n ( 8 ) 

is obtained by maximizing (6) subject to (8) and will turn out to be a 
rational function of jr^|. 

Thus, in general, the maximum entropy principle works if the 
constraints are practically noiseless or error-free and, if they must be 
estimated, they should be based on a large data base of observations 
corresponding to repeated cycles or replications The specification of 
the constraints could also be based on proven prior knowledge arising 
from past and cross-sectional experience. 

We therefore resort to smoothing procedures for the antenna 
gain data. The unacceptability of the maximum entropy method of esti- 
mating the gam at missing points will be demonstrated for the one- 
dimensional case. 
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(9) 

f 0g > 

0 

(10) 

N 



2 

i a i 

= 1 , N > n 

(11) 


We will now estimate f (x } at the missing points by using the maximum 

«J 

entropy principle, i.e., choose f(x ) such that 

J 

N 

H(f) = - J fCx^ log fCx 7 0 02) 

i=l 

is maximized subject to (9), (10) and (11) 

From (9), (10) and (11), it follows that it suffices to maximize 

N 

H(f) = - J fCx,D log fCx.D (13) 

i=n+l 

But maximization of (13) subject to (10) and (11) or, equivalently, 

(11 1 ) where 

N 

2 =1-0 (IT) 

i=n+l 

where c = c, + c 0 +.. + c , is attained by 
1 l n 

fCsD = » i=n+l ,.. . ,N (14) 


i.e., uniform distribution for' the missing points. 
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This solution, while mathematically correct, is totally 
unacceptable for the specific problem at hand since it does not take into 
account the smoothness properties of the gain function. Furthermore , it 
is easy to demonstrate by examples that, in most cases, the missing gain 
distribution between two measured points is not a uniform distribution. 

In the next communication, we specifically address the problem 
of smoothing the gain data' in the presence of noise, and establish the 
rationale, obtain acceptable procedures and develop algorithms for 
implementation. 

3.0 CONCLUSIONS 

The maximum entropy filtering technique is not suitable for 
processing the Shuttle antenna gam measurements. Therefore, it is rec- 
ommended that no further investigation of this technique be undertaken. 

Axiomatix recommends that a more "ad hoc" approach to process- 
ing the raw gam measurements be utilized An example of such a technique 
is the nearest neighbor interpolation to fill in missing data combined 
with a "pseudo-Kalman" filter for smoothing data from multiple passes and 
multiple stations. This is the approach which Axiomatix is now pursuing. 
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1.0 INTRODUCTION 

At the S-Band Antenna Preliminary Design Review (PDR) held at 
Rockwell /Missile Space Division in Anaheim on July 15-16, 1980, the 
main issue was the proposed change from the initial switched two-beam 
configuration to a switched three-beam configuration for the S-band 
quad antennas. The other major issue of note was the poor antenna 
patterns of the lower S-band quad antennas. Further work on these 
antennas cannot proceed until a firm decision is made whether or not to 
pursue this new approach because of the required alterations in the hard- 
ware and software of the system. The meeting concluded with a plan to 
resolve this issue within ten days after the pertinent system personnel 
are queried as to the ramifications of the changes. 

The other aspects of the S-Band Antenna Program appear to be 
progressing satisfactorily except for some specification deviations which 
must be resolved by either waiver or additional engineering adjustments, 
such as the axial ratios of the quad antennas. However, in light of the 
major impact of a new three-beam configuration and the obvious inadequacy 
of the low quad antennas, these latter problems warrant immediate attention. 

2.0 TECHNICAL SUMMARY OF THE PDR 

From the data supplied during this PDR, it appears that the three- 
- beam approach was analyzed within two weeks of the last technical exchange 
meeting held May 14, 1980. The implication is that the two-beam approach 
was deemed unsatisfactory very soon after the preliminary antenna pattern 
measurements were made and, since completed three-beam patterns were 
presented, that the decision has been tentatively made and the two- beam 
approach has been essentially abandoned. 

The two-beam approach resulted in a 65% 4.0 dB spherical coverage, 
for the upper quad antenna (versus the coverage specification of 85%) and 
a modified two-beam approach which moved the aft beam forward to cover the 
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actual roll plane (by using a -66° instead of a -90° phase shift) gave 
a 64% 4.0 dB coverage. The lower quad gave a 35.1% 4.0 dB coverage, 
indicating serious measurement and/or test antenna design problems. 

The corrective measures necessary to improve the two-beam 4.0 dB 
coverage include broadening the beam to redistribute the peak 6.5 to 7.0 dB 
gam more evenly over the specified coverage area. The technique attempted 
by Rockwell MSD was to place a ‘'parasitic' 1 ring on the radome to enhance 
generation of a third-order mode to put more radiative energy away from the 
axis. This parasitic ring was intended to function similar to the directional 
reradiators of the Yagi antenna, but the alternate phasing of the third-order 
mode greatly complicates the design. The modification did not produce the 
desired results and it appears that at this point an alternate concept 
(e.g. three beam) was seriously evaluated. 

3.0 POSSIBLE ALTERNATIVE DESIGN MODIFICATIONS 

Some possible approaches to broaden the two beams might be considered 
in the future if the three-beam scheme cannot be realized. One is using a 
concave-lens configuration on the radome. It was mentioned that Rockwell 
MSD has experimented with this concept earlier, but the concave surface was 
on the opposite side away from the dipole, which is not the ideal situation 
if it is assumed that the dipole radiates quasi-spherical waves and the 
radome surface is to be in intimate contact with the Thermal Protection 
System (TPS) tiles. Further, the argument of high cross-polarization levels 
for broader beams must be evaluated more closely if the beam is selectively 
broadened, primarily in the array plane where the beam is narrowed by the 
two-element array. 

Another reason for using an inner radome aperture concave 
surface is the reduction of the cross polarization. The high dielectric 
constant U r =4) radome behaves as a polarizer in the present dipole config- 
uration»and therefore unequal magnitudes of orthogonal linear polarization 
(equivalently orthogonal circular polarization) are transmitted through the 
radome. The reflected components of polarization create the higher VSWR 
mismatch measurements, which also appear to exceed the desired performance 
specifications. Impedance matching to decrease the VSWR by design is dif- 
ficult to achieve within the limited area of the aperture, but the concave 
inner radome surface will most probably decrease both the VSWR and the 
cross-polarization levels. 
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The reason for the selective polarization reflection is the 
optical principle of reflection at a dielectric interface or, in elec- 
tromagnetic theory, impedance matching. 



Fraction of incident power reflected from the surface of an infinite 
dielectric slab vs. angle of incidence: (a) electric vector perpendicular 
to the plane of incidence, (b) electric vector in the plane of incidence: 
n = 1.6. 


Because of the quasi-spherical phase fronts radiating from the dipole 
element, the incidence angle varies from the central axis and selective 
polarization transmission occurs. The same phenomenon also occurs at 
the outer radome surface where it is in contact with the TPS tiles and 
should be considered in determining cross polarization effects. 

The increased path lengths through the outer edges of the radome 
-also distorts the quasi-spherical phase wave front which is the basis for 
array theory. In summary, the radome should be analyzed as a polarizing 
lens system. 

By broadening the individual dipole pattern in the array direction, 
power redistribution can be achieved since the excess power m the peak 
gam (6.5 to 7.0 dB) is actually being selectively rearranged to reduce 
the pattern-multiplication factor in the array plane. Therefore, the beam 
broadening essentially tries to re-establish the symmetrical single-dipole 
pattern, but uses the two-element phased array to electronically tilt the 
resultant pattern in the desired directions by switching the appropriate 
phase shifts. 
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In order to reduce the gain and broaden the antenna coverage, 
a concave aperture interface is proposed. Snell's law, which bends 
the ray closer to the interface normal in a denser medium, is now mani- 
pulated to spread the energy outwards since more of the incident spherical 
wavefront energy is transmitted into the dielectric at divergent angles. 



PRESENT PLANAR RADOME 


CONCAVE SURFACE 



PROPOSED PLANO-CONCAVE RADOME 



LOWER QUAD CONVEX OUTER SURFACE RADOME 

The dielectric radome, presently polyimide E-glass, is machined 
such that the dielectric interface to the dipole antenna cavity is flush 
with the cavity aperture. The shape of this radome, since the aperture 
is indented to accommodate the radome ,, resembles a convex lens which focusses 
incident radiation. This effect is even more apparent when the dielectric 
constant of the dielectric material is higher than anticipated, as the case 
when the polyimide E-glass replaced the preliminary fiberglass material, 
especially in the case of the lower quad antenna where the outer radome 
surface is slightly curved. The net result of this configuration is a 
plano-convex lens which focusses radiation, or equivalently, increases the 
gain. This focussing effect of the convex radome surface is further enhanced 
by the quasi-spherical waves radiating from the dipole element which is located 
m close proximity to the radome. In order to compensate for this curvature, 

(if this convex lens shape is indeed found to contribute to the poor lower quad 
pattern since it is the only difference from the upper quad antenna), a concave 
inner radome surface will reduce the amount of focussing. Geometric ray 
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tracing can be useful for determining the overall lens effect. This 
technique can be extended to produce divergent beam broadening in the 
lower quad with increased curvature of the inner concave lens surface. 

This beam broadening technique need not be symmetrical. In 
fact, the concave curvature can be tailored to the broadening of the 
beam in desirable directions. The advantage of this type of modification 
is that the existing quad antennas can be used and only the dielectric 
radome, which is readily replaced, is modified. 

The lens effect of the radome slightly affects the effective 
location of the dipole element phase centers. The convex outer radome 
surface of the lower quad decreases the dipole spacing since the dipoles 
are not located on the longitudinal axis of the lens, which might contri- 
bute to the poor results. For the upper quad with the planar outer radome 
surface, this effect is less dramatic since, although the "depth" of the 
phase centers into the Orbiter varies, the dipole spacing remains constant, 
and the earlier analytical models are still valid. 

4.0 RF FENCE DESIGN 

Although mutual coupling between dipole elements is not obviously 
a major problem, it is axiomatic that it exists for adjacent radiators with 
a common radome. During the recent performance tests, the mutual coupling 
is indicated by the frequency-swept input return loss measurements where 
-10 dB was not uncommon. A desirable design goal then is to attempt to 
minimize this effect as much as possible, especially since the measured 
performance is unsatisfactory This problem is complicated in the quad 
antenna design because of the limited space available A technique to 
surmount these difficulties with minimal effort is suggested here which 
can be utilized for the switched two- or three-beam case (and even the pro- 
posed switched four beam case to be discussed later). 

The optical basis for this RF fence is the physical principle of 
the threshold-critical angle for total reflection which exists for electro- 
magnetic radiation propagating from a dense (high dielectric constant) 
medium to a less dense one, as is the case of the radome/free- space interface 
of the S-band antennas The threshold-critical angle e c is derived from 
Snell's low of refraction to be 
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where the angle e c is measured from the boundary normal and e r is the 
relative dielectric constant of the dense medium; in this case., polyinride 
E glass with £^4. Thus any radiation incident on this boundary at a 
greater glancing angle ( e >_e_ ) is totally reflected, constituting an RE fence 
in the classic definition. 

The near-field radiation pattern is difficult to characterize. If 
the radiated waves can be approximated by quasi -spherical waves and geo- 
metric ray tracing is employed, this type of concept of RF fence to isolate 
the two dipole-radiating elements becomes very attractive since there are 
few readily implemented isolation mechanisms that are frequency independent. 

The design criteria for this particular configuration is to have all 
incident rays greater than 30° from the boundary normal so that total re- 
flection occurs. The separation between the two dipole cavities is small, 
but the use of the smaller diameter iris or circular aperture provides suf- 
ficient room to mill the triangular groove into the radome, restricted pri- 
marily by the thickness of the radome and the requirement to retain structural 
integrity since the groove is prone to fracture. To minimize crack initiation, 
the production model should have a minimum milling radius at the tip to reduce 
stress concentrations. The depth of this milling groove is the primary factor 
in determining the degree of isolation, but even a shallow groove will disrupt 
surface waves at the radome/ground plane boundary that result in mutual coup- 
ling. 

5.0 SWITCHED FOUR-BEAM CONFIGURATION 

The switched three-beam configuration proposed by Rockwell/MSD 
attempts to alleviate the apparent design problems associated with obtain- 
ing adequate 4.0 dB coverage. By adding the third beam to cover the gap or 
"cusp" between the initial switched two-beam configuration, the coverage 
problem is essentially solved, at the expense of added complexity. At the 
PDR the consensus appeared to favor this approach as the one with the high- 
est probability of successfully meeting the coverage specification. The 
impact of such an abrupt design change, however, has not yet been fully 
assessed, and a firm decision to change to the switched three- beam system 
was delayed until the responsible system and software engineers could be 
queried. 
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The implications of such a design change require further examina- 
tion, at least in principle, to ensure that engineering requirements are 
optimally satisfied, especially for the added complexity. Besides the 
appropriate phase shifts, an additional identical RF switch is inserted 
in series with the switched two-beam configuration to provide the zero 
phase-shift differential for the broadside array generating the center 
beam. However, one RF switch position is redundant and unused The 
question arises whether or not this switch position can be employed in 
a useful manner, which brings up the point of the feasibility of a switched 
four- beam configuration with the identical hardware and similar software 
alterations. 

The proposed switched four-beam would have the primary advantage 
of a larger area coverage with higher gam, taking full advantage of the 
peak gam which is approximately 6.5 to 7.0 dB (even higher at broadside). 

The differential gain over the 4.0 dB coverage specification of 2.5 to 3.0 dB 
over the additional fourth switched beam area may be critical for some margi- 
nal communication links, and, if it was available for little extra software 
development (which has to be revised anyway), the consideration of the switched 
four-beam approach might be justifiable. 

A logical switched- beam configuration might use an arbitrary center 
line slightly offset from the roll plane (90°) so that one of the two central 
beams provides adequate coverage in this critical region. The other central 
beam would then be electronically and mechanically tilted slightly forward 
Since these inner two beams are closer to broadside, they would have higher 
gains and narrower beamwidths in the array plane than the outer beams (yet 
experience some beam broadening for broader coverage than the third 
. broadside center beam of the switched three-beam configuration). The outer 
two beamwidths can therefore be positioned for the desired coverage, with 
a bias for more coverage towards the nose of the Orbiter. The additional 
flexibility provided by the fourth beam simplifies the coverage "filling in" 
process, an important consideration for the lower quad pattern results where 
obvious problems exist. 

Two switch settings of one RF switch then might be "fore" and "aft" 
of the imaginary center line and then the other two switch settings of the 
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other RF switch would be "inner" and "outer", denoting the relative 
locations of the beam. The appropriate phase shifts for the various 
beam settings can be readily determined. The availability of the four 
distinct switched beams from the two serial RF switches is not conceptu- 
ally unreasonable and makes maximum use of the proposed hardware modifi- 
cation. 

6.0 CONCLUSION 

There are two possible solutions to the switched-beam S-band quad 
antenna problem. One is to continue to pursue the original switched two- 
beam quad- antenna concept which was initially selected for its straight- 
forward simplicity. The problem with the two-beam approach is that the per- 
formance specifications have not been achieved since only 65% of the 85% 

4.0 dB coverage requirement has been measured. Without new design concepts 
to broaden the dipole patterns and therefore the coverage, the two-beam 
effort was abandoned in favor of the more promising three beam scheme which 
which provides the necessary flexibility at the cost of added system complex- 
ity, which is presently being evaluated. Obviously the three beam scheme 

is more likely to be successful and has nearly attained the desired 85% 
coverage for the upper quad, according to the 78.5% coverage of the initial 
three-beam measurements. 

This paper has attempted to illuminate the dilemma by suggesting 
first some concave lens design concepts for the existing radome which would 
broaden the dipole beams in the array plane, an essential feature to correct 
the present two beam coverage deficiencies. This concave surface also reduces 
the reflection and polarizing properties from a dielectric boundary, which 
minimizes polarization effects and therefore decreases the axiaili ratio. In 
.the likely event that the more promising proposed switched three-beam scheme 
is selected instead, it is recommended that a switched four- beam configura- 
tion be seriously considered as an alternative system since the hardware 
and software for the three-beam scheme conversion have to be extensively 
modified anyway. The addition of another beam position provides higher gain 
levels m that particular coverage area (2 to 3 dB of communications link 
margin) and therefore obviously increases the 4.0-dB coverage area. And 
finally, a novel RF fence design using the total reflection concept of criti- 
cal angle has been outlined which should substantially reduce the amount of 
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mutual coupling and deleterious interaction effects between the two 
dipole-array cavities for all switched multiple-beam configurations. 
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A recent presentation of more breadboard evaluation configurations 
for the Orbiter S-band quad antennas considered the mechanically switched 
three-element/four-beam system and two electronically phased eight-el ement/15- 
beam and 16-el ement/20-beam systems. -The primary reason for proposing these - 
systems is that the current two-element/ two-beam system gives only 51%, 4-dBci 
coverage instead of the 85% coverage specified. 

As stated earlier, Axiomatix feels that the present two-el ement/two- 
beam configuration has not been optimized since 2 dB excess mainbeam gain indi- 
cates that the radiated energy is not uniformly distributed to the 4-dBci level. 
However, since Rockwell deems any further suggested modifications to the present 
inadequate system as infeasible, by definition, a configuration change is there- 
fore mandatory if further performance improvement is desired. 

Using the available experimental pattern measurements of the existing 
dipoles, the Axiomatix-proposed three-element/four-beam mechanically switched 
configuration satisfied the 85%, 4-dBci coverage. However, it appears that 
even higher gam coverage is desired since the phased-array approach has been 
suggested. "Are the phased-array configurations proposed truly satisfactory 
from a practical viewpoint 7 " then becomes a reasonable question. 

Since hardware and software complexity are directly related to the 
number of elements in the array, we first consider the simplest two-element 
case. If one element is directly fed and the other has the phase-shifting ca- 
pability, the array will scan. The two-beam position is implemented by a me- 
chanical switch interchanging a discrete phase shift between the two elements. 
This is the present S-band quad configuration. 

If more than two positions were achievable, higher percentage cover- 
ages would be attainable. A multiple-value phase shifter capable of high-power 
operation was therefore required. But, if more positions are available, it be- 
comes more desirable to have higher gam since little performance improvement 
occurs with a wide beam, which has a high degree of beam overlap between posi- 
tions. For this reason, the three-element array was proposed and, with two 
mechanical switches m series with the appropriate phase shifts, analysis indi- 
cates that 85%, 4-dBci coverage is achievable. Since the center element is 
directly fed, some compensation for the additional switch loss occurs. This 
configuration satisfied the quad-antenna-coverage specifications, but was un- 
acceptable due to the higher cost and delayed schedules involved. However, if 
higher gains are desired, it is still not totally impractical to consider using 
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another mechanical switch in series to obtain eight beam positions with this 
configuration. Although the 4-dBci percentage coverage will not be substan- 
tially improved, the higher (> 4 dBci) gain coverages will be greatly increased. 
However, the planar-array approach allows more design freedom in achieving the 
high-percentage-gain coverage. 

The emphasis on higher than 4-dBci gain coverage appears to be the 
motivating factor for suggesting the use of planar arrays. Extending the pre- 
vious arguments, we find that higher gains from more elements and more beam po- 
sitions, in both pitch and roll planes, will result in higher percentage, higher 
gam coverage. There is no question that this is a valid assumption; the ques- 
tions now really are: "has the gain-percentage coverage performance requirements 
changed?" and "what is the feasible limitation on the array complexity?” 

Application of a two-dimensional planar array creates greater flexi- 
bility in gain coverage than simply adding beam positions in a linear array. 

The beamwidth in the orthogonal dimension is decreased, with a corresponding 
increase in gam. The simplest configuration for a planar array might be a 
triangular arrangement, as shown in Figure 1, where the two lower elements are 
identical to those presently being built. 



Electronic 
Phase Shifter 


Figure 1. Triangular Planar Array 
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If two-beam switching can now be incorporated with respect to the 
third element, is is then possible to generate four beams, with the amount of 
beam overlap determined by the degree of phase shifting. Note that this par- 
ticular three-element configuration is most readily implemented by using an 
electronic phase shifter in the third element, with mechanical switching still 
feasible between the other two elements. 

With the mechanical -switching arrangement used here for explanatory 
purposes, it is more feasible to consider using a four-beam/four-element array 
set in a square configuration, as depicted in Figure 2. 

Note that a tradeoff exists between the number of beam positions and 
the number of elements since there is more beam overlap for more beam positions 
using the wider beamwidths from the smaller number of elements. In order to 
increase the gain coverage, another possible configuration would be another 
three-element array superimposed orthogonally on the first three-element array, 
as shown in Figure 3. The center element is still directly fed, with the four 
outer elements determining the beam positions by appropriate phase shifting. 

This configuration is the largest one considered here since any additional ele- 
ments provide illumination taper and higher gain, but at the expense of consid- 
erably higher costs and complexity. 

Using the earlier arguments for the linear array, it is now possible 
to consider a four-beam planar array using the earlier switched two-beam concept 
in both orthogonal directions with two mechanical switches, as illustrated in 
Figure 3(a). The gain coverage is substantially higher than before because of 
the narrower beamwidth in the orthogonal direction. Similarly, since the three- 
element linear array was suitable for four beam positions, it is feasible to 
consider 16 beam positions using a total of four mechanical switches, as shown 
in Figure 3(b). Since the center element is directly fed and each orthogonal 
linear array is fed in parallel, the switch losses are not as significant. 

The phase-shifting mechanism has not really been specified except 
that, since mechanical switches have been utilized in the present configura- 
tion, it has been retained in principle. Mechanical switching is preferable 
to electronic switching because it has low-loss characteristics and possesses 
the required high-power-handl mg capabilities. In addition, being the bistate 
"latching type," no holding current (as m the case of PIN diode phase shifters) 
is required. Since rapid scanning is not a critical factor, the mechanical 
switch is suitable for these applications. The only situation wherein the 
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Figure 2. Square Planar Array 



Figure 3. Five-Element Planar Array 
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electronic-phase shifting might be seriously considered in the planar arrays 
discussed here is for the triangular three-element array. 

In conclusion, if higher gain levels at 85% coverage are desired for 
the S-band quad-antenna system, a five-el ement/16-beam mechanically switched 
system is proposed. The microwave circuit losses would be about the same (1.8 
dB) as the three-element/four-beam linear-array system proposed earlier since 
the switching systems are in parallel and therefore identical. The only addi- 
tional new component is a power splitter used to distribute the power in the 
orthogonal direction, which is easily incorporated into the input strip! me. 
Because the number of beam positions is nearly identical to the eight-element 
(4x2) phased array proposed by Rockwell, it is anticipated that the 85% gam 
coverage would also be about 6.5 dB. However, since the microwave circuit 
losses are significantly lower (1.8 dB versus 3.6 dB) because mechanical 
switches are used, the overall 85% gain coverage should be accordingly higher. 
Also, since there are no electronic-phase shifters, no electrical "holding" 
power (24W - 48W m the el ectromc-phase configurations proposed) is a definite 
advantage in space operations. This configuration also has the advantage that, 
since only three linear elements form the beam, the beamwidth is quite large so 
that extreme pointing accuracy is not required, and there is a great deal of 
pointing margin before a communication link is lost. 

Although it is considered somewhat impractical, by adding another 
switch in series with both feed lines (which adds only 0.2 dB more circuit 
loss), it is conceivable to have a 64-beam/five- element array which will def- 
initely have the highest gain-percentage coverage— if that is indeed the pri- 
mary operational requirement. 
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1.0 INTRODUCTION 

A meeting was held at the Rockwell facility in Anaheim on 
November 24, 1981 to discuss some alternative R&D modifications for the 
Shuttle Orbiter S-band quad antennas. The present 4-dBci gain percentage 
coverage of approximately 50% is substantially below the 85% specifica- 
tion, with the coverage at one frequency, 2.217 GHz, measured to be 38%. 
This frequency sensitivity is somewhat disconcerting since the production 
flight models are fabricated identically, but there is optimism that the 
coverage degradation might be "tuned" out by mechanical adjustments. The 
cause of the frequency sensitivity is not known and therefore cannot be 
easily remedied without extensive testing. The low-percentage coverage 
problem, however, is more amenable to improvement and, after extensive 
discussions, a preliminary modification program, funded internally by 
Rockwell, was formulated, with the antenna measurements to be conducted 
at JSC. This interim report attempts to document the results of this 
meeting such that the approach and justifications are clearly delineated. 

2.0 PROBLEM 

The prototype S-band quad antennas have a 4-dBci percentage 
coverage of only 50% instead of the desired 85% coverage. The selection 
of the two-element/two switched-beam quad configuration early in the pro- 
gram arose from some unrealistic projections on reducing circuit losses 
as well as misunderstanding the interpretation of the data, specifically 
the gain levels. Since the configuration cannot be readily altered, the 
obvious question is: can any suitable modifications be suggested which 

might improve the percentage coverage? The extent of improvement is dif- 
ficult to assess, of course, but some ideas were aired that appeared to 
have a valid basis for improving the percentage coverage, off-axis axial 
ratios and radiation efficiency. 

3.0 RECOMMENDATIONS 

The flush-mounted dipole elements appear to have an intrinsic 
limitation on the off-axis axial ratios because of the small electric 
field component radiated normal to the plane of the aperture. In order 
to circumvent this problem, a protruding-antenna concept was introduced 



2 


which extends the radiating elements to the thermal tile interface, or 
even beyond, by recessing the tile. This radical departure from the pres- 
ent conventional design will indeed improve the axial ratio if properly 
implemented. 

The radiating-element design was addressed next. The existing 
crossed-dipole design was found to be the most suitable because the cir- 
cular polarization characteristic could be readily controlled and the 
criterion of minimal modification was satisfied. Shaping of the radome, 
as described earlier, was deemed unfeasible and the related concept of a 
dielectric rod antenna was similarly dismissed due to the difficulty in 
controlling the propagation modes and, therefore, the radiation properties. 

The "wedge 11 dipole concept, basically a curved dipole as sketched 
with dashed lines on the present dipole configuration in Figure 1, was then 
introduced. This design attempts to broaden the beam of the individual 
crossed-dipole elements and, therefore, the array beamwidth by physically 
altering the radiation pattern of the dipole by geometric means. Gener- 
ally, dipoles are not considered in a three-dimensional manner. The idea 
here is to exploit the maximum radiation intensity normal to the radiating 
current by imposing a symmetric angular relationship to the crossed-dipole 
elements, thereby broadening the beam. The basis for this idea is that 
the peak gain of the array is approximately 6.5 dBci, as illustrated in 
Figure 2, and this excess gain can be usefully employed by spreading out 
the energy by broadening the beam, as depicted in Figure 3. Any increase 
in the beamwidth is substantial with regard to the percentage coverage 
since the area about the perimeter increases rapidly for relatively small 
beamwidth additions. 

The iris aperture restriction, which earlier broadened the beam 
by decreasing the aperture area, was discussed. It was felt that this 
planar iris would affect the wege dipole design by maintaining the same 
aperture fields which determine the final radiation pattern. Incorpora- 
tion of the protruding-antenna concept, however, permitted extension of 
the wedge dipole past the circular iris into the radome, as sketched in 
Figure 4, which appeared to be a reasonable compromise of the conflicting 
requirements. Therefore, the basis of the modification program was 
conceived. 
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Figure 1. Present Dipole-Element Configuration 
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Figure 2. Excess Radiation (> 4 dBci) to be Redistributed 
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Figure 3. Redistributed Radiation by Beam Broadening 
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Figure 4 


Proposed Wedge-Shaped Protruding Dipole 
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The angular relationships of the wedge dipole cannot be easily 
assessed but, using the previous gain values and standard dipole radiation 
equations, an angle of 30° is free space is not unreasonable. If too large 
a wedge angle is used, it is conceivable that the peak gain would be re- 
duced below 4 dBci, which is counterproductive. However, an R&D effort 
could establish the optimum angle using this baseline design wedge-dipole 
concept once it was demonstrated. The effects of the radome, however, 
must be considered since some refraction might occur. The angular inter- 
face, using geometric optics, now has a tendency to further broaden the 
beam and, thus, this lens phenomena must be compensated for by reducing 
the free-space baseline angle. Therefore, an angle of 20° is not unrea- 
sonable for any initial R&D attempt to verify the concept. 

The limitation on the percentage coverage measurement is that 
only JSC has the facilities to measure and analyze the data. Rockwell, 
however, can measure the principal plane dipole-element patterns, then 
analytically use pattern multiplication to predict the array beamwidths. 
This assumption allows Rockwell to assess the modification results in a 
preliminary manner to determine if the percentage coverage was indeed 
improved. 

Since the dipole shape would be modified, the other recommenda- 
tion made to Rockwell was to correct what was felt to be a design defi- 
ciency in the shorting strap of the crossed dipoles. There was a failure 
in the connecting solder joint during earlier vibration acceptance tests 
which was attributed to a "blind" connection where inspection was not 
feasible. By separating the crossed shorting straps onto two separate 
discs, visual inspection was made possible, avoiding a recurrence of the 
problem. However, it was noted that the positions of the shorting straps 
were above the plane of the crossed dipoles, which creates an uneven radi- 
ation pattern since the radiating current is "stepped," especially in the 
critical central region where the current density is highest. Although 
it was felt that the shorting strap functioned only as a bridge to con- 
nect the dipole elements, it would, by definition, affect the radiation 
pattern since it carries current. Therefore, it was suggested that the 
dipole connection points be indented, as shown in Figure 5, to accommodate 
the shorting strap on the same level as the dipoles to retain the classi- 
cal dipole configuration. The wedge design of the dipole complicates this 
indented dipole strapping arrangement, but it might be modified as shown 
in Figure 6 to incorporate the angular configuration. 
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(a) Present Configuration with Stepped Shorting Strap 



(b) Proposed Indented Shorting Strap 


Figure 5 


Indented Shorting-Strap Concept 
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4.0 CONCLUSIONS 

Although it is somewhat late to introduce new concepts which 
attempt to improve the gam percentage coverage performance of the Orbi- 
ter S-band quad antennas, it is still feasible in that the modifications 
required to upgrade the design, once it has been satisfactorily demon- 
strated, has been intentionally restricted to be minimal. At the present 
time, there are definite questions as to the reproducibility of pattern 
coverage results, even though they do not satisfy the existing gam per- 
centage coverage specification. Rather than simply accepting the present 
design and, therefore, the inadequate performance, these possible modifi- 
cations have been introduced which may ameliorate the poor performance. 
Simply tuning to obtain impedance matching for selective frequencies m 
order to reduce losses is often unsatisfactory if large discrepancies in 
predicted performance arise unexpectedly. It should be noted that antenna 
development is not a textbook problem and the antenna performance is often 
dictated by minor modifications to correct obvious deficiencies which 
result from the initial design. 

The wedge-dipole design which has been introduced here will most 
certainly broaden the beamwidth and, therefore, increase the gain percent- 
age coverage. In retrospect, this wedge dipole is a descendent of the 
earlier radome lens concepts since it physically implements the ray bend- 
ing that was to be achieved by refraction effects from a concave lens 
surface. In fact, some design compensation must be incorporated for the 
radome lens effect since it has a tendency to broaden the beam even fur- 
ther It must be recognized that this modification has the specific goal 
of increasing the gain percentage coverage by broadening the existing 
antenna pattern at the cost of reducing the peak gain. If gain was the 
essential criterion, as might be required for some marginal communica- 
tion links, other more complex solutions would be required. 
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There has been a recent emphasis on high-gain antennas for both 
the Shuttle Orbiter S-band quad antennas and the Centaur vehicle. An 
earlier Axiomatix report (R8106-4, June 25, 1981) outlined a two-dimensional 
switched-beam array using five elements. However, since high gam is of 
interest, this report continues the development of larger switched-beam 
arrays. Although the basis of these concepts utilizes electromechanical 
switches, this technique is by no means limited by this constraint since 
PIN diode phase shifters can also be readily substituted. The justification 
for this report is to expand the available options by describing the use 
of switched beams rather than phased arrays for applications where high 
power and no electrical holding current are desired at relatively slow 
switching rates. Since the switches are utilized in parallel, little 
overall loss (0.2 dB for a switch and discrete phase element) occurs. A 
systems trade-off using PIN diode phase shifters can then be made using 
the baseline configurations for comparison purposes. 

The next step in increasing the gain beyond the basic three- 
element linear array described earlier is to add two additional outer 
elements with the appropriate illumination taper. The five-element 
array is sketched in Figure 1 with the additional switching arrangement. 

Note that there is no additional loss compared to the three-element array 
since the switches are incorporated in parallel. The primary difference 
in the new set of switching circuits is that the discrete phase shifts 
are twice those of the corresponding three-element set to account for 
the elements being at the extremes of the array. 

The proposed Centaur bicomcal array is now composed of three 
. elements, but higher gam can be achieved (with no additional circuit loss) 
by adding this second set of switching circuits. The beam position which 
is described by a binary code, as described in an earlier report, is iden- 
tical for both the three-element and five-element arrays since the switch 
positions are the same for both sets of switches. 

This concept of extending the linear array can be continued 
further by adding additional pairs of outer elements but, for explanatory 
purposes, a large planar array will be developed using the five-element 
linear array as a basis. 



ORIGINAL PASS 
OF POOR QUALITY 



Beam 

Element 1 

Element 2 

Element 3 

Element 4 

Element 5 

1 

0 O - 6 0 

00 - 30 

*0 

00 + 30 

00 + 60 

2 

<P 0 - 40 

0 O - 20 

^0 

0 O + 20 

00 + 40 

3 

0 O ” 2 0 

00 " 0 

^0 

00 + 0 

0 O + 20 

4 


^0 

0 O 

0 O 

0Q 

- 

*0 

^0 

0 O 

^0 

00 

5 

0 O + 20 

00+0 

^0 

0 Q “ 0 

0 O - 20 

6 

0 Q + 40 

00 + 20 

*0 


0 O - 40 

7 

00 + 60 

0 O + 30 

*0 

|g 

00 “ 60 


Figure 1 Seyen-Beam/Five-Element Svntched-Beam Array 
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The last planar array described was a five-element array sketched 
in Figure 2. Note that the central element is directly fed, and the four 
outer elements have circuit losses of 0.6 dB (0.2 dB per switch) for a 
49-beam (i.e., 7x7) configuration with redundant broadside beams. A 225- 
beam (i.e., 15x15) configuration would have only 0.8 dB of loss bbt, for 
the present, let us only use three switches in parallel as the baseline 
design. 

There is a definite trade-off between the gain desired and the 
number of beam positions since these two parameters constitute the gam 
percentage coverage. Naturally, the higher the gam, the more complex 
the system, but the number of beam positions can be approximately doubled 
by adding only another switch (and 0.2 dB of loss) to each set of switches 
This feature is especially important when the gam percentage coverage is 
emphasized since, effectively, the percentage coverage is directly related 
to beam positions 

To continue with the large planar array development, if two of 
the five-element linear arrays are orthogonally positioned, a nine-element 
planar array with extremely high gain is developed by the two-dimensional 
pattern multiplication concept. Again, note that no additional circuit 
losses are involved since the radiated power is redistributed and the 
switching sets are in parallel The essential features of this nine- 
element array are sketched m Figure 3 Again, the beam position is deter- 
mined by a binary switching code which, for this two-dimensional case with 
three switches in each arm, is composed of six binary numbers for the 49 
beam positions (assuming redundant broadside beams). The 225-beam config- 
uration is correspondingly characterized by eight binary numbers. 

The next step in increasing the gam involves filling in the array, 
as shown in Figure 4, with another four elements. This 13-element configur- 
ation still has 49 beam positions with no additional circuit loss even though 
four more elements are incorporated The addition of these four additional 
elements, however, greatly increases the complexity since these elements are 
shared by the two orthogonal linear arrays so that the appropriate phase 
shifts must be computed by an algorithm. For example, if the phase shifts 
for a given beam position are characterized by Q>q ± m$0Q>g ± nt{0, the Pf> ase 
shifts for these four elements are then described by £<|>q ± (±m±n)<i>j so that 
the complexity is substantially increased. 
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Figure 4 49-Beam/l 3-El ement Switcfied-Beam Configuration 
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Note that the set of switches required to fill in the array is 
one more than that used for the other sets; in this case, four switches 
instead of three because of the (m+n) factor which doubles the number of 
phase shift increments required. For example, for beam position A1B1 
which, for element 2 of both sets, is Qj>g - 3^3C^g “ mutually ad- 

jacent fill-in element should then have a phase shift of O’g'^'O’ which 
requires four switches for integral values of Because of the added 
complexity, however, a system analysis is required to justify filling in 
the array to obtain the higher gam values. 

In summary, techniques to implement a large switched-beam planar 
array have been described. Higher gain has been achieved by adding more 
switches and elements to Increase the effective aperture. The circuit 
losses remain constant except for any incurred for redistributing the power, 
which should be small. The primary value of this approach is the large 
number of beam positions available, which approximately doubles for every 
switch added to a set, at a cost of only 0.2 dB of additional loss per 
switch and discrete phase element. For two orthogonal three-switch sets, 
for example, 49 beams are achievable. Even more dramatic are two ortho- 
gonal four-switch sets which can generate 225 beam positions. With the 
capability of precise pointing, the peak gain of the array can be fully 
utilized and the net gam percentage coverage is thus accordingly higher. 
Another advantage is that the switches are electromechanical latching 
switches which can handle high-power levels and require no electrical hold- 
ing currents. The discrete phase shifts are lengths of cable or microstrip 
sections which are not temperature-dependent. Finally, the selection of a 
particular beam position can be reduced to a six-bi nary-number sequence 
for the three-switch set (49 beams) and an eight-bi nary-number sequence 
- for the four-switch set (225 beams), thereby simplifying the software 
required for implementation. 

These advantages should be compared with comparable electronic 
phased arrays with PIN diode phase shifters to select beam positions. 

This electromechanical switched-beam approach is used on the Shuttle S-band 
quad antennas for two-beam positions, but its capabilities for large planar 
arrays were not fully investigated or appreciated previously. The intent 
of this series of interim reports has been to discuss and develop the 
necessary techniques to exploit this technology as an alternative option 

i 

in selecting future antenna configurations for the high-gain percentage 
coverage applications required for the new generation of communication 
requirements for space vehicles. 
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"A LARGE SWITCHED-BEAM PLANAR ARRAY" 


Further consideration of the 49-beam/13-element switched-beam 
array has resulted in an improved version which completely eliminates 
the software problem of determining the proper phase settings for the 
four "fill-in" elements. The improved configuration is shown in Fig- 
ure 5, where again there are only two switch sets to consider. The mod- 
ification which permits this simplification is feeding the "fill-in" 
elements off of switch set A, which determines the proper phase shifts 
for one orientation. Switch set B then determines the proper phase 
shifts for the orthogonal orientation. The primary advantage of this 
configuration is that the two switch sets designated B have identical 
binary codes so that the beam position notation is greatly simplified. 
Fewer switches (six versus eight) are required but, because of the power 
distribution, additional circuit loss is encountered (although, because 
of the weighting factors, the actual amount is small). 

This improved version of the 13-element array completes the 
basis for the development of a large planar switched-beam array since 
the basic mechanism is clearly outlined. Further modifications for even 
larger arrays can be developed from the techniques outlined, but the con- 
ceptual feasibility has been demonstrated for this new type of array. 




Figure 5. Improved 49-Beam/13-Element Switched-Beam Array 





A VARIABLE-GAIN MULTIPLE-SWITCHED-BEAM-ANTENNA CONFIGURATION 


Interim Report 


Contract No NAS9-16067, Exhibit B 
Technical Monitor William Teasdale 


Prepared for 

NASA Lyndon B. Johnson Space Center 
Houston, Texas 77058 


Prepared by 


Axiomatix 

9841 Airport Blvd., Suite 912 
Los Angeles, California 90045 


Axiomatix Report No. R8209-3 
September 22, 1982 



ORIGINAL P«®» 

OF POOR QUALITY 

1.0 INTRODUCTION 

The difficulty of achieving the original percentage coverage gain 
specification of 85% 4 dBci, which has been relaxed to 50% 4 dBci , became appar- 
ent during the Shuttle Orbiter S-band quad antenna program. The purpose of this 
specification was to ensure that adequate communication link margins were avail- 
able about the Orbiter. As a result of monitoring the S-band antenna program 
under NASA Contract NAS 9-16067, Axiomatix extended many of the concepts used 
in the switched two-beam quad antenna, taking advantage of the Rockwell -developed 
basic hardware. The outcome of the study has been the evolution of a new gener- 
ation of variable-gain multiple-switched-beam antennas whose principles can be 
applied to the improved Shuttle Orbiter S-band quad antennas and Centaur biconi- 
cal arrays, as well as to proposed future Space Station antenna systems. 

This switched two-beam quad antenna developed by Rockwell was an 
improvement over the earlier single-beam quad antenna built by Watkins-Johnson 
but, due to overoptimistic computer analysis projections of the antenna patterns, 
the actual performance results were disappointing. The primary limitation of the 
gain performance was the selection of only two elements and two beam positions, 
the minimum required for switched-beam operation. Using the same type of hard- 
ware configured m a slightly different manner, Axiomatix developed a two-dimen- 
sional cross array that could be electromechanically scanned in both the roll and 
M N 

pitch planes at 2 x 2 discrete beam positions, where M and N are the number of 
double-pole/doubl e-throw electromechanical latching switches in each plane. The 
advantages of using this type of switch rather than electronic phase shifters are 
(1) high power-handling capabilities, (2) low insertion losses, (3) absence of 
electrical "holding" power, and (4) the binary nature of the double-throw switch 
which can readily be converted to a binary code for a specific beam position. 

Higher gain can be achieved by adding additional radiating elements to 
the basic configuration, the only complexity being a doubling of the differential 
phase shifts of the outer elements. The binary code for the beam positions re- 
mains the same since the switching configuration used for the inner elements is 
identical to that of the outer elements. 

Since it is possible to physically switch out various combinations of 
radiating elements to form an antenna pattern of the desired shape using this 
same electromechanical switch, an even more unique feature arises from this 
switched-beam configuration in that low, medium and high gains can be indepen- 
dently attained in each plane by simply activating two switches, with the beam 
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position remaining constant and without phase transients. This variable-gain 
capability in both planes greatly relieves the tracking requirement since the 
pattern shape can be adjusted to accommodate the relative motion between the 
antenna system and the "target" vehicle. For example, if a free flyer was pass- 
ing the Shuttle Orbiter (i.e., being in a lower orbit), the low/high-gain con- 
figuration combination would result in an elongated pattern which follows the 
path of the free flyer, thus eliminating the tracking function. 

The variable-gain mul tiple-switched-beam antenna takes advantage of 
the trade-off between gain and percentage coverage. If high gain is required, 
the full antenna-array aperture is utilized, and this switched-beam array can 
position the narrow beam anywhere m its sector of coverage. If high gain is 
not required, the outer radiating elements are switched out and the larger re- 
sultant beamwidth associated with the lower gain relaxes the tracking require- 
ment. Finally, if low gain from an omnidirectional antenna is needed for a 
nearby target, only the central radiating element is used by switching out the 
inner two elements, which also automatically removes the outer elements, as 
will be shown later. 
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2.0 MULTIPLE BEAMS 

The Rockwell proposal to develop an improved S-band quad antenna 
using the Axiomatix "cross" five-element configuration prompted a more detailed 
evaluation of the design. As a result, the switched-beam phase relationships 
were reexamined and a significantly better scheme evolved, one which was not 
immediately obvious. Since a similar switched-beam scheme has also been pro- 
posed for the Centaur vehicle, documenting the results herein is considered 
appropriate. 

This modified configuration avoids the "double-angle broadside gap" 
previously described in Axiomatix Report No. R8I12-4, where the two adjacent 
near-broadside beams were pointed ±<j> for an angular separation of 2$. In order 
to avoid this problem, it was suggested that the one-beam position be sacri- 
ficed while having a redundant broadside beam. With this improved phase rela- 
tionship scheme, the beams are equispaced and separated by the same angle, <f>' = 
2<j)/3, without losing a beam position. Note that the near-broadside beams are 
now ± <j>/3. The only design criteria is then the selection of the number of beams 
in each plane and, therefore, the number of switches required. 

In order to dramatically demonstrate the desirability of adding 
another switch in series, which doubles the number of beam positions at the cost 
of only 0.2 dB additional loss, the phase relationships for the two-, three- 
and four-switch configurations are shown in Figures 1-3 for four, eight and 16 
beams in each plane, respectively. It should be noted that a similar switching 
arrangement exists in the orthogonal plane so that various combinations of beam 
positions, up to a maximum of 256, are possible for four switches in each plane 
for each of the four quad antennas. Although this might be considered "overkill 
the points to be emphasized are that the simple addition of another switch in 
each circuit greatly increases the potential advantages of switched-beam technol- 
ogy by doubling the number of beams, and the resultant availability of a large 
number of beam positions easily compensates for the unpredictable performance 
which has plagued the quad antenna program due to the curved Orbiter shape. 

The high-gain version of the switched-beam array is shown in Figure 4, 
where two additional elements create a five-element linear array whose phase 
lationships for the 16 beam positions are observed to be similar to the three- 
element array except that the differential phase shifts for the outer elements 
are doubled. The binary code for the beam position is identical since the 
switching configuration for the five-element array is directly related to the 




Figure 1. Two-Switch/Four-Beam Configuration 
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three-element array. Note that, since the additional four switches for the 
outer array are in parallel, negligible additional switching insertion loss oc- 
curs. And, if four more elements are added in the orthogonal plane, high gain 
can be attained in that plane also, and the total number of possible beam posi- 
tions is then 256. 
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3.0 VARIABLE GAIN 

Further progress in the switched-beam-array technology has resulted 
in the conceptual development of a variable-gam antenna which can be adjusted 
for particular mission requirements. For example, on the Centaur five-element 
bicomcal array, there may be occasions when high gain is unnecessary for a dis 
tant target and it is desirable to have a wider beamwidth with medium gain in 
order to simplify tracking. In this case, a scheme to switch out the two outer 
elements has been developed using the same doubl e-pol e/double-throw electrome- 
chanical latching switches, as shown m Figure 5. Further, extension of the 
variable-gam technique applies for the low-gam or omnidirectional mode where 
only the central element radiates. Thus, by adding two switches, it is possibl 
to have low-, medium- and high-gain capabilities with the same array, with very 
little additional insertion loss. 

For the improved version of the Shuttle Orbiter S-band quad antennas, 
this technique can similarly be used for variable-gain capabilities in both di- 
mensions, as shown in Figure 6. Thus, it is possible to alter the pattern to a 
desirable shape, such as using high gain m one dimension and low gam in the 
other in order to form an elongated beam that accommodates a moving target and 
therefore satisfies the tracking function without any antenna switching. 
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TECHNICAL MEMORANDUM NO. M8209-4 


TO: Don Eggers ' DATE: September 30, 1982 

FROM: Richard Iwasaki FILE: NAS 9-16067 "B" 

SUBJECT: Binary Code Designation of Mul tiple-Switched-Beam Positions 


One of the important considerations in implementing a variable-gain 
multi pi e-switched-beam array is the simplicity of obtaining a desired beam po- 
sition. It was noted earlier that the electromechanical double-throw switches 
have a binary nature which can be described by a binary code. For four switches, 
a beam position in one plane can be characterized by a unique four-digit binary 
number and, for the two-dimensional case, by eight binary digits in sequence 
In actuality, only 16 memory locations need to be allocated for this selection 
process since the orthogonal plane can use the identical binary codes. 

As an example. Figure 1 shows the switching arrangement necessary to 
achieve beam position 4. Using the notation that a "0" represents a "straight- 
through" switching position and a "1" denotes a "cross-over" switching position 
and the binary code begins with the switch closest to the radiating elements, 
the binary code for beam position 4 is 1010. Note that the switching arrange- 
ment is identical for the inner and outer elements so that, if the high-gain, 
five-element array is used, no modifications will be required Table 1 lists 
the remaining binary codes for the other beam positions 

For the two-dimensional cross array, the beam position is specified 
first by the X-plane beam position, then by the Y-plane beam position, as de- 
picted in Figure 2, which has a total of 256 beam positions. If the cross- 
hatch beam position designated by 4,15 is desired, the appropriate binary code 
is simply 10100100, where the first four digits characterize beam position 4 
m the X-plane and the last four digits characterize beam position 15 in the 
Y-plane . 

Copies NAS9-16067"B" Distribution 
H. Ellis, RI 
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Table 1. Binary Codes for the 16 Beam Positions 


Beam 

Position 

Binary Code , 

1 

0 

0 

0 

0 

2 

1 

1 

0 

0 

3 

0 

1 

1 

0 

4 

1 

0 

1 

0 

5 

0 

0 

1 

1 

6 

1 

1 

1 

1 

7 

0 

1 

0 

1 

8 

1 

0 

0 

1 

9 

0 

0 

0 

1 

10 

1 

1 

0 

1 

11 

0 

1 

1 

1 

12 

1 

0 

1 

1 

13 

0 

0 

1 

0 

14 

1 

1 

1 

0 

15 

0 

1 

0 

0 

16 

1 

0 

0 

0 


Plane 


Y-Pl ane 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

1 

















2 

















3 

















4 

















5 

















6 

















7 

















8 

















9 

















10 

















11 

















12 

















13 




ME 

001 













14 



101 

ooj 












15 




If 













16 



















Figure 2. 256-Beam-Position Selection Using Binary Codes 
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9841 Airport Boulevard • Suite 912 • Los Angeles, California 90045 * Phone (213) 641-8600 
TECHNICAL MEMORANDUM NO. M8209-5 


TO: Richard Iwasaki DATE: 
FROM. dim Dodds FILE* 
SUBJECT: Simplification of Antenna Beam-Switching Logic 


September 30, 1982 
NAS 9-16067"B" 


The switched-beam array described in Axiomatix Technical Memorandum 
No. M8209-4 requires that the beam position be mapped into a binary code in 
order to control a set of phase-shifting elements. This can be effected by 
using one of the following methods: (1) a microprocessor to access one of 
16 memory locations, each of which contains a four-bit code, (2) a program- 
mable logic array or, (3) hard-wired logic. The format of the conversion from 
four input bits to four output bits is particularly simple in this application 
and can be implemented with the EXCLUSIVE-OR gates, as can be shown by inspec- 
tion, or by using Karnaugh maps. 

Table 1 of M8209-4 is reproduced below, with the ordinal beam posi- 
tions represented by four-bit binary sequences bg, bp bp bg. The binary 
output is represented as a four-bit code ag, a^, a^, ag. 



INPUT 





OUTPUT 


Beam 

Binary 

Code 

Binary Code 

Position 

b 0 

b i 

b 2 

b 3 

a 0 

a l 

a 2 

a 3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

1 

1 

1 

0 

0 
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1 

0 

0 

1 

1 

0 

4 
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0 

1 

1 

1 

0 

1 

0 

5 
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0 

0 

0 

0 

1 

1 
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1 

0 

1 

1 

1 

1 

1 

7 

0 

1 

1 

0 

0 

1 

0 

1 i 

8 

0 

1 

1 

1 

1 

0 

0 

1 

9 
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0 

0 

0 

0 

0 

0 

1 

10 

1 

0 

0 

1 

1 

1 

0 

1 

11 

1 

0 

1 

0 

0 

1 

1 

1 

12 

1 

0 

1 

1 

1 

0 

1 

1 

13 

1 

1 

0 

0 

0 

0 

1 

0 

14 

1 

1 

0 

1 

1 

1 

1 

0 

15 

1 

1 

1 

0 

0 

1 

0 

0 

16 

1 

1 

1 

1 

1 

0 

0 

0 
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The phase shift logic can be implemented as shown below for both the 
X-plane and Y-plane logics, with a Q =b 3 , a 1 = b 2 ©b 3 , a 2 = b 1 ©b 2 and a 3 =b 0 ©b 1 . 
As can be seen, only three mod-2 (EXCLUSIVE-OR) adders are needed per plane. 



Copies NAS 9-16067 "B" Distribution 
Haynes Ellis, RI 
Sam Wong, RI 
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The major thrust of Axiomatix investigations of payload communications 
centered on the payload signal processor (PSP) bit-synchronizer problem. During 
ESTL tests, it was found that the TRW PSP bit synchronizer would false lock under 
high SNR conditions. This pnenomenon was subsequently duplicated at RI, ADL and 
TRW. 

Axiomatix efforts to analyze this problem included numerous conversa- 
tions with Brett Parrish of JSC, Larry Stein and Rich Helgeson of TRW, and Mike 
Rueterman of RI, as well as attendance at meetings and hardware demonstrations. 
Axiomatix has analyzed the algorithm lock-detector- threshold parameter for op- 
timality over SNR's ranging from low to high as well as carefully evaluating 
the entire algorithm. A recommendation was made to change both the threshold 
parameter and the filter bandwidths, computer simulations later indicated an 
improvement in performance with the new parameter values. 

The technical memoranda which follow document the analyses performed 
and meetings attended by Axiomatix personnel on this problem. 



9841 Airport Boulevard • Suite 912 ® Los Angeles, California 90045 ® Phone (213) 641-8600 
TECHNICAL MEMORANDUM NO. M8206-1 

TO: G. Huth DATE: June 18, 1982 

FROM: J. Holmes 

SUBJECT: PSP Bit Synchronizer NRZ False-Lock Status--Interim Report 

REFERENCES: (1) Telephone Conversations with Brett Parrish (NASA/JSC), 

Larry Stein (TRW), and Rich Helgeson (Consultant) 

(2) Abbreviated Flow Chart for the PSP Bit Synchronizer 
Algorithm, by Brett Parrish 

1.0 INTRODUCTION 

During recent ESTL tests, it was found that the TRW PSP bit synchronizer 
would false lock under high SNR conditions. This phenomenon was subsequently 
duplicated at RI, ADL and TRW Axiomatix has begun to analyze the problem. 

This effort has included numerous conversations with Brett Parrish of JSC, Larry 
Stein and Rich Helgeson of TRW, and Mike Reuterman of RI, as well as attendance 
at meetings and hardware demonstrations. Axiomatix is analyzing the algorithm 
threshold parameter for optimality over SNR's ranging from low to high as well 
as carefully evaluating the entire algorithm. This interim report documents 
the findings to date concerning this problem. 

2.0 STATUS SUMMARY 

(1) Brett Parrish has determined that the summation used in the 
filter computation of the mid-phase filter* overflows at high SNR such that 
the algebraic sign of the phase update control changes and therefore provides 
a stable lock point at nominally 180° away from the true-lock point (a half- 
bit error) Brett has also demonstrated at JSC that, by prohibiting the over- 
flow and associated algebraic sign change, the false-lock condition at -180° 
will not occur This fact has been verified by Rich Helgeson with the TRW 
program simulation 


-?r 


FSTB2+PRTB2 doesn't have overflow protection on page 2 of Brett 
brief algorithm. 


's 
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(2) Rich Helgeson and others have determined that a "lack-of-a- 
transition" detector must be added (on page 2 of the brief flow chart). This 
is because, for example, if the data were turned off, no transitions would 
occur so that it is impossible for the lock detector to enter the out-of-lock 
state and, therefore, to return to the acquisition mode. According to Larry 
Stein, this detector, which counts nontransitions, can be implemented without 
a great deal of trouble. Most probably, it is also a problem when the subcar- 
rier or the carrier drops out temporarily. 

(3) The offset in both the true-lock and false-lock points can be 
attributed to the fact that, at high SNR, the A/D converter outputs its maxi- 
mum value of +3 (in binary-twos complement) so that 64 samples of the maximum 
value limit at about two thirds of the way across a presum, or about 43 samples. 
The remaining samples are, in effect, discarded. Although this problem is dis- 
concerting, the author does not believe that it must be fixed immediately, if 
at all, since the main effect is to induce a static phase (or timing) error at 
either the true-lock or false-lock point. Larry Stein has indicated that it 

is possible to fix the problem by allowing only 16 samples per presum to ac- 
cumulate, then scaling the sum of the presums so that it will not overflow in 
an eight-bit register 

(4) Manchester data does not produce false phase lock Possibly, 
this is because, with NRZ data and a -180° error offset, the sum of the mid- 
phase signal can roll over and change sign, providing a stable ~180° error 
whereas, in the Manchester case, there is always a transition at (true) mid- 
phase so that, at least for transition densities not equal to zero or 100&, 
Manchester data is tested for -180° false lock in the algorithm and would be 
rejected 

(5) The- automatic gain control (AGC) used for the A/D converter is 
done in a 5-MHz bandwidth, but the input is then filtered by a 250-kHz filter 
which makes this signal poorly AGC'd since the signal component will increase 
at high SNR's in such a way as to produce a level larger than the A/D range. 

(6) It appears that the "very small timing error" test performed in 
the bit synchronizer algorithm (which is indicated on page 3 of Brett's short 
algorithm just below the NTSPC transition count from the last phase correction) 
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is not required and could be eliminated if necessary in order to implement the 
changes required in items (1) and (2). This check bypasses a phase update cor- 
rection if the error is small enough; otherwise, it does nothing. 

2.0 CONCLUSIONS 

It is the opinion of the author (and others) that it is important to 
fix only items (1) and (2) in order that the bit synchronizer will not exhibit 
false lock, namely, change the mid-phase summation of FSTB2 and PRTB2 so that 
it will not roll over to the opposite algebraic sign, then add a nontransition 
counter to limit the number of nontransitions which can occur before the algor- 
ithm declares that the loop is out of lock Neither "fix" appears to be a major 
problem. It would be desirable to fix item (3), the saturation in the presum, 
but this item is not related to false-lock avoidance. Also, item (5) is not 
an important item and probably should not be changed if its space in the micro- 
processor is not needed. 

Any additional comments, observations or corrections should be indi- 
cated to the author so that a revised memorandum can be written prior to the 
next meeting at TRW scheduled for early July. 


Distribution * 

R. Helgeson, Consultant 
J. Johnson, NASA/JSC 

S. May, TRW 

B. Parrish, NASA/JSC 

T. Pederson, TRW 

M Reuterman, Rockwell 
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TECHNICAL MEMORANDUM NO. M8206-4 


TO: G. Huth DATE: June 30, 1982 

FROM: J. Holmes 

SUBJECT: PSP NRZ Lock Detector and Loop Control Signal 

Phase Error Range of Operation at High SNR — Interim Report 


1.0 INTRODUCTION 

This is the second technical memorandum/interim report that concerns 
the PSP bit synchronizer false-lock problem. The causes of the false lock are 
considered based on analyzing the bit synchronizer and lock detector, but not 
considering the DC bias removal algorithm It appears that, at least at high 
SNR, the lock detector constant (6.5) should be increased. In order to enable 
us to choose the best value of the lock detector constants, an analysis is con- 
tinuing and documentation will be forthcoming. 

2.0 SUMMARY 

The error control signal and the lock detector test for the PSP S-band 
bit synchronizer has been analyzed at high SNR for the 8-kbps data rate case by 
hypothesizing various fixed phase errors and evaluating the inphase and midphase 
filter outputs. It is to be noted that DC bias removal was not considered in 
this analysis. 

It was found that the lock detector indicates an in-lock condition 
anywhere within a phase error of zero to about ±75° (±21% of a bit, which is 
considered too large) Furthermore, the lock detector will indicate an in-lock 
condition between about 104° and somewhat less than 130° (and between -104° and 
— 130°) and also between 150° and 210°. 

The error control signal of the bit synchronizer loop pushes the loop 
in the wrong direction when the phase error is between 90° and 150°, which is 
due to the unprotected overflow of the midphase addition . 

This control signal reversal causes the false-lock condition The 
m-lock indication between a phase error of 150° and 210° can be removed by pro- 
hibiting the saturation of the presums. 

It is expected that the above two problems will occur at the other 
data rates as well, although they have not yet been analyzed in detail. 
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3.0 RECOMMENDATIONS 

Based on this preliminary study, it is felt that, at least at high 
SNR, the lock detector constant of 6.5 is too small and should be increased 
since this value allows large timing errors (±21%) to be declared in lock. 

It might be desirable to provide different constants for different data rates 
as well as for low or high SNR conditions. Further analysis aimed at recom- 
mendmg values for this constant is currently underway. 

Also, overflow protection must be provided for the midphase summation 
and the presums should be modified to prevent limiting. 

4.0 ANALYSIS 

Vie consider the lock detector test given in [1] as 

AVIN > 6 .5| AVMI j (1) 

where AVIN and AVMI are lowpass filtered versions of the inphase and midphase 
summers, respectively. Basically, the idea of the lock detector test given by 
(1) is to announce when the PSP bit synchronizer loop is in or out of lock. 

In Figure 1, a simplified block diagram illustrates some details of 
the bit synchronizer loop. First note that all presums and adders outside the 
microprocessor are combined with two's -complement arithmetic with overflow pro- 
tection. Internally, the mphase (AVIN) signal is formed with two's-complement 
combining with overflow protection, whereas the midphase combining is done with 
two's-complement combining without overflow protection 1 

The loop uses a filtered version of the output error to correct the 
timing of the loop, with the addition of a rate-aiding component which attempts 
to remove small doppler errors (rate tracking). 

In order to analyze the lock detector, it is necessary to analyze the 
value of AVIN, AVMI and TRIND, which is the variable that indicates whether a 
positive, negative or lack of transition has just occurred We now consider 
various phase error angles Table 1 lists the data rate, samples per presum, 
fractions of a bit for a presum, samples per bit, and whether or not the number 
of samples saturates the presum. 





Figure 1. PSP Data Transition Tracking Loop for NRZ Data 
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Table 1. Samples Per Bit Versus Data Rate Matrix 


R b 

Samples/ 

Sum 

Fraction of 
a Bit 

Samples/ 

Bit 

Presum 

Saturated? 

16 kbps 

16 

1/4 

64 

No 

8 kbps 

32 

1/4 

128 

No 

4 kbps 

64 

1/4 

256 

Yes 

2 kbps 

64 

1/8 

512 

Yes 

1 kbps 

64 

1/16 

1024 

Yes 


4 1 <)> = 0°, 8 kbps 

Denote 0 as two 's-complement addition {or subtraction) with overflow 
protection and 0 as two 's-complement addition without overflow protection. 
Based on Figure 2 and Table 1, we can write 


AVlN (n+1 ) = LjAVIN^ + ] 96 0 96 J 


( 2 ) 


since the two's complement has a maximum of 127 and a minimum of -128. There- 
fore , 

n 

AVIN C n+1 ) = 127 2 l] (3) 

1=0 1 

Solving (3) yields 

(4) 

For large n, AVIN^ converges to 


AVIN 


(n+1) 


127 


n+1 


1 - L, 


1 - L, 


AVIN^ 






Figure 2 


Phase Error Relationship for $=0° (Zero Timing Error) 
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For the lowpass filter defined by (2), Lj = 0.99 so that 

AVIN^ 12700 or 12800 (5) 


or an average of 


AVIN^ = 12750 (6) 

The bit transition decision variable (TRIND) takes on the value of 1, -1 or 0 
according to whether the bit goes from negative to positive, positive to nega- 
tive, or does not change. Now consider the midphase integration (AVMI) for the 
situation depicted in Figure 2. We have (TRIND = 1) 

AVMI^ n+1 ^ = L m AVMI^ n ^ + | FSTB2 + LSTB2 | (7) 

Now, when centered at zero phase error, the value of FSTB2 = 127 and PRTB2 = -128 
so that we have 


AVMI^ n+1 ^ = L m AVMI ^ = (-128 0 127) 


Solving (8) yields (-128 0 127 = -1) 


AVMI 


(n) . . t 1 - QJ) L „ 0 


95 


1 - L 


M 


Hence, for large n, we have 


AVMI ^ -*•— — — = -20 
1 - L„ 


( 8 ) 


(9) 


(10) 


The lock detector makes the comparison 

AVIN > | AVMI j 


( 11 ) 
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POOR QUALITY 

which, in our present case, leads to 

•? 

12,750 > 20(6.5) = 130 (12) 

Consequently, the loop is declared in lock. 

4.2 4 = 45°, 8 kbps 

Consider Figure 3, which illustrates a 45° phase error (1/8-bit 
timing error). We find that 

AVIN^ n+1) = Lj AVIN^ + 1 96 E (96E96)j (13) 


or 


AVIN^ n+1 ^ = Lj AVIN + 127 (14) 

where, for convenience, we have assumed a positive bit in the AVIN interval* 
Solving, we obtain 


AVIN^ — * 12700 


(15) 


Since TRIND = 1, we have 


AVMI^ n+1) = L h AVMI^ + [(-960-96) S (-48+48E96)] 
or 

AVMI^ n+1 ^ = L m AVMI (n) - 32 
M 

Solving, we obtain 

flVMl(n) - 640 


(16) 


(17) 


( 18 ) 


~± 


We use this assumption in what follows. 



CO 

IX) 

0 

01 

I 

■Pa 


AVIN 



AVMI 


Figure 3 Phase Error Relationship for = 45° (1/8 of a Bit) 


OO 
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The bit synchronizer lock detector comparison is given by 

7 

12700 > 640 (6.5) = 4160 (19) 

c o the lock detector again declares an In-lock condition. 

4.3 <j> = 67.5° , 8 kbps 

Now consider Figure 4, which illustrates the case when <|> = 67.5°. 

For the inphase component, we have 

AV I N (n +1 ) = Lj AVIN (n) + |Q-72+24) 096) 0(127) | (20) 


or 


AVIN (n+1) = Lj AVIN^ n) + | 127 ] 


(21) 


or solving 


AVIN (n) 12700 


( 22 ) 


The midphase recursion relationship is given by (TRIND = 1) 


AVIN (n+1) = L m AVMI (n) + [-128 0 (-72+24 ) 0 96] 


(23) 


This yields 


AVMI (n+1) = L m AVMI (n) - 80 


(24) 


Solving, we obtain 


AVMI 


(n) -80 

0.05 


= -1600 


(25) 


The lock detector test is given by 

? 

12700 > (1600) (6 .5) = 10,400 (26) 

Hence, the lock detector declares that the bit synchronizer loop is still in 
lock. 




Figure 4 Phase Error Relationship for <j>= 67 5° (0 188 of a Bit) 
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4.4 <j> = 75°, 8 kbps 

Figure 5 illustrates the case in which the phase error is 75°. Again 
it can be shown that 

AVMI^ — v 12700 (27) 

and, further, that TRIND = 1 so that, for the midphase signal, we have 

AVMI^ n+1 ^ = L M AVMI^ + | -128 0Q -80+16 ) E] 96] | (28) 

which can be evaluated to be 

AVMI^ n+1 ^ = L m AVMI ^ - 96 (29) 

Therefore, for large n, we have 

AVMI (n ^ -»• = -1920 (30) 

0.05 

Checking the lock detection equation, we have 

12700 > 1920 (6.5) = 12480 (31) 

which indicates that the lock detector is almost ready to declare the bit syn- 
chronizer out of lock 1 It can be shown that the out-of-lock condition for a 
positive transition occurs at 76° 

4.5 = 90°, 8 kbps 

Now we consider the case when $ = 90° , as shown in Figure 6. It can 
again be shown that 

AVIM^ — > 12700 (32) 

so we now consider the midphase component We have 

AVMI^ n+1 ^ = L m AVMI (n) + [(-96 0-96) E (-96 096)] (33) 



AVMI 


~o © 
o 3 

O 353 
S3 p 

O 13 

C > 


Figure 5 Phase Error Relationship for <6=75° (0.21 of a Bit) 


AVMI 


Figure 6 Phase Error Relationship for <j> = 90° (1/4 of a Bit) 


CjO 


ORIGINAL PAGE SI 
OF POOR QUALITY 


M8206-4 


14 


so that 


AVMI (n+1) = L m AVMI (n) + [-128] 

Solving, we find that 

AVMI (n) = -2560 

0.05 


Checking the lock detector function, we obtain 

1 

12700 > (2560) (6.5) = 16640 


and the lock detector still declares that the loop is out of lock. 

4.6 <p = 135°, 8 kbps 

We now consider the case when <j) = 135°. Consider Figure 7. We 
that (TRIND = 1) 

AVIN^ n+1 ) = Lj AVIN^ + | (-96 00)0(96 096) | 


or 


AVIN^ n+1 ) 


Lj AVIN 


(n) 


+ 31 


Solving, we obtain 


AVIN^ -+• +3100 


Now AVIN^ satisfies 

AVMI^ n+1 ^ = L M AVMI^ n ^ + C(-960-96) 0 (-9600)3 


or 


AVMI 


(n+1) 


L m AVMI (n) + (-128 0-96) 


(34) 


(35) 


(36) 


have 

(37) 


(38) 


(39) 


(40) 


(41) 



o G 


rO ”0 

C 35 * 


Figure 7 Phase Error Relationship for <j>= 135° (0 375 of a Bit) 
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Simplifying, we have 


AVMI (n+1) = l H AVMI^ + 32 (42) 

Solving (42) produces -128 0-96 = +32 since the subtraction is not protected. 
In the limit, we have 


AVK l<"> * ^ = 640 


(43) 


The lock detector test yields 

? 

3100 > 6.5 (640) = 4160 (44) 

which implies that the loop is declared out of lock at <j> = 135°! Additionally, 
the error control increases the phase error since the algebraic sign of the 
correction is now positive. In fact, a little reflection will indicate that 
the error signal is positive between 90 +o and 135°; furthermore, the lock de- 
tector is satisfied in a region between 104° and -130°. 

4.7 j - 145°, 8 kbps 

Consider the case when <j> = 145°, which is illustrated in Figure 8. 
Considering AVIN, we have that 


A Vim(h +1 ) = Lj AVIN^ 


( -96 E (96 0(96 096) 


(45) 


or 


AVIN 


(n+1) _ 


= L t AVIN 


(n) 


+ 1-117 + 1271 


(46) 


so that 


AVIN 


(n+1) 


L r AVIN^ + 10 


(47) 


Solving, ,we obtain 



AVIN 



AVNI 


Figure 8. Phase Error Relationship for <j> 
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AVIN 


(n) 


10 

0.01 


= 1000 


(48) 


Since TRIND = 1, we have 

AVMI^ n+1) = L m AVMI (n ^ + ^(-96 0-96) 0 (-96 0-96 (49) 


Evaluating, we obtain 


AVMI^ n+1 ^ + l_ M AVMI ^ + 11 (50) 

so that 

AVH! tn UiL = 220 (51) 

The lock detector test produces 

1000 > (220) (6 5) = 1430 (52) 

which indicates that the loop is out of lock. 

4.8 <j> = 150°, 8 kbps 

In the same manner as that computed above, we find that 

AVIN^ 100 , TRIND = 0 ; AVHI* n * -*■ 0 (53) 

so that there is no update signal, however, the lock detector signal is satis- 

fied and declares that the loop is in lock. 

4.9 <j> = 158°, 8 kbps 

In the same manner, we can show that 

AVIN^ 100 ; TRIND = 0 , AVMI^ 0 (54) 

and again there is no update signal and the lock detector test declares a 

locked condition. It appears that AVMI^=0 out to <{> = 180° and beyond to 210° 
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beyond which the lock detector Is not satisfied, but the control is to push 
the error back to 4) = 180°, where the lock detector again declares that the 
loop is in lock. This would explain a false-lock mode in both the conerent 
and noncoherent tests. The coherent test would yield a very stable false lock 
whereas the noncoherent test would not appear to be as stable. 

5.0 CONCLUSIONS 

We have considered the case in which the data rate is 8 kbps and the 
input SNR is very large; therefore, the signal level into the A/D converter 
yields a signal which is either +3 or -3 (seven-level output). The result 
found in this case was that the lock detector indicates an in-lock condition 
up to a phase error of ±75° (±20$ of a bit). Furthermore, the lock detector 
will indicate an in-lock condition between about 104° and less than 130°. 

Again, in the region 150° to 210°, the lock detector declares an in-lock 
condition. 

It was further determined that the error control signal changes its 
algebraic sign when the phase error is between 90° and 150°, so that phase false 
lock potential exists in the region between -150° and -210°. 

Based on the error control and lock detector functions, it appears 
that false lock could occur if the signal was turned off, then turned back on 
when the phase error was between 90° and 250°. By symmetry, the loop is also 
pushed in the wrong direction between -90° and -135°. 

By preventing the midphase adder from overflowing, it appears that 
the loop will no longer produce a control signal in the wrong direction. The 
only remaining point is that, since the lock detector will be satisfied only 
in the region between 150° to 210°, drift would tend to pull the loop in. But, 
by removing saturation this lock detector false indication would be prohibited, 
thus causing the loop to enter the acquisition mode and hasten acquisition. 
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TECHNICAL MEMORANDUM NO. M8207-2 


TO: B. Parrish, J. Johnson DATE: July 23, 1982 

FROM: J. Holmes, P. Nil sen 

SUBJECT: Recommendations on PSP Bit Synchronizer Design Changes 


Axiomatix has been involved in a review and analysis of the PSP bit 
synchronizer false-lock problem. This has entailed the following effort: 

(1) Attending technical discussions, meetings/telecons 

(2) Witnessing laboratory tests 

(3) Analyzing algorithms and hardware performance 

Axiomatix has not yet completed this activity, especially with respect to (3) 
above However, Axiomatix does feel it appropriate to make certain design 
change recommendations in light of its findings to date. The design changes 
which Axiomatix recommends at this time are considered by Axiomatix as the 
minimum changes required m order to achieve a viable, satisfactorily per- 
forming payload communication capability. These changes are as follows 

(a) Correct the microprocessor midphase addition overflow (and 
sign change) 

(b) Modify or reparameterize the lock dectector algorithm 

(c) Add a bit transition counter. 

Without these changes, there appears to be a reasonable probability 
that, given certain operational situations which result in signal strength 
fluctuations, i.e., antenna pattern modulation of signal amplitude, the bit 
synchronizer will false lock and data will be lost. Furthermore, the present 
lock detector does not reliably indicate this condition. Also, the present 
PSP design may necessitate the SSO crew reinitializing the PSP several times 
in order to achieve initial acquisition. At present, the only reliable indi- 
cation of successful PSP bit synchronizer acquisition is the successful frame 
synchronizer indication. 
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In addition to the three design issues C( a )"( c )3 cited previously, 
there are several other related design issues which have surfaced during the 
investigation which have the potential to affect communication performance. 
These are: 

(d) Saturation of the PSP A/D converter due to a gain change made 
between the 250-kHz filter and the A/D converter 

(e) Saturation of the presum integrators 

(f) Potential DC bias-removal algorithm problems. 

Axiomatix has analyzed the second item (e) and feels that the present presum 
design, while being undesirable, does not cause the false-lock problem to occur 
however, correction of this problem is a relatively simple software change. 
Because the DC bias-removal algorithm (f) has been substantially changed, 
we feel that this should be reviewed by all interested parties in order to 
verify the accuracy of the new algorithm. We also believe that the gain be- 
tween the 250-kHz lowpass filter and the A/D converter (d) should be reduced 
to the original design value or changed to a new, more nearly optimum value 
in order to avoid clipping the signal plus noise. 





Axiomatix 
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TECHNICAL MEMORANDUM NO. M8208-1 
TO: P. Nil sen 

FROM: J. K. Holmes 

SUBJECT: Observations of NRZ and Biphase-L 

PSP False Lock at Rockwell 


DATE. August 9, 1982 

DISTRIBUTION: R. Helgeson 

G. Huth 
J. Johnson 
S. Mayo 
B. Parrish 
M. Reutermann 


Both NRZ and biphase-L false lock induced by RF power interrupts were 
observed at the PSP tests conducted at Rockwell ADL. Details of this observa- 
tion are provided herein. 

Testing of the PSP false-lock phenomenon was observed at Rockwell 
(Downey) on August 2-3 by the author, representing Axiomatix, Mike Reutermann, 
Michelle Knowlden, Pete Berndsen and Barry Mathias of Rockwell (who actually 
performed the tests), Sal Mayo of TRW, and Jack Rivers and Bob Vermillion of 
NASA 

Both NRZ and biphase-L false lock was observed at high SNR (-60 to 
-80 dBm). Three types of simulated RF dropouts were used, as follows: 

(1) Fast switching between two antenna ports with negligible 
differential signal level, except for the approximate 3 ms switch time 

(2) Slow switching from antenna port to dummy load and back to 
the antenna dwelling on the load for about one second 

(3) Manual attenuation to simulate an antenna fade, with maximum 
attenuation held for about one second 

The last two switching procedures were used since the first did not 
induce false lock. Biphase-L false lock was different from NRZ false lock in 
that it was data-pattern dependent and NRZ was not Also, both lock detector 
and frame synchronization discretes are low in biphase-L false lock. Practi- 
cally speaking, since data would tend to be more or less random, the biphase-L 
false-lock problem may not be "real" although it should be verified with actual 
payload data. Table 1 summarizes the differences between NRZ and biphase-L 
false lock. Some of the test data are included at the end of this memorandum. 
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Table 1. False-Lock Summary 


NRZ 

1 

1. Only at high SNR (-60 to -80 dBm) 
(At 16 kbps, need > -80 dBm) 

2. Not data-pattern dependent 
(ability to false lock) 

3. In fal se lock: 

Bit synchronization - high 
Frame synchronization - low 

4. Changing data pattern does not 
break fal se lock 

5. False locks at all data rates 


6. False locks on switching from 
antenna to load to antenna and 
when using the attenuator to 
increase attenuation to simulate 
a fade 

7. Not dependent on having wire on 
PSP clock pin to cause false lock 


Biphase-L 

1. Only at high SNR (-60 to -80 dBm) 
(At 16 kbps, may need > -80 dBm) 

2. Data-pattern dependent (ability 
to false lock) 

3. In false lock: 

Bit synchronization - low 
Frame synchronization - low 

4 Changing data pattern breaks 
false lock 

5 False locks only at 1 and 16 kbps 
(2 kbps not yet tested, does not 

false lock at 4 or 8 kbps) 

6. False locks when switching from 
antenna to load to antenna and 
when increasing the attenuation 
by means of an attenuator to 
simulate a fade 

7. Not dependent on having wire on 
PSP clock pin to cause false lock 
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9841 Airport Boulevard • Suite 912 • Los Angeles, California 90045 • Phone (213) 641-8600 
TECHNICAL MEMORANDUM NO. M8208-2 

TO: P. Mil sen DATE: August 27, 1982 

FROM: J. K. Holmes 

SUBJECT: Recommended Changes to the Tracking Mode of the 

PSP Bit Synchronizer Lock Detector for NRZ Data 

REFERENCE: B. Parrish, "PSP Bit Synchronization Algorithm Flowchart Summary 

for the NRZ Mode," JSC document, undated (about July 1982) 


1.0 INTRODUCTION 

This technical memorandum documents the recommended parameter value 
charges which were previously transmitted to TRW by telephone. The recommen- 
dations are to narrow the m-phase and midphase filter bandwidths and to change 
the lock detector test ratio in the NRZ mode of the PSP bit synchronizer lock 
detector By making these changes, it is believed that the problem of the lock 
detector not indicating loss of lock, sometimes not until minutes after the 
signal is removed, will be eliminated. These recommended changes will require 
only that the bit synchronizer algorithm parameters be modified, and no algo- 
rithm or hardware changes will be involved. 

2.0 RECOMMENDED CHANGES 

It is recommended that both the in-phase and midphase filters be re- 
duced in bandwidth so as to increase their respective output signal -to-noise 
ratios (SNR's) and thereby improve the reliability of the lock detector test. 
Specifically, it is recommended that (I) the L-value of the in-phase filter be 
changed from 0.99 to 0.999 (if feasible to implement digitally) so that the 
"gain" goes from 100 to 1000 and, (2) the L-value of the midphase filter be 
changed from 0.95 to 0.995 so that this "gain" goes from 20 to 200. With the 
narrow bandwidth digital filters, it is possible to change the in-lock bit 
synchronizer test from 

AVIN ;> 6 5 1 AVMI | 

AVIN > 50 j AVM I j 


to 
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The result of this change would be to greatly improve the probability 
of declaring that the bit synchronization loop is out of lock when the signal 
is removed. Presently, the loop can display an in-lock condition, sometimes 
until minutes after the signal has been removed. This would be confusing to 
flight operation personnel. 

The analysis which led up to the recommendations of narrowing the 
digital filters and changing the m-lock test is based on a rather involved 
analysis, which will be the subject of a future memorandum. 

DISTRIBUTION 

R. Helgeson, Consultant 
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TECHNICAL 

MEMORANDUM NO. M8209-2 

DATE: 

September 3, 1982 

TO: 

P. Nilsen 

COPIES: 

NAS 9-16067"B" File 
R. Helgeson, Consult. 

FROM' 

J. K. Holmes 


J. Johnson, JSC 
S. Mayo, TRW 

SUBJECT: 

PSP Bit Synchronizer Lock Detector 


B. Parrish, JSC 


Parameter Value Changes for the 


T. Pederson, TRW 


NRZ Data-Tracking Mode 


M. Reutermann, RI 

REF. (I) 

Technical Memorandum No. M8208-2 



(2) 

B. Parrish, "PSP Bit Synchronization Algorithm Flowchart Summary 
for the NRZ Mode," JSC document, undated (about July 1982) 


1.0 INTRODUCTION 

This memorandum makes a "final" recommendation for the PSP bit syn- 
chronizer lock detector parameters m the NRZ data mode based on the current 
state of the microprocessor available memory and processing rate. Due to the 
lack of additional memory storage space. Rich Helgeson has pointed out that, 
with the present schedule and memory constraints, it is feasible only to repar- 
ameterize the lock detector algorithm at this time. 

The recommended recursive L parameters of the respective one-pole 
in-phase and midphase filters are given by 

254 

L I " = °' 9922 

L, = 248 = 0.9688 
M 25F 

and the lock detector test becomes 

AVIN > 13*|AVMI | 

for NRZ data. 

2.0 BASIS OF NEW RECOMMENDATION 

Based on a recent telecon with Rich Helgeson, it was determined that 
it was impossible to implement the desired values of L = 0.999 for the in-phase 
filter and L = 0.995 for the midphase filter of the bit synchronizer lock detector 
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recommended in ref. 1. Basically, the allowable values of the recursive filter 
parameter L must be of the form 


where n is an even integer less than 256. Based on this constraint, the maxi- 
mum recursive filter L value is then L * 254/256 = 0.9922 and, hence, 0.999 cannot 
be achieved. 

In ref. 1, it was recommended that the lock detector for the NRZ case 
have separate in-phase and midphase lock detector filters which are much nar- 
rower than the bit synchronizer loop in-phase and midphase filters; however, 
in light of the above-mentioned limitations, it is recommended that both the 
loop and the lock detector use the same filters. Thus, the in-phase and mid- 
phase filters should be modified from 


L 


I 


L 


M 


252 

256 

242 

256 


0.984375 

0.9453 


with the lock detector test m tracking specified by 

-? 

AVIN > 6 5* | AVMI | 


to 


L 


I 


L 


M 


254 

256 

248 

256 


0.9922 

0.9688 


with the lock detector test in tracking now being specified by 

■? 

AVIN > 13*| AVMI | 


where the coefficient 13 is close to the desired value of 15 and is a simple 
shift of one digit, resulting in a multi ply-by-two of the current value of 6.5. 

The above changes should provide the same probability of detection 
value when in lock as those of the existing parameters but, in addition, they 
will decrease the probability that the lock detector declares that the signal 
is present when it is not. As a result, the lock detector light should not 
remain on for long periods of time when the signal drops out; this is a problem 
in the original algorithm. 
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TECHNICAL MEMORANDUM NO. M8210-1 

TO: P. Nil sen DATE: October 4, 1982 

FROM: 0. Holmes FILE: NAS 9-16067 ,, B" 

SUBJECT: PSP Lock Detector Performance in the Tracking Mode for NRZ and 

Biphase Data 

1.0 SUMMARY 

In this technical memorandum, we analyze the lock detector performance 
of the Shuttle PSP (in tracking only) and recommend improvements and/or modifi- 
cations to minimize the time required to indicate out-of-lock status when the 
input signal is removed. Because of microprocessor limitations in both memory 
size and speed, it was determined that the digital filters used for the inphase 
and midphase bit synchronizer tracking loop and lock detector must be the same, 
forcing a compromise of sorts. Consequently, with this constraint and the fact 
that only a limited number of one-pole digital -filter-recursion parameters were 
available in the microcode, the following results were obtained. 

• NRZ Data 

A. Present parameters 

Inphase parameter: Lj = 252/256 = 0.9844 

Midphase parameter: L M = 242/256 = 0.9453 

Lock detector criteria: AVIN > 6.5 | AVMI | 

B . Recommended parameters 

Inphase parameter. Lj = 254/256 = 0.9922 
Midphase parameter: L^ = 248/256 = 0.9688 
Lock detector criteria* AVIN > 13 | AVMI J 

® Manchester Data 

A. Present parameters 

Inphase parameter. Lj = 252/256 = 0.9844 
Midphase parameter* L^ = 252/256 = 0.9844 
Lock detector criteria: AVIN > 5 | AVMI j 
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B. Recommended parameters 

Inphase parameters: Lj = 254/256 = 0.9922 
Midphase parameter: L H = 254/256 = 0.9922 
Lock detector criteria: AVIN £ 15 jAVMlj 

The mean time to indicate out-of-lock with the new parameters when the 
signal is removed reduces to one-half the time in the NRZ case and about one- 
tenth the time in the Manchester (biphase-L) case, compared to existing param- 
eters. Since the analysis and system modeling were very involved, approximations 
were made in the analysis; however, it is believed that the answers are represen- 
tative of the actual system performance. 

During the course of this work, numerous helpful discussions were held 
with Rich Helgeson (consultant) and Brett Parrish (NASA/JSC) concerning the op- 
eration of this bit synchronizer.* In addition, observations of the false-! ock 
phenomena on the actual hardware were observed at both Rockwell and TRW. 

2.0 INTRODUCTION 

Shuttle-launched payloads will communicate with the Space Shuttle 
Orbiter payload communications system over a range of tens of nautical miles. 
Command data is transmitted to the payloads while telemetry data are received 
at the Orbiter on multiple S-band frequency channels. Figure 1 illustrates the 
Space Shuttle Orbiter receiver block diagram for payload signals, which is part 
of the Payload Interrogator (PI). The PI sends a 1.024-MHz subcarrier PSK mod- 
ulated by either NRZ or biphase signals at date rates of either 1, 2, 4, 8 or 
16 kbps. Signals from the PI are sent to either the NASA PSP, the DOD PSP or 
the Ku-band signal processor (KuSP). This memo is concerned with the perfor- 
mance of the NASA PSP which performs subcarrier demodulation, bit and frame 
synchronization , as well as matched-filter detection of the data, as shown in 
Figure 2. 

This memo concentrates on PSP bit synchronizer lock detector perfor- 
mance during tracking. Performance specifications of the PSP bit synchronizer 
loop require less than 1.5 dB implementation loss over a signal level corre- 
sponding to an E b /N Q range of 2-11 dB at the minimum average bit transition 
density of 12.5% (64 bit transitions in 512 bit periods). Acquisition must 
occur within four seconds with 0.9 probability for E^/Nq values and data formats 
listed in the specification amendment of March 25, 1982, shown as Appendix D in 
this memo. A more detailed (but still simplified) block diagram of the original 
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bit synchronizer is shown in Figures 3 and 4. Basically, the data is demodulated 
by the 1.024-MHz subcarrier (Costas) loop and is sampled at 1.024 Msamples per 
second, with each sample being quantized to three bits. These samples are stored 
in presums that are added together in eight-bit registers to form the last, mid- 
dle, and first half-bit integrals, as depicted in Figure 5. 

The inphase averaging filter is driven by the magnitude of the sum of 
the first and last half bits, and the midphase filter is driven by the sum of 
the first half bit and the previous last half bit (delayed one bit)— all multi- 
plied by the transition indicator [1,2]. 

3.0 THE THREE-BIT QUANTIZER 

In order to proceed with our analysis, certain assumptions were made 
so that the analysis would be manageable. First, the effect of the second-order 
tracking component was neglected because it was designed for a very limited 
correctional capability. Second, the DC restoration algorithm was not modeled 
because the new algorithm was not available at the time of the analysis for this 
writing and, also, it was felt that, in tracking, the DC restoration was not an 
important factor in the loop operation. In addition, the scaling algorithm was 
not incorporated since it does not have a major effect at low SNR's. Finally, 
it was assumed that the revised algorithm for acquisition and tracking will not 
cause either limiting or register rollover due to unprotected registers This 
has been verified by simulation, calculation and laboratory measurements. An 
abbreviated flowchart is presented in Appendix A which was constructed by Brett 
Parrish, NASA/JSC. 

Table 1 lists the number of samples per bit at each data rate with 
the new microcode. 


Table 1. Samples Per Bit at Each Data Rate 


R b 

Samples/Bit 

Samples/Presum 
(Current System) 

Samples/Presum 
(Proposed System) 

16 

64 

16 

16 

8 

128 

32 

16 

4 

256 

64 

16 

2 

512 

64 

16 

1 

1024 

64 

16 
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Figure 3. PSP Data-Transi tion Tracking Loop and Lock Detector for NRZ Data in the Tracking Mode 
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Figure 4 PSP Data-Transition Tracking Loop and Lock Detector for Biphase Data in the Tracking Mode 
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Figure 5 Time Blocks in One Bit Period for NRZ and Biphase Data 
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Let N denote the number of samples per bit which is either 64, 128, 
256, 512 or 1024, as indicated in Table 1. The sample rate is 1.024 Msps and 
the baseband signal plus noise is filtered to 250 kHz. In Figure 6, a model of 
the deadzone A/D converter is illustrated. Notice that, if the signal were 
less than one-half of a cell in peak amplitude, there would be zero output. 
Consequently, AGC action is necessary to keep the peak amplitude greater than 
one-half of the cell size. The analog baseband signal to the A/D converter is 
given by 


y(t) = A d(t) + n(t) 0 < t < T (1) 

where A is the signal amplitude, d(t) is the baseband NRZ data signal, and n(t) 
is modeled as bandlimited white Gaussian noise. Each bit integral is approxi- 
mated by the sampling of y(t) over one bit time to yield 

N N N 

FSTB2 + LSTB2 = I d. + J, n. + £ (2) 

i=l 1 i=l 1 i=l 1 


Based on the new implementation, we assume that saturation cannot occur in any 
of the registers. In (2), d^ is the ith sample of d(t), n- is the ith sample 
of n(t), and q- is the ith sample of q(t), where q(t) is the quantization noise 
introduced by the A/D sampling process. Although q n is not generally indepen- 
dent of d- and n.. , we will assume that it is m order to simplify the following 
analysis. Now we determine the level of the signal and the thermal noise out 
of the 250-kHz lowpass filter He have 



V 


at) 


(3) 


where B is the lowpass filter noise bandwidth and R is the data rate, R = 1/T, 
with T being the bit duration. The signal power to noise power is then given by 
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For example, at R = 8 kHz and E b /N Q = 4 dB, we have 


SNR = 1.58 ■ - 8 X 10 = -13 dB (5) 

250 x 10 J 


Thus, we tend to have rather negative sample SNR's. Now the sample signal volt- 
age to noise standard deviation is given (from (4)} 


a ^N fl B 



( 6 ) 


When the loop is at the correct timing position, (2) becomes 


FST2 + LSTB2 



N 


+ I 

i=l 


q 


i 


(7) 


the ± sign is determined by the algebraic sign of d(t), te(0,T). Now, since 
the samples are taken at a sample rate of 1.024 Msps and since the lowpass filter 
preceding the sampling is 250 kHz wide, the noise samples are not independent. 
Hence, we assume that sample pairs are independent and each sample has the same 
value in the pair. Using this mode, we have 


FSTB2 + LSTB2 



N/2 

+ 2 l q, 
1=1 1 


( 8 ) 


where, since we have assumed that the noise samples are pairwise independent, 
we make the same assumption for the quantization samples 

Letting n denote the noise sample sum in (8) produces 


N/2 

n = 2 Y n 
i=l 1 


E[n] = 0 


N/2 N/2 r 

Var[n] = 4 ^ ][ E n-,n = 2Na 

i=l j=l L 


(9) 

( 10 ) 


( 11 ) 
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If we let n denote the quantization noise, we have 

H 

N/2 

"q ' 2 "l (12) 

Hence, by symmetry, 

= 0 (13) 

N/2 N/2 r -t 

4 1-1 j-1 E L? 1 q 4i = 2N Var( ^ = ff q (14) 

If we model the quantization noise as having a uniform distribution 
and let L denote the width of the quantization cell, then the quantization noise 
is given by 



and 



E 



= 0 


(15) 


(16) 


If we let the load factor* be denoted by LF, then the optimum LF should be 
about 2.28. From Figure 6 and the definition of the load factor, we have 


or 


2.28 


3L 

o 


L = 0.77 a 


(17) 

(18) 


From (15) , we obtain 



0.05 a 2 


(19) 




The A/D converter load factor is the ratio of peak output to RMS 
noise input. 
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For arbitrary load factors, we obtain 


a q 2 = 0.0092 (LF) 2 a 2 

Using (10), (11), (15) and (16) in 


( 20 ) 


V E h p N/2 N/2 

rj— ~ cx + 2 y n + 2 2 q, (21) 
N 0 B 1-1 1 1=1 1 


yields for the mean and variance 


and 


E[FSTB2 + LSTB2] 



( 22 ) 


Var[FSTB2 + LSTB2] = 2Ncr 2 + 0.1Na 2 = 2.1 N a 2 (23) 


Since the magnitude of the inphase term is recursively filtered, we now inves- 
tigate the statistics of [ FSTB2 + LSTB2J that are needed before we can analyze 
the filtered version. Let x denote FSTB2 + LSTB2 Since the major noise is 

thermal noise, we shall approximate FSTB2 + LSTB2 as a Gaussian random variable 

— o 

with mean y and variance 2, • From [3], it can be shown that 



(24) 


where 


and where 


erf (x) 




2 2 = 2.1 N a 2 


(25) 


(26) 


(27) 
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Now consider the variance of |x|. We have 


Var 


M ■ e 0 ’ 


?I . e 2 


or 


Var 


0*0 ■ 


J . 2 2 

2 + y - 2 


Jl e -» 2 A 2 + V f /_H_\ 

’’ J \M) 


When (y/sr =4 dB, we have 


VarQxQ * £ 2 + ii 2 - (1 .05 w ) 


2 e 2 


Hence, for (y/2/ >4 dB, we have that 


Var 


= Var[x] 


= j E[x] | for > 4 dB 


Hence, if we assume that d(t) is +1 m every bit interval*, we have E 
Now note that 


N = 


R ~ R B 


with R g the same rate and R the data rate. Hence, 


or 



Since R s /B = 4, we have 


a 


* 

Henceforth, we assume that y is always positive and drop the 
magnitude signs for simplicity. 



(28) 

(29) 

(30) 

l>Dc]. 

(3D 

(32) 

(33) 

(34) 
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It is worthwhile to compare the performance of an analog matched 
filter with the digital version that we have just analyzed. The digitized SNR 
is given by 


or 


SNR d 


E 2 [FSTB2 + LSTB2] 
Var (FSTB + LSTB2) 




2 


2.1 N a 2 



(35) 


(36) 


Since the ideal analog matched filter would yield an SNR of 2.0C = Eb/^o3> we see 
that, within our approximations, the sampling and quantization noise have pro- 
duced a loss of 0.2 dB. 


4.0 BIT SYNCHRONIZER DIGITAL FILTER OUTPUT STATISTICS 

WITH SIGNAL PRESENT— NRZ DATA 

We now present some of the mam elements of the bit synchronizer and 
lock detector algorithms. An abbreviated flowchart developed by Brett Parrish 
is included in Appendix B. 

Basically, two one-pole recursive filters are used in the bit syn- 
chronizer and lock detector. The first is AVIN, the inphase bit integral fil- 
ter (Figures 4 and 5); the second is AVMI, the midphase bit integral (Figure 4). 
The former is defined by 


(n+1) 
AVIN = 


in) (n) (n) 

L AVIN + | FSTB2 + LSBT2 j 


(37) 


where n denotes the nth bit time The midphase filter has an output given by 


(n'+l) (n 1 ) / (n * ) (n')\ (n') 

AVMI = L AVMI + ( FSTB2 + PRTB2 jTRIND (38) 


where the n' denotes the n 'th bit transition. Figure 4 illustrates FSTB2, 
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LSTB2 and PRTB2. The term TRIND takes on the value 0, 1 or -1 according to 
whether no bit transition occurs or whether a positive or negative transition 
occurs. Therefore, the last term of (38) is a proportional measure of the bit 
timing error whenever a bit transition occurs. 

To determine if the bit synchronizer should be in the tracking or ac- 
quisition mode, a lock detector is utilized by making the following comparison: 

7 

AVIN >6.5 |AVMlj (39) 


so that, when AVIN is larger than 6.5 AVMI, the lock detector declares that the 
bit synchronizer loop is in lock. 

Consider AVIN first. Using the approximations of (30), we have that 


,ln) 


( n+1 ) ~ W (nl (nl 1 
AVIN = LjAVIN + n in '+n q in; + 1 


a with probability 

f 7 


(40) 


with probability (1 - p t ) 


where p^ is the bit transition probability. The term j is the number of samples 
in which the timing is in error so that j/N is the discrete measure of phase 
error. Note that, when no bit transition occurs, there is no signal loss, how- 
ever, when a transition does occur, the signal component is reduced by the fac- 
tor [l - 2|j|/N}, with j being the number of samples (out of N samples) in 
error. For convenience in the analysis, we shall approximate the discrete timing 
error by a continuous phase error and let 


1 = JL 

N 2tt 


(41) 


where $ is the continuous-valued phase error. Using (41) in (40) leads to 


and 






AVIN^ n+1 * = LjAVIN^ + x (n) 


(42) 


(43) 
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(44) 


where p t is the bit transition probability. Since the phase error is normally 

quite small, we will approximate <j> by a zero-mean Gaussian random variable with 
2 

variance . 

First, we compute We have from (44) that 


n - 


2vfU/sr a 


A _ iLlij] 


"t + ( x -p t ) 


(45) 


For a Gaussian random variable, it can be shown that 


E[|*|3 - V| 


(46) 


Using (46) in (45) produces 


[* (nr I ’ 



2^ +i jj— a 
N 0 


1 _ ' l j 


(47) 


To evaluate the variance of X^, we first compute E 
we have 


(*<■)) 


Using (42), 


(x(n)) 2 ] - E [(n{n) + n q (n)) 2 J + 4E (n (n) + n q (n) ) SK a X (n) 


+ 4 N ~ a 2 E 
IN 0 


(48) 


Evaluating, we obtain (since the second term is equal to zero) 


(x<»>) ? 


- 2.1No 2 +4N i?-o 2 ^ll-- 

7T 


V!^ + \)pt + ( l " p t)| (49) 
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2.1 No 2 + 4N J- a 2 (1 - 

% 


(!) 


3/2 a/ 

V^ p t 

ir 


(50) 


Now 


E 2 P n n = 4N a 2 

L J IN n 


,3/2 


, /2T / " n . 2 2 n 

1 ~ (?) % p t + T % p t 


(51) 


Since 


Var[x (n ^J = E(x (n ))^ - E 2 |^(n)J 


(52) 


we obtain 


Varp] . 2 lN„ 2 + 4^N<, 2 (i) [p t -|p t 2 ], l>p t >0 (53) 


It follows from (A-5) and (A-13) that [0 < p, < 1} 


Var 


AVIN 


o E. o/^.\ 2 r 
2.1 Na + 4 J-Ha M- p t - 
(nn _ N o L Ht 


2 rx 2 

v P t 


1 - L 


(54) 


I 


and 


E AVIN^ 


1 


1 - L, 


(55) 


Now we consider the inphase one-pole digital filter The output of 
this filter is characterized by the following recursion, as indicated in (38) 
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AVMI^ n 1+1 ^ = AVMI^ n ^L M + (FSTB2 


^ n * ^+PRTB2^ n ' ^)TRIND^ n ' ^ (56) 


recall that n denotes the nth bit and n* denotes the n'th bit transition. Only 
where p^ = 1 does n = n'. In the following analysis, we will model TRIND^ n ^ as 
a perfect estimator of the algebraic sign of the bit transition*. Bit transi- 
tions from -1 to +1 produce a value of TRIND = 1 and bit transitions from +1 to 
-1 produce a value of TRIND = -l. Again writing the sum of (FSTB2 + PRTB2)TRIND 
in the form of thermal noise, quantization noise and loop phase error control 
signal, we have as our recursion (positive bit transitions considered) 


AVMI (n ‘ +1) = L M AVMI^ n ‘) + n T (n ) 


- 2j 


fSI 

Vn 0 b 


(57) 


where /E^/W q R/B a is the signal level out of the A/D converter and, therefore, 
the input signal level to the imdphase recursive one-pole filter, and we have 
let n-^ n ^ - n^ n ^+n^ n ^ for convenience. 

From (31) and the fact that R s /B=4, we have the alternate form 


AVMI 


^ n ' +1 ) = L M AVMI^ n ') + 

M 


V ‘ «©£ 


vfT a 


(58) 


Further equating the discrete phase error j/N to the continuous equivalent 
phase error <j)/2ir, we can express (58) in the form 


x 


(n 1 ) 


AVMI 


(n ' +1) =L M AVMI (nl) + n T (n,) - 

M T it * Nq 


(59) 


Since the rightmost two terms are of zero-mean value, we see from (A-5) that 


E 


AVMI 


(n '+1 



0 


(60) 


Now consider the second moment We have 


* ~ 

Since we are dealing with low bit-error rates, this assumption 
causes a negligible error in the results. 
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From (61) and (A-13), it follows that 



(62) 


Also, from (60) and (A-5), it follows that 


E 



(63) 


5.0 BIT SYNCHRONIZER DIGITAL FILTER OUTPUT STATISTICS 

WITH NOISE ONLY— NRZ DATA 

In the last section, we characterized the output of the inphase and 
midphase filters when signal and noise were present and E^/Ng > 4 dB. At lower 
values of E^/Ng, the inphase model estimate becomes poorer and one is forced to 
reconsider the problem when the signal is not present. Since the SNR in the 
5-MHz filter bandwidth is, at most, -12 dB at the highest data rate, decreasing 
the signal will have a negligible effect on the noncoherent AGC gain. We ne- 
glect this negligible gam change between no signal and the 4-dB E^/Ng value. 

Since there are no additional approximations m the midphase case, 
we may let E^/Ng equal zero m (62) and (63) to obtain for the NRZ noise-only 
statistics 


and 


Var AVMI^ 


2 1 


1 -L, 


c 

T 


AVMI 


(n'JT - 


(64) 


(65) 


Now consider the mphase digital filter output statistics. From (37) 
with E^/Ng = 0, we have that 
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(n) 


n<"> + n <"> 

q 


21 

( 66 ) 


since it is no longer possible to employ the approximations of (30), i.e., 
Ejn^+n q ^ = E Q n ^ n ^ + n q ^ n ^ Q » Var n^ + n q ^|] = Varjn^ +n q ^ n) (67) 


If we again model the thermal noise plus quantization noise as Gaussian, we 
then have 


E 


x 



2 



x e 


2 2 
-x /2i 


/27 2 


( 68 ) 


where e 2 is the variance of ^n^+n^j. Evaluating (68) and using (27) leads 
us to 


E 







IN a 


(69) 


as is clear from (24). Further, the variance is given by (28) with y = E[x n ]=0 
and 2 defined by (27) so that 


Var n^ n ^+n^ n ^ = (l - ~^(z. 1 N a 2 ) 


It then follows that 


and 


AVIN 


(njl 


VI 


(70) 


(71) 


k\ 


Var AVIN 


(nJT _ 


(‘-f) 


2 1 N o 


1 - L-, 


(72) 
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6.0 BIT SYNCHRONIZER DIGITAL FILTER OUTPUT STATISTICS 
WITH SIGNAL PRESENT— BIPHASE DATA 

The biphase inphase digital recursive filter input when no timing 
error exists is identical to the NRZ data case. When a timing error exists, 
the output differs somewhat. We have that 


AVIN^ n+1 ^ = LjAVIN^ + | FSTB2 - LSTB2 1 


(73a) 


Assuming that E^/Nq > 4 dB allows us to write (when | jJ/N < 1/4} 

,(n) 


AVIN^ n+1 ) = LjAVIN^ 11 ^ + n^ + n ^ + X^ (73b) 


where 


X (n) = 2v¥, 




with probability Q - (74) 


and 


,(n) 


= a - 2 ^ j with probability p t (75) 


Again approximating the discrete variable j/N by the continuous-phase error <f> 
O'.e., 1 j J/N = |4>|/270, we have 


and 


x<"» 

x (n) 




with probability (l - p^) 

(76) 

with probability p t 


Averaging 



of (73) , we obtain 


E 




(77) 
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(78) 


Evaluating E ^ , 


we obtain 


E 



2 1 N a 2 + 4N ™ a Z 
N 0 


4/2" 
irVir °4 


4a 


ir/ + (h ?°d 



(79) 


It follows from (77) and (79) that 


Var 


= 2.1 N cr 2 + 4N J*- a 2 

J 


1JK p 

IT ITT 6 r t 3 


(3ir+8) 2 2 2 2 2 

3”^" a <j> P t ~ G <f< P t 

IT TT 


(80) 


Therefore, 


Var 



2.1Na 2 + 4N rj— o 2 
N 0 


4 !2 

II T 7T °g P t 


(4ir-8) 2 

^ t ' _3 G <f> 


+ - 


1 - L, 


(3iH~8) 2 2 

3 G <j> P t 

TT 


2 

3 

IT 


"A 


2 

t 


^ (81) 


Now we consider AVMI^, the midphase filter output In the biphase 
case, the integral over the transition is for a duration of one-half a bit time, 
not a whole bit time, as in the NRZ case. 

The midphase filter is specified by 


AVMI 


(n+1) .. 


L m *AVMI 

M 


(n) 


x (n) 

t “* \ 

+ MITB2*TRIND 


(82) 


Now the signal-pl us-noise component, again neglecting bit errors, is given by 


x 


(n) 





J-1 ,-l ,2,-2 , . . 



By (31), we have 



M8210-1 


ORIGINAL PAGE SS 
OF POOR QUALITY 


24 


R R s , 
B ~ NB “ N* 


(84) 


so that (83) becomes 


x 


(n) 





(85) 


if we again approximate the discrete timing error j by an approximately equiva- 
lent zero-mean Gaussian phase error <j>. We use (41) in (85) to obtain 


x 


(n) 



+ n (n) - 2 i M 
q ir 



( 86 ) 


Again, since <j> is assumed to be zero mean, we obtain 


E 



0 


(87) 


so that 


E 



(88) 


To obtain the second moment and, therefore, the variance, consider (n^ = n^+ 





a 



+ 4 N 

IT 


-b 2 
r- a 


(89) 


The second term is zero if n-j-^ and <(> are uncorrelated Since <j> depends on 
at least 20 or more bits, the correlation is very small; hence, we assume zero 
for the second term. 
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Var[^(«jj = 1.05 N cj 2 + ^-^0 2 


7F U 


(90) 


From (A-13), we obtain 


k\ 


Var AVMI 


(nil „ 


9 4Na, E, 
1.05 Nff + 2 - sj - a 


1 - L, 


(91) 


Thus, (88) and (91) characterize the discrete one-pole filter output when the 
distribution is assumed to be normal . 


7.0 BIT SYNCHRONIZER DIGITAL FILTER OUTPUT STATISTICS 

WITH NOISE ONLY— BIPHASE DATA 

In this section, we determine the one-pole digital filter output 
statistics for the noise-only case with biphase data. Since the inphase 
j^AVIN ( n filter has the same statistics for noise only with both NRZ and bi- 
phase data, we therefore obtain from (71) and (72) 


and 


AVIN 


(nJT 


. 


2.1 N 


1 -L, 


Var 


AVIN (n ^ 


i 1 



2.1 N a 2 



(92) 


(93) 


Now the midphase filter statistics with noise can be obtained from (88) and 
(91) with E b /N Q set equal to zero: 


and 


: A 


E AVMI 



(94) 



1.05 No 


1-L. 


( 95 ) 



M8210-1 


26 


ORIGINAL PA@P: ® 

OF POOR QUALITY 

8.0 PROBABILITIES OF THE LOCK DETECTOR TEST 

From (29), we see that the current lock detector test is based on 
comparing the inphase filter output with 6.5 times the magnitude of the mid- 
phase filter output. If the inphase filter output is larger, lock is then 
declared; otherwise, an out-of-lock condition is indicated. We consider the 
slightly more general test 


AVIN > r | AVMI j 


(96) 


where r is a real -valued number greater than zero. Let p^ denote the probabil- 
ity that the inequality of (64) is satisfied so that 

p n = p(AVIN^ > r AVMI^) (97) 

Denote the probability density* of AVMI^ and AVIN^, respectively, by 
P n ^AVMI^ n ^) and Pj|AVIN( n )j; then we obtain 

p(|avmi^|) J p 

r | AVMI ^ 

Let u = AVMI and y = AVIN, then, from (54), (55), (62) and (63) and 
using the Gaussian assumption, we have 


(avin^) 


d AVMI 


(n) 


d AVIN 


(n) 


(98) 



and 


P(u) 




p(y) 


l 

^*IN 



Vu (99) 


Vy (100) 


*We assume that AVMI^ n ) and AVIN( n ) have continuous-density functions 
by virtue of the assumption that the filter outputs are Gaussian. 
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where VarjAVIN^J = a IN 2 , Var[AVMl( n )] =a HI 2 . E[AVIN( n )] =y JN and E 
If we let x - [u] , we then have 


p(x) = 


2 ; , 2 

2 ~ X / 2ct MI 


JS cr, 


MI 


x > 0 


so that x=|AVMl|. Using (96) and (97) in (94) yields 


/' 


n J /27 cr 
0 ^ ff MI 


2/o 2 

_2 e ' x / 2ct MI 


r_i_ 

rx ^ °IN 


-(■‘'-■'in)' 


2a 


IN 


A change of variables leads to 


Pn - r-i-. 

" / ^«HI 


-Xi 2 


ERFC 


rx - ]i 


IN 


/2 


IN 


dx 


where 


ERFC(t) 



dt 


First consider the NRZ case. The original parameters for 
digital filter were 



252 

256 


0.9844 


_ 242 
M " 256 


0.9453 


AVIN I 6.5 | AVMI | 


and the new recommended values, to be discussed below, are 


L i = Mf = 0 9922 ; L M = S = 0 9688 ; AVIN > 13 l AVMI l 


27 

[avmi ( n )J=o 
( 101 ) 


dydx (102) 

(103) 

(104) 

the one-pol 

(105) 

(106) 


The most desirable fix to the lock detector would utilize new narrower band 
filters for the lock detector which would be distinct from the loop inphase 



M8210-1 


28 


and rmdphase filters; however, the modified version (corrections, etc*) of the 
microprocessor does not have excess memory or the speed necessary to implement 
the additional lock detector filters. Consequently, single filters will be a 
compromise between the lock detector and the bit synchronizer loop requirements. 
It should be noted that Richard Helgeson has tried these NRZ parameters values 
(from (106)3 and found them to be quite satisfactory. The Manchester mode of 
the bit synchronizer had original parameters of 

«r 

4 = Mf = 0-9844 ; [_ M = ||| = 0.9844 ; AVMI > 5 | AVMI { (107) 

The new recommended values are 

4 = IIf = °- 9922 1 4 = Sb 0 - 9922 • AVIN > 15 |AVMI| (108) 

Note that Lj applies to both the AVIN and AVMI filters. 

As mentioned previously, these filter L values were compromises since 
they were constrained to be of the form N/256, where N is an even positive inte- 
ger. Various values of r were tried for the new parameter values. Due to the 
approximations in the analysis and the computer program involved, it was deter- 
mined that, for given "L" values, p^ would saturate at around four Q.9999z, 
where z is the fifth digit and is between 1 and 9. Consequently, by narrowing 
the digital filters and increasing r, it was possible to keep the probability p^ 
high (0.9999z) when the signal was present and simultaneously reduce the proba- 
bility p n when only noise was present. This will reduce the time to indicate 
lost lock due to signal dropout. In other words, the detection probability could 
be maintained and the false alarm reduced, thereby reducing the mean time to 
indicate the out-of-lock condition when the signal drops out 

The results for the NRZ and biphase cases are indicated in Tables 2 
and 3, with P^ denoting the value of p n when the signal is present and Pp^ is 
the value of p n when the signal is absent for all data rates. Equation (103) 
was used to evaluate these probabilities. 
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Table 2. 


NRZ Data 


r- < dB ) 

Performance 
wi th 
Original 
Parameters 

Performance 

with 

Recommended 

Parameters 

f 4 

0.999921 

0.999940 

p rf 6 

0.999921 

0.999940 

d \ 7 

0.999921 

0.999940 

Pp^j Noise Only 

0 98850 

0.9478 


Table 3. Biphase Data 


Jr fdB > 

Performance 

with 

Original 

Parameters 

Performance 

with 

Recommended 

Parameters 

( 4 

0.999940 

0.998923 

P H \ 6 

0 999940 

0.999921 

d \ 7 

0 999940 

0.999954 

Pp^| Noise Only 

0.98783 

0.769926 


In order to determine the mean time to indicate the out-of-lock state 
when the signal is not present, we must evaluate 



1 ' P FA + 2 


( 1 - p fa) p fa + 4- p fa) p fa 2+ ••] < 109) 


where Pp^ is the false-alarm probability and T M is the "memory time" of the 
inphase digital filter output*. 

From (108), 


t ol - ( 1 - p fa) t m[ 1 + 2p fa + 3P fa 2 + 


( 110 ) 


* 


The memory is approximately T^/Q ~ 





M8210-1 
and, since 


QR1GIWAL PASS IS 

OF POOR QUALITY 

1 + 2ct + 3a 2 + 4 a 3 + • ♦ • = — 

(l-«) 


we have the simple result: 


'M 
1 - P 


FA 


30 


(111) 


( 112 ) 


Evaluating T^ for the four cases of Tables 2 and 3 yields the results presented 
in Table 4 for both the NRZ and biphase cases. 


Table 4 Mean Time to Indicate Out of Lock (Noise Only) 


Case 

T ql (Seconds) 


16 kbps 

8 kbps 

4 kbps 

2 kbps 

1 kbps 

NRZ - Original Parameters 

0.348 

0.696 

1 392 

2.784 

5 568 

NRZ - Recommended Parameters 

0.154 

0.307 

0.614 

1 228 

2 456 

Biphase - Original Parameters 

0 329 

0.657 

1.314 

2.628 

5.256 

Biphase - Recommended Parameters 

0.0348 

0.070 

0 139 

0.278 

0.556 


Thus, we see that the recommended parameters for the filters and lock detector 
reduce the time to indicate the out-of-lock condition to one-half of the origi 
nal case for the NRZ mode and about one-tenth for the Manchester mode 

It should be emphasized that the loop performance could have been 
considerably improved if separate filters were utilized for the lock detector 
since they would have been made narrower Due to the changes and additions in 
the microprogram for saturation, etc., however, the additional room for the 
extra filters was not available. Furthermore, the constants relating AVIN and 
| AVMI | were constrained to be of the form 2 n x6.5, where n=l, -1, 2, -2, etc 
for the NRZ mode and, for the biphase mode, were constrained to be 5, 10, or 
15. Consequently, only approximately optimal values could be used It should 
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also be emphasized that, because the DC bias-removal and scaling algorithms 
were unavailable, their effects could not be included in the analysis pre- 
sented here. It is believed they will not affect the analysis during tracking; 
however, when the signal is removed, the DC bias-removal algorithm could con- 
ceivably cause the lock detector to indicate the out-of-lock state more quickly 
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APPENDIX A 


EVALUATION OF THE MEAN AND VARIANCE OF 
THE RECURSIVE ONE-POLE DIGITAL -FILTER OUTPUT 


In this appendix, we derive the steady-state mean and variance of the 
recursive one-pole digital filter used for both the inphase and midphase filters. 
Let be the value of the filter output at time n. We then have 


-(n+1) = L z (n) + x (n) 


0 < L < 1 (A-l) 


where x^ 'is the nth bit sample of the input digital matched-filter output. 
In (36), we showed that its SNR is 1.9 E b /N Q . Solving (A-l), we obtain 


: (n+1) = I x. L n_i 


(A— 2) 


To obtain the mean E z ^jL we average both sides of (A-l) to achieve 

E |z n+1 J = Lj e|z^ + E x <"> 


(A-3) 


In steady state, 


E [>|] = E [z<" + 1] 


(A-4) 


so that (A-3) yields 


;-(nr = 


:[x ( ^ 

1 - L 


(A-5) 


Now we consider the variance of z^ nj . First, we square and take the ensemble 
average of (A-l) to yield 


,(n+l) 


= L* E (z<">) 2 + 2 L E p>z (n >] + e[[^ n >) Z T (A-6) 
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OF POOR QUALITY 


EUz <n+i) ) i - E rf z ‘">) 2 


so that 


:p>] Ep] + e[(x<">) 2 ~ 


1 - L 


By definition. 


Therefore, 


Var z 


= E(z <n) ) ? - E 2 [z (n ^] 


2L E[x( n )jE[z( n ) + E (x(h)) 2 ] E 2 [x^ n ‘ 


1 - L 


(1 - L)‘ 


■M 


2(1-L) L E 2 jx( n )j + (1-L) 2 E (xtn))^ - (l-L 2 ) E 2 [x (n ^J 


ci-oa -iy 


After simplifying, we obtain 


z ^ n |J = 


: (x(n))^J - E 2 [x ( n 


or, finally, 


i- -> Varfx^"] 
. z (n) = LJ 

L -1 1 - l d 


Equations (A-5) and (A-13) constitute the main results of this appendix 
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APPENDIX B 

ORIGINAL PSP BIT SYNCHRONIZER ALGORITHM 


The following documentation is a slightly updated form of an abbre- 
viated flowchart of the PSP bit synchronizer algorithms for NRZ and biphase-L 
data. They are not considered to be exact and can be compared to the "final 
version" given in Appendix C. The NRZ algorithm includes a transition check 
which was not in the original algorithm, plus the changes of the lock detector. 
The biphase-L algorithm has not been updated due to the lack of a complete al- 
gorithm; however, the lock detector portion does have the new parameters. This 
flowchart was developed at Johnson Space Center by Brett Parrish independently 
of any flowcharts that TRW might have used in the development of their bit syn- 
chronizer program. The PSP bit synchronizer "machine-code" program (as supplied 
by TRW) was disassembled and flowcharts were generated from it. 

B 1 NOMENCLATURE 

Each bit of incoming data is sampled at a rate of 1.024 MHz and inte- 
grated (summed) in hardware. At the end of the bit period, the hardware inter- 
rupts the bit synchronizer computer (which is caught in a loop at program loca- 
tion PEXEC) and the "bit-time interrupt" routine is performed. 

The bit-time interrupt routine polls the hardware for three partial 
integrals, as shown in Figure B.l (assume that the bit synchronizer is locked). 
When the bit synchronizer is locked, the hardware takes 1 .024 x 10^/DATA RATE = 

Ns samples per bit. 

The current microcode program utilizes DMPNM, as described below. 

The new program was modified to avoid saturation but, due to time limitations, 
has not been introduced in the discussion to follow. Whenever a phase correc- 
tion is required, the software changes the nominal dump count by a small amount 
±e. 


DMPCU = DMPNM + DMPCP + DMPCD 

(nominal) (a$) (rate tracking) 


where 

DMPCU = current value of dump count to be output to hardware (number 
of samples/4 to be taken during the next bit period) 
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NRZ Input 


+ 

0 


f l 

H — Bit Period — 

! MITB2 1 



FSTB2 LSTB2 


NOTES FSTB2 
MITB2 
LSTB2 
PRTB2 


first half-bit integral 

middle half-bit integral (used in biphase mode only) 

last half-bit integral 

last half-bit integral of previous bit 



PRTB2 FSTB2 LSTB2 


NOTES For NRZ data formats, the software implements two filters: 


] FSTB2 + LSTB2 1 ->- inphase filter 
(PRTB2 + FSTB2)*TRIND midphase filter 


TRIND is the transition indicator and has a value of 
+1 or -1 or 0 The MITB2 value is not used in the 
NRZ mode 


Figure B 1 Half-Bit Integrals Used in the Algorithm 
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DMPNM = nominal dump count = Ns/4 

DMPCP = dump count correction due to phase-tracking loop 

DMPCD = dump count correction due to data-rate tracking loop 

Bit-time interrupt routine (executed after hardware integrates over 
one bit period): 

1. Outputs previous bit decision to frame synchronizer circuit 

2. Takes current LSTB2 and puts it in PRTB2 

3. Inputs FSTB2 from hardware 

4. Inputs MITB2 from hardware 

5. Inputs LSTB2 from hardware 

6. Outputs new DMPCl) to hardware 

7. Sets DMPCU equal to DMPNM 

8. Set NDATR flag 

9. Returns to PEXEC loop in main routine. 



Loop here while 
hardware is .X 
sampling and. 
computing 
integral s 


Before bit-time' 
interrupt occurs, 
the computer is 
looping here 


( PEXEC ) 


B.2. NRZ Algorithm 


B4 
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•NDATR flag is set during bit-time 
interrupt routine (after hardware 
has integrated one bit, it inter- 
rupts the computer) 


For 16 kbps data rate, 
skip DC offset routine 


Check frame 
synchronizer 
and data rate 
correction 



DC offset 
routine 


~ ^PC0FF ) 


offset routine 

=□ 


NOTES- 


-label only 



Frame synchronization 


Check Hamming 
Computation 




J 


2. Executes once 
each bit period 


No doppler 
correction 
in acquisition 


, Counts to 4096 when CUDRC = 0 

[CCID - 8 -> CCID 
^-XTime for correction 
\ (if underflow 
Yes 



Signal format* ( SGFCK \ 
.check f v — t — — y 



Every 4096 bits set 
RTIFL flag 
(New rate-tracking 
interval required) 


CUDRC is 
set to 0 
or 1 or -1 
during 

NRTIN routine 


( NR ZIP ) 

NRZ branch 


( BIPIP ) 

' Biphase-L branch 
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Inphase averaging filter / — * — :* T*. rTO o i 

3 3 \ ( NRZIP J I I FSTB2 and LSTB2 are integral: 

V t ^ t t computed by the hardware and 

L I _ !‘!n j Uri . ng acquistionr— 1 ■ — — - lupdated every bit time 

L T = 0 99 during steady AVIN - Lj*AVIN+ FSTB2 + LSTB 2 1 K 

state — * * " i ^ 


I I FSTB2 and LSTB2 are integrals 

f f computed by the hardware and 


Data ] 

Rate 

(kbps); 


3-d 

B point (Hz) 


Aqq. 

Gain 

ss 

(TRK) 

Max. 

Gain 

8.16 

20 

1.6 

100 

16 33 

20 

3.2 

100 

32 65 

20 

6.4 

100 

65.31 

20 

12.8 

100 

130.6 

20 

25.6 

100 


Set Bit 
Decision 
- 1 


Set Bit 
Decision 
= 0 


Bit decision 
also stores 
bit for TRIND 
cal culation 


/TranX. 

sition^ 
5 os or Neg. 


Midphase 

Averaging Filter: 

S M = 10 

L„ = 0.90 during acq 
M 

L„ = 0.95 during loci 
M 

Gain = 10 during acq 
= 20 during lock 


■TRIND 


/iran-N. 

' sition\!' 
since last/” 

X. bit / 

\Neg 

k==i — 

-1 TRIND 


Bypass if no 
^ transition 


► |AVMI * L M *AVMI + (FSTB2+PRTB2) * TRIND 

TRK AC Q 

♦ /%? \ No , 

x “ A i mp H pI Tu * 


Filtered update 
signal triggered 
by transitions onl; 

Lock flag test 


Steady-State 
Bit Synchronizer Light 


Set/n del ay r use k ^ acquisition 
.1 oc k flag 4096 \ 

Yv qc > transitions ^Ov 


Test to see “are we 
still in lock 7 " 


Lost lock 


LSTLK 


Clear Lock Flag 
Clear- Hardware Lock Ind 
Set to Acquisition Mode 


Phase-tracking 

1 -> PTC IN 1 correction interval 
:q+ mftCv Vn TTnTTTI ™ acquisition _ 


Stea djn 
State/ 

N4 


V IN >13 *] AVMI 

NSSSs Ty?$ // ^ 


.XAVIN 

(0^10 * j AVMI | forces 
>.40 transitions^ 


in a row. 

\ 


Set MFTCV as follows 
Rate (kbps) MPTCV 

L6 
8 
4 

8 2 

16 1 


Clear AVIN, AVMI 


PTRAK 


Avoids hang-up in v 
lock mode when out 
of lock * 

C PTRKl) 


Set Lock Flag 
Clear AC0R 
4096 -► CNTR 


/1024V 
nontran-\ 
itions ny 
\a rrn/ 


See Next Page 











/ Enter once 
\every 4096 bits 


Switch to 
steady- state 
parameters — ^ 

l~ 

Set loop to SS 
Set hardware lock ind 


'>} ( NRTIN ) 

r ' 

Clear RTIFL Flag 

No v^Stea dy- \ Yes 
\ Sta te 


Rate tracking 
inhibited during 
acquisition 


Bit synchronizer lock 
indicator to 0/1 MDM 


Set to lock mode: 

0.99 for in phase filter 
0,95 for midphase filter 


original page Ig 

OF POOR QUALITY 


Set phase-tracking correction 

interval as follows: 

kpbs 

Interval (see NTSPC) 

1 

1 

2 

1 

4 

2 

8 

4 

16 

8 


[Set phase- tracking correction value 
MPTCV to 1 (if 1 kbps, set it to 2) 


= 5 * AVIN 


= 2 * AVMI 


Filter "gain has changed, 
so normalize current value 

Makes acquisition parameter 
"gains" the same as tracking 
J'gains" (avoids transient) 


( NRTT1 
| ACOR | + / 


( Magnitude of accumulated 
\ dump count over 4096 bits 


1 CUDRC 


Pos/ AC0R \Nec 
— - ^is pos 

X.neg/^ 


-1 CURDC 


Don't correct 
rate 'v. 

r-^ 

0 CUDRC 
32768 CCID 


Clear ACOR 


ABOR > NSPB8 


NSPB8 


32768 

TOT 


+ CCID 


1 

128 

2 

64 

4 

32 

8 

16 

16 

8 


pcypr ’N-a-Go to PEXEC and 
..J wait for interrupt 















( PTRK1 ) 


i 
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( PTRAK ) I ~- 

v — — y — ' ^During Acquisition, NTSPC = 1 

/ 


B7 

Phase Updating 



During SS 

NTSPC 
(# Tran- 
sitions) 

(kbps) 

1 

1 

1 

2 

2 

4 

4 

8 

8 

16 


MPTCV - Maximum 

Phase-Tracking Correction Value 


i ^ /Update dr retard 


ACQ 

kbps 

MPTCV 

1 

16 

2 

8 

4 

4 

8 

2 

16 

1 


SS TRK 

kbps 

MPTCV 

1 

2 

2 

1 

4 

1 

8 

1 

16 

1 


to data rate 


ACOR = ACOR + DMPCP + DMPCD 



ACOR is the accumulated dump count 
correction over 4096 bits 


•Hardware dump count is made up 
of nominal dump count plus 
"dump count correction" (phase) 
plus dump count correction 
(data rate) 


r Go back to* 
PEXEC and wait 
for interrupt 


New rate-tracking 
jlp'jX^nterval required flag 

Flag 
Set 4 

es 

C mrtin — ' Do once every 4096 bits 


(To Next Page) 









ORIGINAL PAGE Eg 
OF POOR QUALITY 


Bit Time Interrupt Routine (Executed after hardware integrates one bit) 











Loop here while 
hardware i s ^ 
sampling and 
computing 
integral s 


Before bit-time' 
interrupt occurs, 
the computer is 
looping here 


PEXEC 


B.3. Biphase-L Algorithm 

(Executes once each bit period) 


f - NDATR flag is set during bit-time interrupt 

routine after haredware has integrated one bit, 
it interrupts the computer. 


16 kbps) 


For 16 kbps data' 
rate, skip DC 
offset routine 


DC offset 
routine^ 


'4 or^ 
J3 kbps 


DCOFF 


DC offset routine 


Integration by parts 
for 1 or 2 kbps 


Check frame 
synchronizer 
and data rate 
correction 


FSYNC 


Check Hamming 
Computation 


^Counts to 4096 when CUDRC = 0 


p-No doppler 
correction in 
acquisition 


'Time for 
Data Rate 
Correction, 
\ ? / 


CCID - 8 + CCIP 
Time for correction 
if underflow 
Yes 


Put contents of 
CUDRC in DMPCD 


Signal 

format check 


SGFCK 


Every 4096 bits 
set RTIFL flag 
(New rate-tracking 
interval required) 


CUDRC is 
set to 0 
or 1 or -1 
during 

NRTINE routine 


NRZ / Data \ Biphase 
\ For 


NRZIP 


NRZ branch 


BIPIP 


Biphase-L branch 






Biphase Mode 

Inphase Filter* 

Lj = 0.99 SS or Acq. 


i. 

( BIPIP) 


TB2 


AVIN = L*AVIN + 

L ' 

FSTB2- LSTB2 

1 


r \ 


BIO 


FSTB2 and LSTB2 are 
integrals computed by 
the hardware 


Data 

Rate 

(kbps) 

: 3 - dB Point 
(Hz) 

(Hz) Gam 

1 

1.5 

100 

2 

3.2 

100 

4 

6.4 

100 

8 

12 8 

100 

16 

25 6 

100 


Pos / FSTB2" 


Neg 



'^\LST^ 

1 Bit Decision = 1 


Bit Decision = 0 

1 - 1 - 

TRIND 


+1 TRIND 


— — **- — 


— * — i 


TR IND =#* ± bit, or 
no bit transition 



AVMI = L * AVMI + MITB2 * TRIND 


Midphase Filter 

(L = 0.95, Gain = 20) - Acq 
(L = 0 99, Gain = 100) - SS 


Test to see 
are we still 
in lock? 



( LSTLK) 


Yes 


(L j = 0.99, Gain - 100_ 


A1 ternate 
inphase filter 


ALIN = L*ALIN + [ PRTB2 - FSTB2 


Bame as for NRZ Mode 


Delay Hardware] 
Integrator 
1/2 Bit 


Switch AV IN 
and ALIN 


pet Lock Flag 
Clear ACOR 
4095 CNTR 


I 

(ptrAK ) 
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c.-v 











PTRAK 
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NRZ 

.During Acquisition, NTSPC 


During SS 
NTSPC kbps 


Have NTSPC X 
transitions occurred 
.since last phase > 
\ correction 


NTSPC 

kbps 

1 

1 

1 

2 

2 

4 

4 

8 

8 

16 


>65.0*|AVMl}> J^<AVIN > 156* | AVMl] 


Biphase- 


Biphase 

Lock SS 

NTSPC 

kbps 

4 

1 

8 

2 

16 

4 

32 

8 

64 

16 


PHCOR 


Pos >^AVMrv Neg / 

[ - <^os or Ne<^> — ] p 

-MPTCV + DMPCP \ / MPTCV -> DMPCP 


PTRK1 


ACK ^ 

SS 

kbps 

MPTCV 

kbps 

MPTCV 

1 

16 


1 

2 

(1) 

2 

8 

(1) 

2 

1 

(1) 

4 

4 

(1) 

4 

1 

(1) 

8 

2 

(1) 

8 

1 

(1) 

16 

1 

(1) 

6 

1 

(1) 


PEXEC 


ACOR = ACOR + 

1MPCP + DMPCD 



DMPCU = DMPCU 

i- DMPCD + DMPCP 



Clear DMPCD, DMPCP 


fio back to PEXEC 
and wait for 
interrupt No 


(ACOR is the accumulated dump count 
I correction over 4096 bits 

/Hardware dump count is made up of 
Jnominal dump count plus "dump count 
\correction" (phase) plus dump count 
/correction (data rate) 


MPTCV Maximum phase-tracking 
correction value 


KTIFL 

Flag 


New rate-tracking 
interval required flag 


NRTIN 


■Do once every 4096 bits 


(To Next Page) 












PEXEC 


-Go to PEXEC and 
wait for interrupt 
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Bit Time Interrupt Routine (Executed after hardware integrates one bit) 
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APPENDIX C 

NEW PSP BIT SYNCHRONIZER ALGORITHM 


This new algorithm was made available to the author by Brett Parrish. 
The changes in the original algorithm leading to this algorithm were made pri- 
marily by Rich Helgeson and Brett Parrish. This program is current as of 
October 22, 1982. 
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i . C3 



Set up the 
following constants 


00001111 DMPNM 


Will be changed to 14 
( 00001110 ) 


[Rate j 

NSPB8 MPTLA A'BTIPB BTIPQ 


10000000 oO'oiOOoo 00010000' 00001000 


01000000 00001000 00001000 00000100 

61 

00100000 00000100 00000100 00000010 

99 ^ 

00010000 00000010 00000010 00000001 

mmmm 

00001000 00000001 00000001 00000000 


TEN*KBPS -»■ RSMPY 


t 

Clear BTICT 
Clear BTIQT 
Clear QBICT 


Set ud’ 

r these constants* 

Mask Number of Frame Sync Words From FSYLC 


’ 

No of Frame Sync Words 

FSYNL FSYNP * (8*FSYNL) 

1 

2 

3 

4 

oooooooi 00001000 
00000010 00010000 
00000011 00011000 
00000100 00100000 
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C4 


Set up these constants. 1 


^timber of Frame Sync Bit Errors 

, Number of Frame Sync Bit Errors 
for Inverted Data 


L 


No. of Frame Sync Words 

- 

NALDT^ 

^ NALDI 

1 

0 

'8 ! 

2 

1 

15 

3 

1 

23 

4 

2 

30 


Point HMGWA at Hamming Word 1, 2, 3 

or 4 according to the number of 1 

Frame Sync Words 

No of Frame Sync Words 

HMGWA Points to 

1 

HMGW1 

2 

HMGW2 

3 

HMGW3 

4 

HMGW4 



SYNC 1 -> PSYN4 


SYNC 1 4 - PSYN3 
SYNC2 -> PSYN4 


SYNC1 


PSYN2 

SYNC2 


PSYN3 

SYNC 3 

4 

PSYN4 


SYNC1 4 

PSYN1 

SYNC2 4 

PSYN2 

SYNC3 4 

PSYN3 

SYNC4 4 

PSYN4 


( INIP4 ) 


Out 

DMPNM 

■f 

Hardware 

13 


DMPNM 

-»• 

DMPCU 


Out 

NMLEN 

to 

Hardware 

12 

Out 

NULEN 

to 

Hardware 

12 
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INIPM (Cont't) 


C5 


Biphase —Sro rma 


121 

126 

126 

126 

5 

1 

KB PS/2 or 
1 if KBPS = 1 
4*KBPS 


AMIDC - 
SMIDC - 
AINPD - 
SINPD - 
ACDRA - 
SSDRA - 

PTCIA - 

PTCIS - 

MPTCA - 


1 LPDRA 

9 CUDRP 


Midphase Average Filter Decay Constant (ACQ) 115 

Midphase Average Filter Decay Constant { SS) 121 

Inphase Average Filter Decay Constant (ACQ) 121 

Inphase Average Filter Decay Constant (SS) 127 

Filter Decay Constant Ratio (ACQ) 2 

Filter Decay Constant Ratio (SS) 5 

Phase-Tracking Correction Interval (ACQ) 1 

Phase-Tracking Correction Interval (SS) j l^if^KBPS°=l 

Max. Phase-Tracking Correction Value (ACq)j^j t ^KB PS°= G l° r 

ACQ. to SS Normalization Multiplier 5 

Current Filter Decay Constant Ratio 5 



1 -> PIFLG (Enable Interrupts) 
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( CFGTO ) 

r 

75 ys Delay 

: ; * 

CFGWA = Ho. of CFGWD 

( pixie ) 

— I 

1 PIFLG 
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Point CFGWA at CFGWD 


■Enable Interrupts 


. SCLIN ) 

X t _ 

Clear NDATR 


/KBPs X Yes 

V 16/ 

TNo 

PRTB2 = LSTB2 


►RSMPY = 


FSTB2 

LSTB2 


FfIBIU, FHBIL *2* RSMPY 
256 

LHBIU, LHBIL *2* RSMPY 
256 


PRTB2 

= LSTB2 

,r 

LSTB2 = 

CULHI 


FSTB2 = FSTB2 + OFFCO 


SGFMT 

\<4y 


Maximum value 

of FSTB2 or Cl 

LSTB2 or MITB2 |f STB 2 = 

at this point r— 

is ±120 *“ 


< Dver\ pos 1 

Flow/ 

X/ (fSTB2 = +127 
Ino _ n 


LSTB2 = LSTB2 + OFFCO 
Neq /h^r\ Pqs 


Biphase 


LSTB2 = -128 


LSTB2 = +127, 


NRZ Format 
Yes 


SCC03 


| MITB2 _ MHBIU, MHBIL *2* RSMPY 
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C7 
















( SGFCK ) 


NRZ 


/ l=T 

C liRZip ) 

x — 



Biphase 


AVINUsL x 20NPDC + 


256 


FSTB2+LSTB2 I 


c 


AVINU ,L 

zr 


NRZTR 

X 


3 


During acquisition 
INPDC = 121 0 

la =a i _ 1 (0 95)+| FSTB2+LSTB2] 

| During steady state 
INPDC = 126.0 
,a 


7 = 1 — >, 

r ipmp ) 




00 

ro 


INPMU ,L = FSTB2 - LSTB2 

( BIPTR ) 

.. . * 


Put current bit decision into 
previous bit decision 


i 


a i _ 1 (0.99)+|FSTB2+LSTB2 


CUDEC PRDEC 


X 

Clear transition indicator 



Set NLBIT=0 


Set current bit decision 


Set current bit decisior 

>- 

A 

X 

f 
1 — 1 

t — f 
1! 


= 0 


Pos 

/\ FSTB2\ 
/- LSTB2 \ 
N^ositi ve / 

Neg 

1 CUDEC 

-1 -*■ CUDEC 

-1 -1- TRIND 

N. / 

1 -> TRIND 

1 NLBIT 


1 ■> NLBIT 

1 


i 


* 1 



f BDCVT ) 

Check 



^■Assembly 
Yes Code 


Counts 

Transitionless 

Bits 



( BIPMP ) 


O O 
"n Xj 

"D 

O “5 

O rfi i 

AJ |f*« 

rO "© 
C > 
> Q 
r; m 

9$ 


O 

CO 


( To NRZMP ) 


( RTRAK ) 


C nrfil ) 



















CO 

PO 


o 

r 


o 

LO 


OftiGlMAL PAGE 

OF POOR QUALITY 
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Cll 


Note 


If lock is lost, 
set parameters to 
acquisition mode. 
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PTRAK 


C12 


NTSPC = 

NTSPC - 1 

^0 /fTTCn 

r - n\ 

"C n l o r 

L» ” U ^ 

^Yes 

NTSPC 

= PTC IN 
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Pos / AVMI\ Neg 


-MPTCV+DMPCP ■* 

— 

"^Positi ve^>—- — 

-►MPTCV-* DMPCP 


( RTRAK ) 

T 

CITDUjL = CITDLIjL + AB0R 


'''Overflow'''' 
on CITDU ,L 


RTRAK1 


CUDRC -»• DMPCD 


4096 -v CNTR 
-1 -»■ RTIFL 


■"''Done > 
4096 
.Times 


PTRK1 



ACORU ,L = ACORU 

.L+DMPCP+DMPCD 

\ 


DMPCU = DMPCU + 

DMPCP + DMPCD 


Clear DMPCD 
Clear DMPCP 


/'RTIFN S 
Flag Set 
\ ? ^ 


NRTIN 
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INP1U,INP1L 

+• INP2UJNP2L 



INPHUjINPHL 

* INP1U, INP1 L 


(FHBIU,FHBIL + LHBIU,LHBIL) -»■ INPHU,INPHL 


No transition 


TRIND = 0 

\ , / 


DCMNU, DCMNM, DCMNL 
+ INP2U,INP2L 
+ INP1U,INP1 L 
DCMNU , DCMNM , DCMN L 

DCOBC = DCOBC + 1 


DCOBC = 64 

\ » / 


OFFCO 


^DCMNU, DCMNM, DCMNLj ^ p <; MpY 
256 


Clear DCOBC 

I — 

Clear DCMNU, DCMNM, DCMNL 
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C BTIML ) Bi t time interrupt 
i ' y 1,2, 4, 8 Kbps 


C15 


XR + SAVXB , A + SAVB1 , B -> SAVB2, CY SAVB3 


original page is 

OF POOR QUALITY 



Clear BLFLG 

i — 


NLBIT LDNRZD 

1 


Input Current Last-Half Interval from Hardware 1 CULHI ' 
Input Current First-Half Interval from Hardware ‘CUFHI * 


I 


j Out Current Dump Count^to Hardware 1 DM^CU 1 ] 


C 


DMPNM -> DMPCU 


: 


X -s- SAVB4, Y SAVB5 

i —i . = — Z 

1 Set Quarter-Bit Clock Line to 0 1 LDQBCL I K 1 



No/ Mes First 
.BTICT=g> j 


ICUFHI -> AQB1U ,L ] 


|CULHI-?-AQB2U,L 


X 


-1 -»■ BLFLG 

X 


Set Quarter-Bit Clock Line to 1 ‘LDQBCLK. 1 

X 


BTIQT = BTIQT + 1 1 


cWr 

+ CULHI 
AQB4U,L 



AQB1U,L 
+ AQB2U,L 
FHBIU,L 


■1 + BTICT 

i . — 


LDCKFF 


Clear BTICT 

t 


1 + NDATR 


Clear BTIQT 
1 -> NDATR | 
T 



Clear SHFTH 
Clear BTICT 



CUFHI 
+ CULHI 
AQB2U,L 


FHBIU ,L = AQB2U,L +AQB1U.L 
MHBIU,L = AQB2U,L + AQB3U,L 
LHBIU,L = AQB3U,L + AQB4U,L 


Yesr~ 



Clear SHFTH 
Clear BTIQT 


- C BTI15 y *- 


SAVB5 4, SAVB4 -> X, SAVB3 CY , 

SAVB2 B, SAVB1 * A, SAVXB + XR 


Q Return 
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Bit time interval for 16 kbps only 
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Echo Configuration to Hardware 


Hardware issues PROGM Interrupt 
six times to set configuration 
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APPENDIX D 

AMENDMENTS TO BIT SYNCHRONIZER ACQUISITION TIMES 


The bit synchronizer shall acquire in four seconds with a probability 
of at least 0.9, with the E^/N Q value specified in Table D.l below for the given 
data format. 


Table D.l. Data Format Versus Data Rate E^/Nq Values 


Format 

E b /N 0 ( dB ) Required 

1 kbps 

2 kbps 

4 kbps 

8 kbps 

16 kbps 

NRZ 

U M, S 

9.5 

6 5 

4.5 

2.5 

5 5 

1 

m 

■ 

9 5 

mm 

5 5 

3 5 

5 5 

















10.0 


SHUTTLE/TDRSS AND SSTDN COMPATIBILITY ANALYSIS 


While Task 5, Space Shuttl e/TDRSS and GSTDN Compatibility Analysis, 
was directed at the overall system compatibility, many of the analyses were 
performed as part of the individual communication subsystems, such as the $- 
band network equipment and the Ku-band communication system. Under this task, 
Axiomatix attended all TDRSS design reviews, evaluated review material and, 
together with NASA JSC personnel, generated Review Item Dispositions (RID * s ) . 
The final TDRSS design review which Axiomatix attended was the TDRSS Ground 
Segment Final Design Review at TRW, September 22-24, 1981. At this design re- 
view, Axiomatix was still concerned that the ground segment decoder could not 
achieve branch synchronization with the all -zero sequence. An action item was 
submitted at the design review to resolve this problem. The final resolution 
was to inform the users in the payload ICD (NASA JSC ICD No. 2-19001) that 
branch synchronization is not achieved until the number of transitions is 
greater than 64 and the maximum number of consecutive bits without a transi- 
tion is less than 64 within any sequence of 512 bits. 

During the design review, TRW stated that the Ku-band open-loop 
pointing budget would be revised. While relaxation of the open-loop pointing 
budget could be troublesome for Shuttle Ku-band acquisition, it was learned 
that TRW intends to tighten this budget. Axiomatix reviewed the current point- 
ing budget and found that no changes had been made during the past two years. 

An action item was written in order to prompt TRW to update the pointing budget 
in light of current knowledge regarding the system parameters. 

This section discusses two of the analyses performed to establish 
compatibility between the Shuttle Ku-band communication system and the TDRSS. 
These analyses involve (1) the TDRSS Antenna Scan for Shuttle Acquisition and 
(2) the Power Spectral Density of Staggered Quadri phase PN With Identical Se- 
quences. For additional analyses concerning Shuttle/TDRSS and GSTDN compati- 
bility, refer to sections of this report dealing with the individual communi- 
cation subsystems. 

10. 1 TDRSS Antenna Scan for Shuttle Acquisition 

Axiomatix was asked to review and comment on an investigation funded 
by Goddard Space Flight Center regarding the feasibility of effecting a TDRS 
scan during Shuttle acquisition. This study was conducted in response to a 



concern voiced by NASA/JSC and Axiomatix as to the ability of the TORS Ku-band 
to open-loop point to within 0.22° of the Shuttle. 

The study explores the possibility of using externally generated 
(NASA) or internally generated (TDRSS) scan patterns. In the former case, a 
NASA-supplied processor inputs updated state vectors to the TDRSS from the NASA 
Network Control Center (NCC). The disadvantage of this approach is that, in or- 
der to account for TDRSS/NASA propagation delays of the state vectors and signal 
status, the scan must be slowed down. The latter approach uses a TDRSS computer 
to generate a scan, which conceivably could result in a faster scan since the 
NASA/TDRSS interface is eliminated. Since it corrects on-track errors only, the 
use of a At adjustment is not considered an acceptable solution. 

Axiomatix finds no fault with the conclusions reached in the Goddard 
study; however, the fundamental systems question— how to implement an acquisi- 
tion procedure— was not addressed, i.e., "what signal is the TDRSS supposed to 
find?" It is unlikely that the Shuttle will be able to open-loop point to the 
TDRS accurately enough to provide a strong signal, and we cannot expect both 
the TDRS and the Shuttle to scan simultaneously. 

10.2 Power Spectral Density of Staggered Quadriphase PN 
with Identical Sequences 

Axiomatix Report No. R8010-7, "The Power Spectral Density of Stag- 
gered Quadriphase PN With Identical Sequences," dated October 27, 1980, anal- 
yzes some of the TDRSS waveforms to determine the power spectral density and 
the effect on processing gam against interference. 
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THE POWER SPECTRAL DENSITY OF STAGGERED QUADRIPHASE PN 
WITH IDENTICAL SEQUENCES 

By 

Marvin K. Simon 

1.0 INTRODUCTION 

Staggered quadriphase PN is a direct-sequence spread spectrum 
technique which employs a complex spreading waveform. Such a communica- 
tion system offers twice the processing gain of the corresponding system 
employing a real spreading waveform since, in the former, the interfer- 
ence energy will be distributed equally in the two quadrature channels 
while, in the latter, all the interference energy appears in the data 
channel . 

More often than not, the two PN sequences will be independent 
of one another and, thus, the power spectral density (PSD) of the spread 
spectrum waveform is determined by the waveform of either sequence In 
this report, we examine the effect on the PSD using the same PN sequence 
on both channels (except for the 1/2-chip delay between the two). 

2.0 CALCULATION OF THE POWER SPECTRAL DENSITY 

A staggered quadriphase PN spread spectrum signal may be written 
in the form 

x(t) = /2S Rcjm(t)s(t) expCjw 0 t3j- (1) 

where S is the average signal power, m(t) is a unit power digital data 
modulation, Wq is the carrier frequency in rad/sec, and s(t) is the com- 
plex spreading waveform, namely, 

s(t) = PN^t) + j PN 2 (t) k (2) 

Here we consider the case where PN-j (t) = PN(t) and PN^ft) = PN-j (t - A/2) , 
where A is the PN chip interval. Assuming for simplicity that m(t) is 
a real modulation, e.g , BPSK, combining (1) and (2) then gives 



2 


x(t) *» JS m(t)|pN(t) cos » 0 t- PN(t-A/Z) sin u 0 tj (3) 

Furthermore, since the data rate of m(t) is typically much lower than the 
PN chip rate (large processing gain), we may ignore m(t) in our power 
spectrum calculation. Thus, we are interested in the PSD of the normalized 
signal 


y(t) = PN(t) cos to Q t - PN(t - A/2) sin a> 0 t 


(4) 


The autocorrelation function of y{t) is easily seen to be 


Ry(-r) = <y(t) y(t+r)> 


" ^os w q t 2 


R pN (t+A/2) - R dw (t-A/2) 


PN' 


sin UqT (5) 


where < > denotes time average, the overbar denotes statistical average, 
and Rp^(r) is the autocorrelation function of the PN sequence. Taking 
the Fourier transform of (5), we get 


if” 

Sy(w) * 2 S PN^“" W (P + S PN^-“ + “(p 


+ \ y* R pN (t+ A/2) sin « 0 t e -Ja,T dx 

— co 

- \ J R pN (t-A/2) sin w q t e"' ]WT dt 


( 6 ) 


Letting x = x + a/2 or x- a/2 in the second and third terms of (6), the 
sum of these two terms becomes 
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h / F - e-j^“,>)(x- A / 2 J] dx 

-00 

- h f ?S pN (x)[e- jC “-“(P (x+i/2) - dx 

-CO 

= F S PN^“"“Cp S 1 n {^“ _t0 Cp A>r ' F S PN^ + “(P sin [C“ +w o) A//2 


(7) 


Finally, using (7) in (6) gives the desired result, namely. 



The corresponding result for independent PN sequences PN-j(t) and P^Ct) 
would be simply 



(9) 


Note that (8) and (9) are valid independent of the shape of the PN pulse. 

As an example, consider the common case of a rectangular PN 
pulse. For long codes, we may then describe the PSD of the sequence by 
its envelope, namely, 



S pN (w) = A 


( 10 ) 





Substituting (10) into (8) yields 


ORIGIN^ 

QF POOR Q !J 


s y U) = f 


n (to - ujq}a/2 
Go - Uq3 A/2 


-,n 2 


1 + si 


n jC“ - (Oq^a/^J 


4- A 

2 - 


in|((>>+ coq3a/2 

&> + “(p A/2 


-f \ 2 


1 - SI 


n C“ +(t) {p A ^ 


( 11 ) 


Figure 1 is an illustration of S y («) as described by (11). For simplicity, 
we show the result for positive o> only with the understanding that S y (w) 
is an even function of w Also shown in dotted lines is the corresponding 
result for independent PN sequences obtained by substituting (10) into (9). 




